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A rapid phase of accelerated expansion in the early universe, known as inflation, dilutes all matter 
except the vacuum induced quantum fluctuations. These are responsible for seeding the initial 
perturbations in the baryonic matter, the non-baryonic dark matter and the observed tempera- 
ture anisotropy in the cosmic microwave background (CMB) radiation. To explain the universe 
observed today, the end of inflation must also excite a thermal bath filled with baryons, an amount 
of baryon asymmetry, and dark matter. We review the current understanding of inflation, dark 
matter, mechanisms for generating matter-anti-matter asymmetry, and the prospects for testing 
them at ground and space based experiments. 
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I. INTRODUCTION 



This review aims at building a consistent picture of the 
early universe where the three pillars of modern cosmol- 
ogy: inflation, baryogenesis and the synthesis of dark 
matter can be understood in a testable framework of 
physics beyo nd the Stand ard Model (SM). 

Inflation (lGuthl . [T98l . which is a rapid phase of ac- 
celerated expansion of space, is the leading model that 
explains the origin of matter; during this phase, primor- 
dial density perturbations are also str etched from sub- 
Hubble to super- Hubble length scales ( Mukhanov et al] . 
,1992 ). A strong support for such an inflationary sce- 
nario comes from the precision measurement of these per- 
turbations in the cosmic microwave background (CMB) 
radiation, e.g. by th e Cosm ic Background Explorer 
(COBE) (ISmoot et oA Il992h and the Wilkinson Mi- 
crow ave Anisotropy Probe (WMAP) (jKomatsu et al\ . 
I2OIII ) satellites. However, one of the most serious 
challenges faced by inflationary models is that only a 
few of them provide clear predictions for crucial ques- 
tions regarding the nature of the matter created after 
inflation and the mode of exiting inflation in a vac- 
uum that can excite the SM degrees of freedom (d.o.f) 
( Mazumdar and R,ochej . [201lll . 

From observations we know that the current universe 
contains 4.6% atoms, 23% non-relativistic, non-luminous 
dark matter, and the rest in the form of dark energy. 
While some 13.7 billion years ago it was 37% atoms. 
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photons and neutrinos, and 63% non-relativistic dark 
matter ( Komatsu et al\ . l201ll ). Therefore, it is manda- 



tory that the inflationary vacuum must excite these SM 
baryons, and create the right abundance of dark mat- 
ter. Si nce the success o f Big Bang Nucleosynthesis 
(BBN) (jlocco et al\ . l2009[ ) requires an asymmetry be- 
tween the baryons and anti-baryons of order one part in 
10^*^, it is necessary that the baryonic asymmetry must 
have been cre ated dynamicany in the early universe be- 
fore the BBN (|SakharovL 119671) . 

The prime question is what sort of visible sector be- 
yond the SM would accomplish all these goals - infla- 
tion, matter creation, and seed perturbations for the 
CMB. Beyond the scale of electroweak SM (at energies 
above > 10 — 1000 GeV) th e re are plethora of candi- 
dates, e.g. ( Bustamante et al\ . l2009t) . However the low 
scale supersymmetry (SUSY) provides an excellent plat- 
form, which have been built on the success of the elec- 
troweak phys i cs (Chung et al\. 2005t iHaber and Kand . 
Il985t iMartinl . llQQTt iNilled . Il984 . The minimal super- 
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symmetric extension of the SM, known as MSSM, or 
its minimal extensions, pro vides many testabl e imprints 
at the collider experiment ( Nath et al\ . |2010() . In par- 
ticular, the lightest SUSY particle (LSP) can be elec- 
trically neutral, and will be an ideal candidate for the 
weakly interacting massive particle (WIM P) as a dark 
matter (|Ellis et a/.l . ll984al:lGoldberid.ll983l ). whose abun- 
dance can now be calculated from the direct decay of the 
inflaton, or from the decay products of the inflaton, as 
shown in Fig. [TJ 

If such a visible sector, with the known gauge inter- 
actions, can also provide us with an inflationary poten- 
tial capable of matching the current CMB data, then 
we would be able to identify the origin of the inflaton, 
its mass and couplings, and the vacuum energy den- 
sity within a testable theory, such as the MSSM. The 
inflaton's gauge invariant couplings would enable us to 
ascertain the post-inflationary dynamics, and the exact 
mechanism for particle creation from th e inflaton's coher- 
ent osc illations, known as (p)reheating ( Allahverdi et al\ . 
l2010aD . We would be able to precisely determine the 
largest reheat temperature, Tr, of the post-inflationary 
universe, during which all the MSSM d.o.f come in chemi- 
cal and in kinetic equilibrium for the first time ever. Once 
the relevant d.o.f are created it would be possible to build 
a coherent picture where we will be able to understand 
the origin of baryogenesis and the dark matter in a con- 
sistent framework as illustrated in Fig. [T] 

This review is divided into three parts. In the first 
part we will discuss the origin of infiation, and how to 
connect the models of infiation to the current CMB ob- 
servations. We will keep our discussions general and pro- 
vide some examples of non-SUSY models of inflation. We 
will mainly focus on SUSY based models and its gener- 
alization to supergravity (SUGRA). We wih discuss the 
epoch of reheating, preheating and thermalization for an 
MSSM based models of inflation. In the second part 
of the review, we will focus on baryogenesis. We will 
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FIG. 1: An illustration of a visible sector model for the 
early universe. EW stands for the electroweak. 



state the conditions for generating baryogenesis. We will 
discuss electroweak baryogenesis, baryogenesis induced 
by lepton asymmetry, known as leptogenesis, and the 
MSSM based Afflck-Dine baryogenesis which can create 
non-topological solitons, known as Q-balls. In the third 
part, we will consider general properties of dark matter, 
various mechanisms for creating them, some well moti- 
vated candidates, and link the origin of dark matter to 
the origin of inflation within SUSY. We will briefly dis- 
cuss the ongoing searches of WIMP as a dark matter 
candidate. 



II. PARTICLE PHYSICS ORIGIN OF INFLATION 

There are two classes of models of inflation, which 
have bee n discussed extensively in the literature, e.g. see 
revie ws ( Lvth and Riottol . Il999l : [Mazumdar and Rocheil . 
l201lh . In the first class, the inflaton field belongs to the 
hidden sector (not charged under the SM gauge group). 
The direction along which the infiaton field rolls belongs 
to an absolute gauge singlet, whose couplings to the visi- 
ble sector - such as that of the SM or MSSM fields are not 
determined a-priori. A singlet infiaton would couple to 
the visible and hidden sectors without any biased - such 
as the case of gravity which is a true singlet, and a color 
and flavor blind. In the second class, the inflaton can- 
didate distinctly belongs to the visible sector, where the 
inflaton is charged under the SM or its minimal extension 
beyond the SM gauge group. This has many advantages, 
which we will discuss in some details. 

Any inflationary models are required to be tested by 
the amplitude of the density perturbations for t he ob - 
served large scale structures ('Mukhano v et all Il992h . 
Therefore the predictions for the CMB fluctuations are 
the most important ones to judge the merits of the mod- 
els, which would contain information about the power 
spectrum, the tilt in the spectrum, running in the tilt, 
and tensor to scalar ratio. These observable quantities 
can be recast in terms of the properties of the poten- 
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tial which we will discuss below. From the particle ori- 
gin point of view, one of the successful criteria is to end 
inflation in the right vacuum - where the SM baryons 
are ex cited naturally for a successfu l baryogenesis before 
BBN f|Mazumdar and RocheilbOllI) . 



A. Properties of inflation 

The inflaton direction which leads to a graceful exit 
needs to be flat with a non-negligible slope provided 
by a potential V{(j)) which dominates the energy den- 
sity of the universe. A completely flat potential, or 
a false vacuum with a very tiny tunneling rate to a 
lower vac uum, would r ender inflatio n future eternal, but 



not p ast (iBorde et all 2 003: Bor de and Vilenkinl . Il994l: 
iLindd . Il983lll986l: iLinde et all Il994ll996l) . A past eter- 



nal inflati on is possible only if the n ull geodesies are past 
complete (iBiswas 011 . 120101 [20061) . The slow-roll infla- 
tion assumes that the potential dominates over the ki- 
netic energy of the inflaton 0^ ^ and (p V'{(j)), 
therefore the Friedmann and the Klein-Gordon equations 
are approximated as: 



ij2 V{(l))/3Mj 
« -V'{^), 



p ' 



(1) 
(2) 



where prime denotes derivative with respect to 0. There 
exists slow-roll conditions, which constrain the shape of 
the potential, are give by: 



^{V'/Vf«l, 



\i^{4>)\ = M2|y"/y|«l. 



(3) 
(4) 



These conditions are necessary but not sufficient for in- 
flation. The slow- roll conditions are violated when e '-^ 1, 
and rj 1, which marks the end of inflation. 

However, there are certain models where this 
need no t be t rue, f or instance in hybrid inflation 
models (iLindd . I1994D . where inflation comes to an 
end via a phase transition, or in oscillatory mod- 
els of inflation where slow-roll conditions are satis- 
fied only on average (iDainour and Mukhanovl Il998l : 
iLiddle and Mazumdaii 1998), or in f lation happens in 
oscillations ("Bisw as and Mazumdad . l2009l) . or in fast 
roll i nflation whe re the slow-roll conditions are never 
met dLindd . l200ll) . The K-inflation where only the ki- 
net ic term dominates where ther e is no potential at 
all ( Armendariz-Picon et aLl . ll999t ). 

One of the salient features of the slow-roll inflation 
is that there exists a late time attractor behavior, such 
that the evolution of a scalar fleld after sufficient e- 
foldings become indepen dent of the initial conditions 
(ISalopek and Bondl . Il990l) . 

The number of e-foldings between, i, and the end of 
inflation, tend, is defined by: 



iV^ln^%^ 
a{t) 



Hdt 



1 

14 



^d(b, (5) 



where (pend is defined by e{(pend) 1, provided infla- 
tion comes to an end via a violation of the slow-roll 
conditions. The number of e-foldings can be related 
to the Hubble crossing mode k — akHk by comparing 
with the present Hubble length uoHq. The final result is 
(|Liddle and Leachl . 120031 ) 



N{k) = 62 - In 
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lO^^GeV 
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where the subscripts end (R) refer to the end of inflation 
(end of reheating). Today's Hubble length would cor- 
respond to Nq = N{k — aoHo) number of e-foldings, 
whose actual value would depend on the equation of 
state, i.e. w = p/p (p denotes the pressure, p denotes 
the energy density), from the end of inflation to radi- 
ation and matter dominated epochs. A high scale in- 
flation with a prompt reheating with relativistic species 
would yield approximately, Nq « 50 — 60. A significant 



low (lArkani-Hamed et al. 


^200C 


; Green and Mazumdar 


2OO2I: Lvth and Stewart!. 


19951 


19961: iMazumdar. 199S 


Mazumdar and Perez-Lorenzanal 200ll). 



B. Density Perturbations 

1. Scalar perturbations 

Small inhomogeneities in the scalar field, (t>{x, t) = 
(j){t) + S(j){x, t), such that S(f> <^ </>, induce perturbations in 
the background metric, but the separation between the 
background metric and a perturbed one is not unique. 
One needs to choose a gauge. A simple choice would be to 
fix a gauge where the non-relativistic limit of the full per- 
turbed Einstein equation can be recast as a Poisson equa- 
tion with a Newtonian gravitational potential, The 
induc ed metric can be written as, e.g. ( Mukhanov et all 
I1992D : 

ds^ = a2(r) [(1 -I- 2<i>)dT^ - (1 - 2'^)5ikdx'dx''] , (7) 

Only in the presence of Einstein gravity and when the 
spatial part of the energy momentum tensor is diagonal, 
i.e. ST^ = (5j, it follows that $ = 

During infiation the massless infiaton (with mass 
squared: ^ V <^ H"^) perturbations, (5(/), are 

stretched outside the Hubble patch. One can track their 
perturbations from a sub-Hubble to that of a super- 
Hubble length scales. Right at the time when the wave 
numbers are crossing the Hubble patch, one finds a solu- 
tion for 56 as 



{\8^k?) = {H{t.,fl2k^)., 



(8) 



where denotes the instance of Hubble crossing. One 
can define a power spectrum for the perturbed scalar field 



27r2 



\H{t.)] 
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' H' 


2 


2tt 




2^ 





(9) 
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Note that the phase of 54>k can be arbitrary, and there- 
fore, inflation has generated a Gaussian perturbation. 
Now, one has to calculate the power spectrum for the 
metric perturbations. For a critical density universe 



5p_ 
P 



k 



(10) 



where ^^(i) « {'i / {5(l)k / (j))\k=aH ■ Therefore, one ob- 
tains: 

where the right hand side can be evaluated at the time of 
horizon exit k = aH . The temperature anisotropy seen 
by the observer in the matter dominated epoch is propor- 
tional to the Newtonian potential, ATk/T = -(l/3)$fe. 

Besides tracking the perturbations in the longitudinal 
gauge with the help of Newtonian potential, there ex- 
ists another useful gauge known as the comoving gauge. 
By definition, this choice of gauge requires a comov- 
ing hypersurface on which the energy flux vanishes, 
and the relevant perturbation amplitud e is known as 
the comoving cu r vature perturbation, (^k ( Lukashl . 1 19801 
iMukhanov et all . Il992f ). For the super-Hubble modes, 
/c — ^ 0, the comoving curvature perturbation, is a 
conserved quantity, and it is proportional to the New- 
tonian potential, (k — — (5/3)$fe. Therefore, Sk can 
also be expressed in terms o f curvature perturbations 
(|Liddle and Lvthl . ll993l[2000h 



2 / k 
5 Va^ 



Ck , 



(12) 



SI 



and the corresponding power spectrum 
{4/25)Vc{k) = (4/25)(H/(^)2(ff/27r)2. With the 

help of the slow-roll equation 3H<j) = —V, and the 
critical density formula 3H^Mp ~ V, one obtains 



Si 



V 



757r2M| y2 



1507r2M4 e 



Vdk) 



1 



V 



247r2M^ e 



(13) 



where we have used the slow-roll parameter e = 
{M^/2){V'/V)^. The COBE satellite measured the 
CMB anisotropy and fixes the normalization of 'P(;{k) on 
very large scales. If we assume that the primordial spec- 
trum can be approximated by a power law (ignoring the 
gravitatio nal waves and th e fc— dependence of the power 
Us) (Komatsu et all |2009[ ) 



k 

Vcik) ~ (2.445 ± 0.096) x 10"^ ( — 

ko 



(14) 



where ng is called the spectral index (or spectral tilt), 
the reference scale is: fco = 7.5aoi?o ^ 0.002 Mpc~^, and 
the error bar on the normalization is given at la, and 



n,{ko) = 0.960 ±0.13 



(15) 



It is important to stress that these central values and 
error bars vary significantly when other parameters are 
introduced to fit the data, in part because of degenera- 
cies between parameters (in particular with Ofjft.^^ tJ^g 
optical depth r, its running, the tensor-to-scalar ratio, r, 
and the fraction of cosmic strings). The spectral index 
n(fc) is defined as 



i{k) 



dink 



(16) 



This definition is equivalent to the power law behavior 
if n{k) is close to a constant quantity over a range of 
k of interest. One particular value of interest is Hs = 
n{ko). If Us = 1, the spect rum is flat and known as 
Harr ison-Zeldovich spectrum ( Harrisonl Il970t [Zeldovichl . 
|1970D . For Us ^ I, the spectrum is tilted, and > 1 
(ris < 1) is known as a blue (red) spectrum. In the slow- 
roll approximation, this tilt can be expressed in terms of 
the slow-roll parameters and at first order: 



-6e + 2f] + Oie^,v\er,,e) , 



(17) 



where 



Ml 



V'{d^V/d(l)^) 



y'2(dV/d(/)4 



(18) 

The running of t h ese parameters are given 
by ( Salopek and Bondl . Il990l ). Since the slow-roll 
inflation requires that e <C 1, |r;| <C 1, therefore naturally 
predicts small vari ation in the spectr al index within 
A In fc « 1 ( Kosows kv and Turned . 119951) 



dn(fc) 
dlnfc 



-1667; + 24e2 + 2^2 



(19) 



It is possible to extend the calculation of metric pertur- 
bation beyond the slow- roll approximations based on a 
formalism similar to that developed in Refs . ( Kolb et all . 
119951: iMukhanovl I1985L Il989t ISasakiL I1986D . 



2. Multi-field perturbations 

Inflation can proceed along many flat directions with 
many light fields. Their perturbations can be tracked 
conveniently in a comoving gauge, on large scales = 
—HS(j)/(f> remains a good conserved quantity, provided 
each field follow slow-roll condition. The comoving cur- 
vature perturbations can be related to the number of e- 
foldings, N, given by ( Salopekl . ll995l : [Sasaki and Steward . 
I1996D 



C = 5N= {dN/dc^a)Ha , 



(20) 



where N is measured by a comoving observer while pass- 
ing from fiat hypersurface (which defines 5(j)) to the 
comoving hypersurface (which determines C). The re- 
peated indices are summed over and the subscript a 
denotes a component of the inflaton ( Lvth and Liddlel . 
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120091: iLvth and Riottd . [l999 ') . If the random fluctuations 
along S4>a have an amplitude {H/2tt)'^, one obtains: 



V 



dN dN 



757r2 Af 2 



(21) 



For a single component dN/dcf) = [M^'^V/V), and then 
Eq. (PT|) reduces to Eq. ([T^ . By using slow-roll equations 
we can again define the spectral index 



M^VqVa 



, M^N,aN,bVab 



(22) 

where = dV/d4)a, and similarly TVq = dN/d<j)a- For 
a single component we recover Eq. p7)) from Eq. (P^ . 
In the case of multi-fields, one has to distinguish adi- 
abatic from isocurvature perturbations. Present CMB 
data r ules out pur e isocurvature perturb a tion spec- 
trum (|Beltran et all 12004 iKomatsu efd\ . l2009l) . al- 
though a mixture of adiabatic and isocurvature pertur- 
bations remains a possibility. 



3. Gravitational waves 

During inflation stochastic gravitational waves are ex- 
pected to be prod uced similar to the scalar perturbations 
Ullcn, 1988; Grishchuy, 119751: iGrishchuk and Sidorovl 
[l989; Salmi, 1990). For reviews on gravitat ional waves, 
see (Maggiore. ,2000HMukhanov et ay.l . Il992ll . The grav- 
itational wave perturbations are described by a line ele- 
ment ds^ + 5ds^ , where 

ds^ = a2(T)(dT2 - dxMxi) , Sds^ = — a'^{T) hi jdx^dx^ . 

(23) 

The gauge invariant and conformally invariant 3-tensor 
hij is symmetric, traceless = 0, and divergenceless 

'Vihij = (Vi is a covariant derivative). Massless spin 2 
gravitons have two transverse degrees of freedom (d.o.f) 
For the Einstein gravity, the gravitational wave equa- 
tion of m otion follows that of a massless Klein Gordon 
equation ( Mukhanov et aLl . ll992l) . Especially, for a fiat 



universe 



h] + iHh) + (A:Va2) h)^Q, 



(24) 



As any massless field, the gravitational waves also feel 
the quantum fiuctuations in an expanding background. 
The spectrum mimics that of Eq. ([9|) 



''grav('s) — ^_^2 ( o_ ) 



'P v27r) 



(25) 



k=aH 



Note that the spectrum has a Planck mass suppression, 
which suggests that the amplitude of the gravitational 
waves is smaller compared to that of the scalar perturba- 
tions. Therefore it is usually assumed that their contribu- 
tion to the CMB anisotropy is small. The corresponding 



spectral index can be expanded in terms of the slow-roll 
parameters at first order as 



grav , „ 



nt = 



dln^grav(fc) 

dlnfc 



-2e, 



(26) 



Note that the tensor spectral index is negative. It is 
expected that PLANCK could detect gravity waves if 
r ^ 0.1, however the spectral index will be hard to mea- 
sure in forthcoming experiments. The primordial gravity 
waves can be generated for large field value inflationary 
models. Using the definition of the nu mber of e-foldings 
it is p ossible to d erive the range of Ad (iHotchkiss etall . 
120081 : iLvthl . [l997t iLvth and Liddld . l2009lv i 



16e = r < 0.003 (50/7V)^ (A(/)/Mp) 



(27) 



Note that it is possible to get sizable, r, for 
A(j) <C Mp in assisted inflation (discussed be- 
low), and in inflection point inflation discussed in 
Ref. (jBen-Davan and BrusteiiJ . |2010|) . If the tensor- 
to-scalar ratio r and/or a running are introduced, 
the best fit for Us and error bars (at la) = 
1.017^^:1^43 , as = -0.028 ± 0.020 n, = 0.970 ± 



1.089 



-1-0.070 



r < 



0.015 , r < 0.22 (at 2cr), n, 
0.55 (at 2a) , = -0.053 ± 0.028 (|Komatsu et all 
l2009l) . These data therefore suggest that a red spectrum 
is favored (rig — 1 excluded at 2.5cr from WMAP and 
at 3.1(7 when other data sets are included) if there is no 
running. 



C. Generic models of inflation 

1. High scale models of inflation 

The most general form for the potential of a gauge 
singlet scalar field 4> contains an infinite number of terms, 



V = Vo + 



Xq. 



OL—2. 



(28) 



The renormalizable terms allows to prevent all terms with 
a > A. By imposing the parity Z2, under which — >■ — 0, 
allows to prevent all terms with a odd. Most phenomeno- 
logical models of infiation proposed initially assume that 
one or two terms in Eq. (|28)) dominate over the others, 
though some do contain an infinite number of terms. 

a. Power-law chaotic inflation: The simplest inflation 
model by the nu mber o f free parameters is perhaps the 
chaotic infiation (jLindd . Il983l ) with the potential domi- 
nated by only one of the terms in the above series 



V 



(29) 



with a a positive integer. The first two slow-roll param- 
eters are given by 



•q — a{a — 1) 



Ml 

A2 



(30) 
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Inflation ends when e = 1, reached for (pe = aMp/^/2. 
The largest cosmological scale becomes super-Hubble 
when (pQ — ^j2NQaMp, which is super Planckian; this 
is the first challenge for this class of models. The spectral 
index for the scalar and tensor to scalar ratio read: 



n., = 1 



2 + a 



4q; 



2Nq + a/2 



Nq + a/4 



(31) 



The amplitude of the density perturbations, if normal- 
ized at the COBE scale, yields to extremely small cou- 
pling constants; Aq ^ 1 (for i.e. A4 ~ 3.7 x 10^^^). 
The smallness of the coupling, Aa/Afp"**, is often con- 
sidered as an unnatural fine-tuning. Even when dimen- 
sion full, for example if a = 2, the generation (and 
the stability) of a mass scale, ^/J^Mp ~ 10^'^ GeV, 
is a challenge in theories beyond the SM, as they re- 
quire unnatural cancellations. These class of mod- 
els have an interesting behavior for initial conditions 
with a large phase space distribution where there ex- 
ists a late attractor trajectory leading to an end of in- 
flation when the slow-roll conditions ar e violate d close 
to the Planck scale (,Brandenberger and Kunel Il990t 
iKofman et all . l2002HLindel . Il983lll985[) ! 

Note that the above mentioned monomial potential 
can be a good approximation to describe in a certain 
field range for various models of inflation proposed and 
motivated from particle physics; natural inflati on when 
the i nflaton is a pseudo-Goldstone boson (Freese et al\ . 
I199C1I ). or the Landau-Gin zburg potential when the in- 
flaton is a Higgs-type field ( Bezrukov and Shaposhnikovl 
I2OO8). The necessity of super Planckian VEVs represents 
though a challenge to such embedding in particle physics 
and supergravity (SUGRA). 

b. Exponential potential: An exponential potential 
also belongs to the large field models: 



V(c^) = exp 



pAfp 



(32) 



It would give rise to a power law expansion a{t) cx i^, so 
that inflation occurs when p > \. The case p — 2 cor- 
responds to the exactly de Sitter evolution and a never 
ending accelerated expansion. Even for p ^ 2, viola- 
tion of slow- roll never takes place, since = 1/p and 
inflation has to be ended by a ph ase transition or grav- 
itational production of particles ( Copeland et aLl . [2OOII : 
iLvth and Riottd . 119991 ). 

The confrontation to the CMB data yields: Ug = 1 — 
2/p and r — 16/p; the model predicts a hight tensor to 
scalar ratio and it is within the one sigma contour-plot 
of WMAP (with non-negligible r) for p e [73 - 133]. 



2. Assisted inflation 

Many heavy fields could collectively assist infiation by 
increasing the eff ective Hubbl e frict ion term for all the 
individual fields ( Liddle et al. I I1998D . This idea can be 



illustrated with the help of 'm! identical scalar fields with 
an exponential potentials, see Eq. p2p . where now (j) — >■ 
where i = 1,2,--- ,m. For a particular solution; 
where all the scalar fields are equal: 0i = 02 = • • • = 4'm- 

1 



Mil 



m[y(0i) + 0f/2] 



-'iH^i-dV{cj)i)/d(i)i 



(33) 
(34) 



These can be mapped to the equations of a model 
with a single scalar field (j) by the redefinitions (t>\ — 
m(j)\ ; V — mV ; p — mp, so the expansion rate 
is a (X iP, provided that p > 1/3. The expansion becomes 
quicker the more scalar fields there are. In particular, po- 
tentials with p < 1, which for a single field are unable to 
support inflation, can do so as long as there are enough 
scalar flelds to make mp > 1. 

In order to calculate the density perturbation produced 
in multi-scalar field models, we recall the results from 
Eq. jUl). Since N ^ - f H dt, we have J2^ §^^t = ~H, 

we yield: Vc = {H/2Tr)^{l/m){H^ /4>l). 

Note that this last expression only contains one of 
the scalar fields, chosen arbitrarily to be (f>i. The 
estimation for the spectral tilt is given by : n — 
1 = —2/mp, which matches that produced by a sin- 
gle scalar field with p — mp. The more scalar fields 
there are, the closer to scale-invariance is the spec- 
trum that they produce. The above calculation can 
be repeated for arbitrary slopes, pi. In which case the 
spectral tilt would have been given by n = 1 — 2/p, 
where p — "^Pi- The above scenario has been gen- 
eralized to study arbitrary exponential potential s with 
coup l ings, V = y^"' Zgexpiy^ "^ asi(j)i) (jCopeland et all 
119991 : iGreen and Lidsevl . l2000l) . 

a. Assisted chaotic inflation: Multi-scalar 

field s of chaotic type has interesting proper- 
ties ( Jokinen and Mazumdail , |2004t ): 



(35) 



(for 



> 4). 



driven at VEVs, 



The chaotic inflation can now be 
< M p, below the Planck 



scale ( Kanti and Olivd . Il999al lbl). The effective slow- 
roll parameters are given by: eeff = e/m <^ 1 and 
I'^e/zl = 1^1/"^ <?C 1, where e, 77 are the slow-roll parame- 
ters for the individual flelds. Inflation can now occur for 



flcld VEVs (| Jokinen and MazumdaA |2004|) : 



Afp \ m / \ 60 / V 



600 



Nc 
60 



1/2 



< 1 



(36) 



where Nq is the number of e-foldings. Obviously, all the 
properties of chaotic inflation can be preserved at VEVs 
^ Afp, including the prediction for the tensor to scalar 
ration for the stochastic gravity waves, i.e. r = 16ee//. 
For eeff ~ 0.01 and m ^ 100, it is possible to realize 
a sub-Planckian inflation, the spectral tilt close to the 
flatness: Ug — 1 = —Qeeff + 2rieff, and large tensor to 
scalar ratio, i.e. r = 0.16. 
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b. N-flation: Amongst various realizations of assisted 
inflation, N-flation is perhaps the most interesting one. 
The idea is to have N - 300(Afp//) - 10"^ number of 
axions, where / is the axion decay constant, of order 
/ ^ O.lAfp ^ dri ve inflation simultaneousl y with a leading 
order potential ( Dimopoulos et ad 120051 ): 



F = yo+^Atcos(0,//,) + 



(37) 



where (f>i are axion flelds correspond to the partners of 
Kahler moduli. The ellipses contain higher order con- 
tributions. In a certain Type-IIB compactification, it is 
assumed that all the moduli are heavy and thus stabi- 
lized by prior dynamics, including that of the volume 
modu lus. Only the axions of Tj — 4>i/fi + iM^R"^ are 
light ()Dimopoulos et all |2005[ ). The assumption of de- 
coupling the dynamics of Kah ler modulus fro m the axions 
is still a debatable issue, see (iKalloshlliooi . After rear- 
ranging the potential for the axions, and expanding them 
around their minima for a canonical choice of the kinetic 
terms, the Lagrangian simplifles to the lowest order in 
expansion: 



(38) 



The exact calculation of rrii is hard, assuming all of 
the mass terms to be the same m.i ~ lO^'^ GeV, and 
> (Mp//)^, it is possible to match the current ob- 
servations with a tilt in the spectrum, n ^ 0.97, and 
large tensor to scalar ratio: r ^ 8/Nq ^ 0.13 for 
Nq ~ 60. T here are also r ealiza t ions of assisted inflation 
via branes (iBecker et al\ . 120051 : ICline and Stoical . l2005t 
iMazumdar et all l200lh 



3. Hybrid inflation 

The end of inflation can happen via a waterfall trig- 
gered by a Higgs (not necessarily the SM Higgs) field cou- 
pled t o the i nflaton, first discussed in ( Copeland et all 
11994 iLindd . Il99ll ll994l) The mod el is based on the 
potential given by (|Lindel . Il99ll 11994 



where </> is the inflaton and is the Higgs-type field. A 
and A' are two positive coupling constants, m and M 
are two mass parameters. It is the most general form 
(omitting a quartic term X" of renormalizable poten- 
tial satisfying the symmetries: ip ^ —ip and cj) —cf). 
Inflation takes place along the ip = valley and ends 
with a tachyonic instability for the Higgs-type held. The 
critical point of instability occurs at: 



0c = 



(40) 



The system then evolves toward its true minimum at V 
0, (0) = 0, and (ip) = ±M. 



The inflationary valley, for (■0) = 0, where the last 50 — 
60 e-foldings of inflation is assumed to take. This raises 
the issue of initial conditions for (</>, il ^) fields and the fine 
tuning required to initiate inflation (IClesse and Rocher . 
20081: iLazarides and Vlachosl. fl997l: iMendes and Liddle . 



20001 : iPanagiotakopoulos and Tetradisl. 119991: iTetradis . 
1998h . In Ref. (|Clesse and Rocheil . l2008( ) it was found 
that when the initial VEV of the inflaton, (p <C Mp, a 
subdominant but non-negligible part of the initial condi- 
tions for the phase space leads to a successful inflation, 
i.e. around less than 15% depending on the model pa- 
rameters. Initial conditions with super-Planckian VEVs 
for (f) ^ Mp automatically leads to a successful inflation 
similarly to chaotic inflation. In the inflationary valley, 
(tp) — 0, the effective potential is given by: 



K=ff(0) 



AM" 



1 



2 J,2 



-m <p' 



(41) 



The model predicts a blue tilt in the spectrum, i.e. Us > 
1, in the small field regime, (pQ < Mp, which is slightly 
disfavored by the current data. 

Two variations of the hybrid inflation idea were pro- 
posed assuming that the term (p"^ is negligible. The two- 
field scalar potentials are of the form: 



ArV 



(42) 



They share the common feature of having an infla- 
tionary trajectory during which {ip) is varying and 
not vanishing. For (p, q) = (1, 2), t he mo del is 
known as Mutated hybrid ( Stewartl . Il995b[ ). and 
(jo, q) = (4, 6) corresponds to Smooth h ybrid inflation 
( Lazarides and Panagiotakopoulosl Il995( ). The latter in- 



volves non-renormalizable terms of order Mp ^ to keep 
the potential bounded from below. 

The potential is valid in the large field limit cp ^ cp^,, 
since in the small field limit, the potential is not bounded 
from below and should be completed. For mutated, the 
model predicts a red spectral index and negligible ten- 
sor to scalar ratio, — 1 ~ —3/{8Nq) ~ 0.97, and 



3m/(2AAr^/2) 



3/(8A^^) 



10- 



if we assume 



Nq ~ 60. For smooth, the end of slow- roll infiation hap- 
pens by a violation of the conditions; e,r] <^ 1, since no 
waterfall transition takes place. This allows the predic- 
tions for the spectral index to be ns — 1 —5 /(3No) ~ 
0.97 ( Lazarides and Panagiotakopoulosl Il995l ). and the 
ratio for tensor to scalar is found to be negligible. 



4. Inflection point inflation 

One of the challenges for inflation is to realize inflation 
at low scales, preferably below Mp, with the right tilt and 
the amplitude of the power spectrum. Inflection point 
inflation admits a large amount of flexibility in the held 
space - similar to the analogy of a ball rolling on an 
elastic surface following the least action principle. With 
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the help of two independent parameters, A and B, it is 
possible to obtain a large range of tilt in the spectrum, 
while keeping the amplitude of the perturbations intact. 
Let us consider a simple realization of such a potential: 



where C — f{A, B) in order to obtain a point of inflec- 
tion suitable for inflation. The VEV at which inflation 
occurs is intimately related to the two independent pa- 
rameters and can happen at wide ranging scales below 
Mp, and for wide ranging values of (A, B). 

Here we will generalize this potential to any generic 
potential V which can be written in the following 
form (here / den otes differentiati on with re spect to 
(j)) (lEngvist et a/] . i2010a; .Hotchkiss eTaLl . l201l[) : 



V = Vo + a( 



)+ 2(</'-0o)' + g(</'-0o)' + -(44) 



where = l/(0o), a = F'((/)o), h = F"(0o), c = 
^"'("/"o), which is the Taylor expansion, truncated at 
n = 3, around a reference point which we choose to be 
the point of inflection where F"((/)o) = 0, or saddle point 
where F"((/)o) — V^'((/>o) = 0. The higher order terms in 
Eq. (|44|) can be neglected during inflation, provided that 



I0e - (/-ol 



m > 4 . 



(45) 



where </)e corresponds to the field value at the end of 
inflation. Assuming that the slow-roll parameters are 
small in the vicinity of the inflection point and that 
the velocity is negligible, the potential energy Vq gives 
rise to a period of inflation. 

Inflation ends at the point (/)e where I77I ^ 1. By solv- 
ing the equation of motion, the number of e-foldings of 
inflation during the slow-roll motion of the inflaton from 
to 0e, whe re (^n - (<?!'o ~ fe) < ^ < 0o + (0o - 0e), is 
found to be (|Enavist e< all l2010a^ 



^ = ^\/-[^o(0e)-f^o(0)] 

Ain V ac 



Fo(z) = arccot W — (z-(/)o) 
la 



(46) 



aM] 



^ and Y 



It useful to define the parameters X — —rrry 

V 2 Vo 

^J^N M-pX . Note that X is the square root of the slow- 
roll parameter e at the point of inflection. The slow-roll 
parameters can then be recast in the following form: 



2Vi 



Y 



-^|(^/r-Xcosr- 



Xsiny 



2 

iV2 



(47) 
(48) 
(49) 



where S = y/1 — X sin Y + Va cosF. One can solve 
Eqs. (|47M9j) . for X, Y and N in terms of the slow- 
roll. The equations are non-linear and in general can- 
not be solved analytically. However, since e ^ \r]\, ^, 
one can flnd a closed form solution provided that 



(43) y^'^ < IQie GeV and X < ^/i <r 1 dAllahverdi et~ 



2007c-~'Bucno SanchezeFoH 120071: lEnqvist et all l2010a ; 
Hotchkiss et aL . ,20111 ): 



V, 



1/2 



1 



1/2 







Vr 



1/2 







sin^r, (50) 

4 

l + 2r]-6e=l FcotF, (51) 



724^ ei/2M2 ~ 2ttV6M^X 
4 



Y 



svaY 



(52) 



One can derive the properties of a saddle point inflation 
provided Y/svaY 1, and Y coiY — >■ 1. The model 
favors the current observations by matching the COBE 
normalization and the spectral tilt ranging from Ug e 
[0.93, 1.0]. For instance, the lowest value corresponds to 
the saddle point inflation for Nn ~ 60. 



D. Supersymmetric models 

One of the most compelling virtues of SUSY is that 
it can protect the quadratically divergent contributions 
to the scalar mass, which arise in one-loop computation 
from the fermion contribution and quartic self interac- 
tion of the scalar field. Such corrections generically spoil 
the fiatness of the inflaton potential. The quadratic di- 
vergence is independent of the mass of the scalar field 
and cancel, exactly if \g = Aj, where A/ is the fermion 
Yukawa and As is the quartic scalar coupling. However 
this procedure fails at 2-loops and one requires fine tun- 
ing of the couplings order by order in perturbation the- 
ory. In the case of the SM Higgs, a precision of roughly 
one part in 10^^ is required in couplings to maintain the 
Higgs potential, often known as the hierarchy problem or 
the naturalness problem. The electroweak symmetry is 
still broken by the Higgs mechanism, but the quadratic 
divergences in the scalar sector are absent. In the SUSY 
limit the fermion and scalar masses are degenerate, but 
the SUSY has to be broken soflly at the TeV scale in such 



a way that it docs not spo i l the solut ion to the hier archy 
problem, see (IChung et aLl . 120051: iHaber and Ka'ia . ll985l: 



iMartinl . 119971 : iNillesl . Il98" 

The matter fields for = 1 SUSY are chiral super- 
fields = -I- \/29'ip + 96F, which describe a scalar (/), 
a fermion -ip and a scalar auxiliary field F. The SUSY 
scalar potential V is the sum of the F- and _D-terms: 



dW 



(53) 
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where W is the superpotential, and 4'i transforms under 
a gauge group G with the generators of the Lie algebra 
given by T°. Note that ah the kinetic energy terms are 
included in the D-terms. For inflation, the effects of su- 
pergravity (SUGRA) becomes important. At tree level, 
= 1 SUG RA potential is given by the sum of F and 
D-terms, see (|Nillesl . Il984 



V 



w 



Ml 



(54) 
(55) 



where we have added the Fayet-Iliopoulos contribution 
f to the D-term, and Z)" = D°- /g"" , where g"" is gauge 
coupling. Here K{(j)i, (j)*^) is the Kahler potential, which 
is a function of the fields 4'i: a-nd = dK/d(j)i. In 
the simplest case, at tree-level K = 4'*^(j)i (and = 
{K^^)j — Sf). In general the Kahler potential can be 
expanded as: K = + [k"^ (i),(j)jcj)*^ + c.c.)/Mp + 

(/c;>,(/.,(/.*V*V*V*' + A;;^'V»<^j0fc0*' + c.c.)/M2 + ...). 
The kinetic terms for the scalars take the form: 



(56) 



The real part of the gauge kinetic function matrix is given 

by Re fab- In general, fab = 6ab{i/gl + />j/Mp H ). 

The gauginos masses are typically given by nixa = 
Re[/^](^i)/2Mp. For a universal gaugino masses, /* are 
the same for all the three gauge groups of MSSM. In the 
simplest case, it is just a constant, fab — Sab/9ai a-^d the 
kinetic terms for the gauge potentials, AJ^, are given by: 



\{Refab)F^.Fr ■ 



(57) 



SUGRA will be broken if one or more of the Fi obtain a 
VEV. The gravitino, spin ±3/2 component of the gravi- 
ton, then absorb the Goldstino component to become 
massive. Requiring classically {V) = 0, as a constraint 
to obtain the zero cosmological constant, one obtains 



'■3/2 



{K]F^F* 
3M2 



\{W)? 
Ml. 



(58) 



1. F-term inflation 



The most well-known model of SUSY inflation driven 
by i^-terms is of the hybrid type and based on t he su- 



perpotential (ICopeland et all Il994l: iDvali et al\ Il994l: 
iLinde and Riottol . Il997l ) 



W = kS'($$ - M^) 



(59) 



where, S is an absolute gauge singlet, while $ and $ are 
two distinct superfields belonging to complex conjugate 



representation, and k is an arbitrary constant fixed by the 
CMB observations. It is desirable to obtain an effective 
singlet S supe rfield arising from a higher gauge theory 
such as GUT (Langackei |l98l|), 

however to our knowl- 
edge it has not be en possible to implement this i dea, see 
the discussion in ( Mazumdar and Rocheil l201l[ ). Typ- 
ically S would have other (self)couplings which would 
effectively ruin the flatness required for hybrid inflation. 

This form of potential is protected from additional 
destabilizing contributions with higher power of S, if S, 
$ and $ carrying respectively the charges +2, a and —a 
under R-parity. Then W carries a charge +2 so that the 
action S — J (P9 W + . . . is invariant. 

The tree level scalar potential derived from Eq. ((59|) 
reads 



Utrcc(5, 0, (/.) = k2|M2 



(60) 



where we have denoted by S, 0, 4> the scalar components 
of 5",$,$. Note the similarity between Eq. and 
Eq. ([Bn]). where m — 0, and both A and A' are replaced by 
only . We will also assume (j)* — <j> along this direction, 
and the kinetic terms for the superfields are minimal, i.e. 
with a kahler potential: K = \S\'^ + |$p + 

Let us define two effective real scalar fields canonically 
normalized, a = V2Re{S), and ip = 2Re(<I>) = 2Re($), 
the overall potential can then be recast as: 

Utrce(a, ^) ^ (^M' - t) ' + T^''^' ■ ^^^^ 

The global minimum of the potential is located at S* = 0, 
(f^ = aP. At large VEVs, S > Sc = M, the potential 
also possesses a local valley of minima (at (■0) — 0) in 
which the field a, now rolls on with Utrcc = Vq = k^AI'^. 
This non-vanishing value of the potential both sustain 
the infiationary dynamics and induces a SUSY breaking. 

This induces a splitting in the mass of the fermionic 
and bosonic components of the superfields $ and $, 
with m|(5) = K^l^p ± k'^M^ and mj, = k^\S\^. 
Note that this description is valid only as long as S 
is sufficiently slow- rolling such that K^jSpl^p can be 
considered as a mass term. There fore radiatiye cor- 
rections do not exactly cancel out (jDvali et all Il994l : 
iLazarideslboOOD , and provide a one-loop potential which 
can be calcula ted by using the Col eman- Weinberg for- 
mula (Coleman and Weinberd . fl973() . 



l^l-loop(0) = V,nfi4)+^V 



,(62) 



where Vinf is now the renormalized potential, A(0) is 
the renormalization mass scale. The sum extends over 
all helicity states i, Fi is the fermion number, and M((j)) 
is the mass of the i-th state. One obtains: 



— loop iS) 



K^UM^ 

327r2 



21n^ + (z + l)2ln(I + z-i) 



-f(z-l)2ln(I-z-i)] m 
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where z = |S'p/M^ = cc^, A represents a non-physical 
energy scale of renormalization and N denotes the di- 
mensionality. Note that the perturbative approach of 
Coleman and Weinberg breaks down when close to the 
inflection point at z ~ 1. For small coupling k, the slow- 
roll conditions (for 77) are violated infinitely close to the 
critical point, z = 1, which ends inflation. 

The normalization to COBE allows to fix the scale M 
as a function of k. If the breaking of G does not pro- 
duce cosmic strings, the contribution to the quadrupole 
anisotropy simply comes from the inflationary contribu- 
tion (see Eq. (HH)) and the observed value can be ob- 
taine d even with a coupling k close to unity ( Dvali et al\ . 
I1994D . Small values of k can render the scale of infla- 



tion very low, as low as the T eV scale (Bastcr o-Gil et al 
20031: iBastero-Gil and Kinp . [T998; .Randall et all 11996 : 



Randall and Thomasl . ll995li 

However it has been shown that the formation 
of cosmic strings at the end of _F-term inflation is 
highly probable when the mo del is embedded in SUSY 
GUTs jjeannerot eFoZI . |2003[ ) . In this case, the normal- 
ization to COBE receives two contrib utions, one from 
inflation and other from cosmic st rings ( JeanneroRll997l : 
iRocher and Sakellariadoull2005al ). which affects the rela- 
tion M{k) at large k, and i mposes new stringent bounds 
on M < 2 X 10^^ GeV, and (iJeannerot and Postmall2005l: 
iRocher and Sakellariadoul l2005b( ) 



< 7 X 10-^(126/7Vq) 



(64) 



by demanding that the cosmic strings cam at best 
contri bute less than < 10% of isocurvature fluctua- 
tions ( Bevis et all l2008l ). Once M is fixed, the spec- 
tral index rig can be computed as the range is found 
to be: n.s € [0.98, 1] whether cosmic strings form 
or not ( Jeannerot and Postmall2005lSenoguz and Shafil . 
I2OO3I ). and by including the soft-SUSY breaking terms 
within minimal kinetic terms in the Kahler potential, the 
spectr al index can be brou ght down to 0.928 < n^. < 
1.008 (|Rehman et a/.l . 120091 ). 



2. SUGRA corrections and solutions 

For inflaton VEVs below the Planck scale, the SUGRA 
effects can become important and may ruin the flatness of 
the potential. The iV = 1 SUGRA potential is now given 
by Eq. (j54p . where the _F-terms containing an additional 
exponential factor. Various cross terms between the 



mass term for the light scalar flelds (iBertolami and Ross 


1987; ICoDeland et all 1994; 


Dine et all 1984 1995b 


1996bl: Linde and Riotto. 1997 


) 



2 _^ V 

''^SUGRA ^susy + r^jyp 



0{1)H^ , (65) 



Once the inflaton gets a mass ~ H , the contribution to 
the second slow-roll parameter 77 becomes order unity 
and the slow roll inflation ends, i.e. \t]\ = M'^V" /V ^ 
"^sugra/H^ - 0(1). This is known as the SUGRA-r; 
problem. 

When there are more than one chiral superfields, as 
in the F-term hybrid model, it can be possible to cancel 
the dominant 0{1)H correction to the inflaton mass by 
choosing an appropri ate Kahler term (jCopeland et all 
11994 IStewaTtl . Il995a^ ■ For non-minimal Kahler poten- 
tials, such as 



K=\S\^ + \^\^ + \<^>\^ + Ks\S\yA4 



(66) 



where rususy ~ 0(100) GeV contains soft-SUSY breaking 
mass term for the low scale SUSY breaking scenarios. 



the kinetic terms if are non-minimal because 

K^i ^ 5^3. One obtains: {dss'K)~^ - 1-Aks\S\'^/MI + 
. . . One again obtains a problematic contribution to the 
inflaton mass, i.e. k,s^O{1)H. Several mechanisms have 
been proposed to tackle this ry-problem. One can impose, 
Ks ~ 10"'^, which is sufficient to keep the slow roll infla- 
tion safe, but without much physical motivation. For a 
generic inflationary model it is not possible to compute 
at all. 

a. Shift and Heisenberg symmetry: Safe non- 
minimal Kahler pot e ntials h ave also been pro- 
poscd (Antusch et al., '2009b'; 'Bastero -Gil and Kind . 
1999; Brax and Martin, 2005; PaUis, 200i) making use 
of the shift symmetry. Under this symm etry, a superfield 
S — » S + i C, where C is a constant. ( Kawasaki et all 
I2OOOI . l200l[ )_ to protect the Kahler potential of the 
form K{^,^) — ^ + $). This symmetry generates 
an exactly flat direction for an inflaton fleld and a 
non-invariance of the superpotential induces some slope 
to its potential to allow slow-roll at the loop level. 
Another symmetry - the Heisenberg symmetry - has 
also been invoked to pr o tect th e form of the Kahler 
potential ( Antusch et all l2009a[) . where the effective 
Kahler is a no-scale type of the form K = ln($i). 
This solves the SUGRA-77-problem by canceling the 
exponential term exp{K). However note that Kahler 
potentials generically obtains quantum corrections 
unlike the non-renormalization th eorem which c an only 
protect the superpotential terms ( Grisaru et all [197^). 
Such corrections are hard to compute without knowing 
the ultra-violet completion, and the exact mat ter sector 
for the inflationary model (|Berg et all l2005al lbl. [2006.) . 
Note that none of these papers considered MSSM matter 
sector. 

b. Inflection point inflation: For any smooth poten- 
tial, it is possible to drive inflation near the saddle 
point, V = 0, V" = 0,U"' ^ 0, or near the point 
of inflection, V ^ 0, U" = 0, U'" ^ 0. These 
are special points where the effective mass term of 
the inflaton vanishes and the potential does not suf- 



fer through SUGRA-77 proble m (jAllahverdi et all l2007d. 
I2OO6I : iMazumdar et all l201l[ ) . In the saddle point case 
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FIG. 2: Regions of parameter space (green) for the 
potential in Eq. dST]) that satisfy the WMAP 7-year 
constraints on ((5jy, ris). Contour lines of |/3cmb| are 
shown in black, for the values of |/3cmb| indicated for 
n = 6 case. For 0(1)7? co rrection to the inflaton mass 
the value of /3 - 10"^ (jMazumdar tt al\ . l201ll ). 



the potential ca n be made so flat that inflatio n can be 
driven eternally ([AUahverdi et al\ . l2006l l2007d[ ). 

From the low-energy perspective, the most generic and 
dangerous SUGRA corrections to the inflaton potential 
( with minimal and non-minimal Kahler potentials for 
(/> ) would have a large vacuum energy contribution. To 
complicate further, one may even assume that the flat- 
ness of (/) is lifted by n on-renormaliza ble contribution to 
the potential (Mazumdar et al\ . l201ll) : 



V(</>) - K + 



101 



2("-l) 



2(n-3) 



M 



(67) 



where Vc — 3H^Mp. As mentioned above the interesting 
observation is that, in fact, there always exists a range 
of field values, A(f>, for which a full potential admits a 
point of inflection with all known sources of corrections 
taken into account. Now, all the uncertainties in the 
corrections to the Kahler potential can be absorbed in the 
full potential, such that th e flat r egion admits a slow roll 
inflation with At/) » (/)o (jMazumdar et all i201ll) . The 



condition for this inflection point is ajj 



where we characterize the flne-tuning by (3 defined as: 



1)CH 



= 1 - 



(68) 



When 1/3 1 is small, a point of inflection (/)o exists such that 
V" ((/)o) = 0, with 



2{n- 1) 



(69) 



We can Taylor expand the potential about (pQ as discussed 
in section III. C. 41 and analyze the CMB constraints as 
shown in Fig. [51 



3. D-term inflation 



In Re fs. (iBinetruv and Dvalil . Il996l: iHalvol Il996t 
[steward . Il995a ) , it was noticed that a perfectly flat in- 



flaton potentials can be constructed using a constant 
contribution coming from the Z?-term. In addition, the 
SUGRA- ?7 problem arising in i^-term models does not ap- 
pear for ZJ-terms driven inflation because the £)-sector of 
the potential does not receive exponential contributions 
from non-minimal SUGRA. The model however requires 
the presence of a Fayet-Iliopoulos (FI) term ^, and there- 
fore a 1/(1)^ symmetry which generates it. For a Kahler 
potential _ft'($m, $„), the ZJ-terms 



-ga[Da = MTaTjK^ 



(where = dK/d^„i) give rise to a scalar potential: 
V^(</),r) = i[Re/((/.)]-i ^i^'^A,. + F - terms (70) 

where ga and T" are respectively the gauge coupling con- 
stants and the generators of each factors of the symmetry 
of the action, 'a' running over all factors of the symmetry, 
and /((/>) is the gauge kinetic function. If this symme- 
try contains a factor f7(l)^, the most general action then 
allows for the presence of a constant contribution ^. 

The simplest realization of Z?-term inflation repro- 
duces the hybrid potential with three chiral superfields, 
5, (/)+, and </)_ with non-anomalous U{\)^ (an abelian 
theory is said to be anomalous if the trace of the 
generator is non- vanishing Y^.Qn. # 0) charges g„ = 
0,-^1,-1 (iBinetruv and Dvali [l99a iHalvol . Il996[) . The 
superpotential can be written as 



W 



D 



(71) 



In what follows, we assume the minimal structure for 
/($i) (i.e., f{^i)—l) and take the minimal Kahler po- 
tential, i.e. K — |(/)-P + 
potential reads 



D- SUGRA 



+i(/._5n 1+ 



= A'^ exp 
. , 1^1' 



1 5 p.. Then the scalar 



Ml 



M4 



Af4 



10- 



M4 



Ml 



(72) 



I (10.1^ 



where is the gauge coupling of U{1)^. The global 
minimum of the potential is obtained for (S) = and 
(<&_) = which is SUSY preserving but induces the 
breaking of U{1)^. For S > Sinst = g^V^/^ the potential 
is minimized for \(f>+\ — 10-1 — and therefore, at the 
tree level, the potential exhibits a fiat infiationary valley, 
with vacuum energy Vq = 

The radiative corrections depend on the splitting 
between the effective masses of the components of 
the superfields $_|_ and because of the tran- 

sient Z?-term SUSY breaking. The radiative cor- 
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recti ons are Kiven by the Colema n- Weinberg expres- 
sion ( Coleman and Weinberg) . Il973() and the full poten- 
tial inside the inflationary valley reads 



ic2 



D-SUGRA 
ff 



167r2 



21n^^exp 



\sl 



(z -I- 1)2 ln(l + z-^) + {z- if ln(l - z~i) 



(73) 



with z = (A^l^p/^l^el'S'l'/^'^p. Inflation ends when the 
slow- roll conditions break down, that is for z^nd — Ij a-nd 
the predictions for the inflationary parameters are very 
similar to the previous discussion on i^-term inflation. 



E. MSSM flat direction inflation 

So far we have discussed inflationary models where the 
inflaton sector belongs to the hidden sector (not charged 
under the SM gauge group), such models will have at 
least one SM gauge singlet component, whose couplings 
to other fields and mass are chosen just to match the 
CMB observations. These models are simple but lack 
proper embedding within MSSM or its minimal exten- 
sions. 

In order to construct a predictable hidden sector model 
of inflation, one must know all the inflaton couplings to 
the hidden and visible matter. One such unique model 
has been constructed within type IIB string theory, where 
it was found that all the inflato n energy is transferred 
to exc iting the hidden matter ( Cicoli and Mazumdail . 
I2010al lbl). and the universe could be prematurely dom- 
inated by the hidden sector dark matter. Such obstacles 
do not arise if the last phase of inflation occurs within 
MSSM. 



1. Introducing MSSM and its flat directions 

In addition to the usual quark and lepton superfields, 
MSSM has two Higgs fields, Hu and Hd- Two Hig- 
gses are needed because W is forbidden in the super- 
potential^ The supcrpotcnti al for the MS SM is g iven by , 



see llChung et all 120051: iHaber and Kanel . ll985l:lMart 



I1984D 



Wmssm = KQHuU + XdQHdd + XeLHdC + nHuHd , 

(74) 

where Hu, Hd,Q, L^u,d,e in Eq. (j74p are chiral super- 
fields, and the dimensionless Yukawa couplings A^,, A^, Ag 
are 3x3 matrices in the family space. We have sup- 
pressed the gauge and family indices. The Hu, Hd,Q,L 
fields are SU{2) doublets, while u, d, e are SU{2) sin- 
glets. The last term is the fj, term, which is a SUSY ver- 
sion of the SM Higgs boson mass. Terms proportional to 
H*Hu or H^Hd are forbidden in the superpotential, since 



Wmssm must be analytic in the chiral fields. Hu and Hd 
are required not only because they give masses to all the 
quarks and leptons, but also for the cancellation of gauge 
anomalies. The Yukawa matrices determine the masses 
and CKM mixing angles of the ordinary quarks and lep- 
tons through the neutral components of Hu — [Hu , H^) 
and Hd = {H^H^). Since the top quark, bottom quark 
and tau lepton are the heaviest fermions in the SM, we 
assume that only the third family, (3,3) element of the 
matrices A„,Ad,Ae are important. 

The II term provides masses to the Higgsinos 



CD-i^iH+H^ -H^H'd) + c.c, 



(75) 



and contributes to the Higgs (massf terms in the scalar 
potential through 



CdVd \ii\^i\H 



0|2 



\H 



+ |2 . 



\H'd\' 



\Hd\') 



(76) 



Note that Eq. ([7^ is positive definite. Therefore, it can- 
not lead to electroweak symmetry breaking without in- 
cluding SUSY breaking {massf soft terms for the Higgs 
fields, which can be negative. Hence, |/Ltp should almost 
cancel the negative soft {massf term in order to allow 
for a Higgs VEV of order - 174 GeV. That the two 
different sources of masses should be precisely of same 
order is a puzzle for which many s olutions has been sug- 
geste d (ICasas a nd Munoj. Il993l : ICiudice and Masierol 
119881 : iKim and NiUea . 1981 . 

Within MSSM one can construct gauge invariant D- 
and f -fl at direction s , for the list of MSSM f iat di rec- 
tions see (|Dine adll996bl : lGherghetta et aLl . ll996l ). A 
flat direction can be represented by a composite gauge 
invariant operator, Xm, formed from the product of 
k chiral superfields $i making up the flat direction: 
Xm = ^1^2 ■■ ^m- The scalar component of the su- 
perfield is related to the o rder parameter (p through 
Xm = c(j)"^ (iDine et aLl . [l996bD . 

A n example of a D- a nd f - f lat direction is provided 
by ( Dine and Kusenkcl |2004 lEngvist and Mazumdail , 
120031) 



''''-Til' 



(77) 



where (/> is a complex field parameterizing the flat di- 
rection, or the order parameter, or the AD field. All 
the other fields are set to zero. In terms of the com- 
posite gauge invariant operators, we would write X^ — 
LHu (to = 2). Note that a flat direction necessarily car- 
ries a global U{\) quantum number, which corresponds 
to an invariance of the effective Lagrangian for the or- 



der parameter cf) under phase rotation 



I. In the 



MSSM the global U{1) symmetry is B — L. For example, 
the Lff„-direction has B — L = — 1 . 

From Eq. ([80)l one clearly obtains F^^ = XuQu + 
fiHd = Fl — XdHdC = for all (/). However there ex- 
ists a non-zero F-component given by F^^ = [iHu- Since 
/i can not be much larger than the electroweak scale 
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M\Y ^ C(l) TeV, this contribution is of the same order 
as the soft SUSY breaking masses, w hich are going to hf t 
the degeneracy. Therefore, following ( Dine et aZl . ll996bl ). 
one may nevertheless consider LHu to correspond to a F- 
flat direction. The relevant I?-terms read 

^51/(2) = Hlr.H, + L^T.L = i|<^p - i|<^|2 ^0. (78) 

Therefore the LH^ direction is also _D-flat. 



2. Gauge invariant inflatons of IVISSM 

A simple observa ti on w as fir st made in 
(jAllahverdi et all l2007d. l2006l l2007e[ ). where the 
inflaton properties are directly related to the soft SUSY 
breaking mass term and the A-tcrm of the MSSM. 
Within MSSM, it is possible to lift the flatness of the 
gauge invariant combinations of squarks and sleptons 
away from the point of enhanced gauge symmetry by the 
i^-term, while maintaining the Z?-flatness. 



a. Squarks and sleptons driven inflation: Let 

us c onsider a non-renormalizable superpotential 
term ( Dine and Kusenkol 12004 lEngvist and Mazumdail . 
I2003D : 



E 

n>3 



n M" 



(79) 



Where ^ — <j> exp[i6'], while 6 is the phase term is a 
gauge invariant superfield which contains the flat direc- 
tion. Within MSSM (with conserved i?-parity) all the 
flat directions are lifted by the non-renornializa ble op- 
erators with 4 < n < 9 (|Gherghetta et all Il996i ). Two 
distinct directions are: udd and LLe, up to an overall 
phase factor they are parameterized by: 



(80) 



1 



where 1 < a, /?, 7 < 3 (a 7^ /3 7^ 7) are color indices, and 
^ l£ h ji k < S (j k) denote the quark families for udd, 
and \ < a,h < 2 [a ^ h) are the weak isospin indices 
and 1 < i, j, k < 3 {i ^ j ^ k) denote the lepton families 
for LLe. Both these dire ctions are lifted by n = 6 non- 
renormalizable operators (jGherghetta et al\ . \l99(i) . 



We D jj^{LLe){LLe) 



M3 



{udd)[udd). (82) 



Rest of the directions within MSSM are lifted by hybrid 
operators of type, W ^ (1/Mp~^)^'$"^^, which does not 
lead to c osmoloqically flat potential via ble fo r slow-roll 
inflation (jAllahverdi et all 2008, 2006, 2007eD . 

The scalar potential along these directions includes 
softly broken SUSY mass term for </> and an A-term gives 



rise to a specific potential ( AUahverdi et aLl . l2008l [2OO6I 
l2007eD 



F(0) 



.\4>?~A 



6Mf 



A- 



101 



10 



(83) 



The A-term is a positive quantity with dimension of 
mass. Note that the first and third terms in Eq. ((83)) 
are positive definite, while the A-term leads to a nega- 
tive contribution along the directions whenever cos(n^? -|- 
9a) < 0. The above potential is similar to Eg. di^ . It 
is possible to find a point of inflection, (/)o, provided that 
A^/40m'^ =1-1- 4a^, where ^ 1, and at the lowest 
orders in 0{a'^), we obtain: 

ViM = ^ml4,l + • • • , V'[c^o) = 4a2m20o + • • • , 

l/"(0o) = O, V"'{cj)o)=i2ml/cj)o + --- . 

(j)o = (m^Ml/ATTo) ^ * + 0{a^) . (84) 

In the case of gravity-mediated SUSY breaking scenarios, 
~ A - m3/2 - (100 GeV - 1 TeV). Therefore the 
condition ^ 40771^ can indeed be satisfied. Inflation 
occurs within an interval 10 — (/)o| ~ 0o/6OA/p <C Mp, in 
the vicinity of the point of inflection, 0o ~ 0(10^^ GeV). 
Within which the slow-roll parameters, e,rj <^ 1. The 
Hubble expansion rate during inflation is given by 



MSSM 



1 m^4>o 
"H Mp 



(100 MeV- 1 GeV). (85) 



The amplitude of density perturbations 5h (see 
Eqs. (IT3I I50I51I) a nd the scalar s pectra l index 
Us are given by (lAllahverdi et all l2006l l2007e : 



AUahverdi and Mazumdar . 2006al: Bueno Sanchez et al. . 
20071) : 



5h = 



8 m^Mp 1 



Us = 1-4%/A2 cot[iVQ%/A2], 



sin^tA^QV^] (86) 
(87) 



where 2 x 10"^ < = 9OOa27Vg^(Mp/0o)^ < 5.2 x 
10^^, and Nq ~ 50. Running in the tilt is very small. In 
this case the universe thermalizes in to MSSM radiation 
instantly in less than one Hubble time after the end of 
inflation ((Allahverdi et al. I. I2011H) . see the discussion in 
Sect. lH.F.3l 

For 00 ^ 10^^ GeV, there is an apparent fine-tuning in 
the parameters A/m^ = a ~ 10"^*^, which may look un- 
pleasant. However note that this fine tuning between the 
two MSSM parameters in the ratio is energy dependent 
and valid only at the scale of inflation at 10^* GeV, but 
at the TeV scale where the soft masses would be mea- 
sured at the co Uider there is no apparen t fine tuning in 
the parameters ( AUahverdi et all l2010ct) . 

As shown in Sect. III.D.2.b[ see Fig. H the SUGRA 
corrections will ameliorate the tuning down to /? = 
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FIG. 3: {6h, ris) is plotted for different values of 
and A = 1. The 2a region for Sh is shown by the blue 
horizontal band and the 2a allowed region of Hs is 
shown by the vertical green band. The la allowed 
region of i s within the solid ve rtical 
lines (|Allahverdi et aLl . l2008l ). 



a ^ 10^^, virtually addressing any fine tuning re- 
quired for the success of MSSM inflation. It was 
shown in Refs. (jAllahverdi et al ]. l2008l . l2007dl) . that the 
inflection point for the MSSM inflaton is an attrac- 
tor solution, provided there exists a phase of inflation 
prior of that of the MSSM with N > 10^° e-foldings. 
Such large e-foldi ngs can be gene r ated w ithin string 
theory landscape (lAUahverdi aTl . l2007dl ). or within 
MSSM (lAllahverdi et a^.l . l2008^ . 



b. Renormalizable superpotential: The left handed 
neutrinos can be of Dirac type with an appropriate 
Yukawa coupling. The simplest way to obtain this would 
be to augment the SM symmetry by, SU{3)c x SU(2) l x 
U{1)y X U{1)b-l, where U{1)b^l is gauged. The rele- 
vant superpotential term is 



W D hNHuL. 



(88) 



Here A'', L and are superflelds containing the RH 
neutrinos, left-handed (LH) leptons and the Higgs which 
gives mass to the up-type quarks, respectively. Note that 
the NHuL monomial represents a D-flat direction under 
the U{1)b-l, as weU as the SM gauge group. 

The value of h needs to be small, i.e. h < 10~^^, in 
order to explain the light neutrino mass, ^ 0(0.1 eV) 
corresponding to the atmospheric neutrino oscillations 
detected by Super-Kamiokande experiment. The poten- 
tial along this direction, after the minimizatio n along 
the angular direction, is found to be (jAllahverdi et all 
I2007by i. 



Ah 



(89) 
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V{(j>o) 7^ 0, V"{(f>Q) = 0, where inflation takes place 

'0.05 eVN 



V3- 



6 X 10^2 mw, 

V m,, / 



^00) = ^ ^3 



X 10^4 ml 



/0.05 eV\2 
\ m„ J 



(90) 



The amplitude of density perturbation s follows from 
Eqs. (Il3l50l5ll) (lAllahverdi et aLl . l2007bB 'l . 



1 f 

— ^ ~ 3.5 X 10-27 
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0.05 cV 



7710 



(91) 

Here denotes the loop-corrected value of the infla- 
ton mass at the scale 0o in Eqs. (|90|91|) . The spectral 
tilt as usual has a ran ge of values 0.90 < n., < 1.0 
(lAllahverdi et all 120071 



c. MSSM Higgses as inflaton: The MSSM Higgses are 
another flne example of a vis ible sector inflaton provided 
some restrictions are met (jChatteriee and Mazumdail . 
l201lf ). The required superpotential is given by 



Afc {Hu-Hd) 



2fc-3 ' 



(92) 



This is the ^-term which were considered an 
ideal candidate to genera te the density perturba- 
tions (|Enavist et all l2004allbh . but now they can also 
provide the required vacuum energy to inflate the 
universe ( Chatteriee and Mazumdail l201ll) . The scalar 
potential along the HuHd I?-flat direction is given by, 

V{^,e) = im2(^?)^2^(-l)(^-i)2AfcAicos((2ft-2)0)^2fc 

+ 2A;V'^''"'\ (93) 

where ^ = V^l^je*", and = (l/%/2)((/), 0)^, Hd = 
(l/V2)-i(0,</))^, and 



ni^ie) = -{m%^^ + mjj^ + 2^^ -23(1 cos2e), (94) 

(95) 



For A w 4m0, there exists an inflection point for which 



For simplicity, we may assume fi and B to be real. 
This choice is compatible with the experimental con- 
straints, ma inly from the Elect r on Di pole Moment mea- 
surements ( Pospelov and Ri"t3 , l2005l ). The inflection 
point can be obtained for 9 = 0, ±7r/2, for simplicity 
let us consider the case for 6 = 0, and when the fol- 
lowing condition is satisfied, ttiq = k'^fi'^/{2k — 1) -I- A^, 
where A is the tuning required to maintain the flatness 
of the potential. Although, this tuning could be harsh 
at the inflationary scale, if ^ 10^^ GeV, but the ra- 
tio evolves to m^/ ^"^ ^ C'(l) at the electroweak scale 
by vir tue of running of the renormaliza tion group equa- 
tions (I Chatteriee and Mazumdail l201ll) . The amplitude 
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of the CMB perturbations can be obtained for A2 ~ 10 ® 
and A3 ~ 0{1), it is possible to obtain a similar plot like 
Fig. dll) for Higgs mass m{e = 0) - 100-250 GeV, which 
yields the spectral tilt in the range 0.93 < Ug < 0.98. 



F. Preheating, reheating, thermalization 

Reheating at the end of inflation is an important 
aspect of inflationary cosmology. Without reheating 
the unive rse would be ernpty o f matter, for a re- 
view see (|Allahverdi et all l2010al) . Reheating occurs 
through coupling of the inflaton field (p, to the SM 
matter. Such couplings must be present at least via 
gravitational interactions. In particular, if the infla- 
ton is a SM gauge singlet, the relevant couplings to 



1. Non-perturbative particle creation 



SM are: ^cj>{Hqi)qR , ^c^F^uF^''' , g^4>''HH, where 
M is the scale below which all these effective opera- 
tors are valid, X, g ^ C'(l), H is the SM Higgs dou- 
blet, and q;,(7R are the left and the right handed SM 
fermions ( Allahverdi and Mazumdail l2007b() . 



Similar couplings would arise if is replaced by right 
handed sneutrinos, axions, moduli, or any other hidden 
sector field. Being a SM singlet, can as well couple 
to other hidden sectors, moduli, axions, etc. Since the 
hidden sectors are largely unknown, it becomes a chal- 
lenge for a singlet inflaton to deca y solely into the SM 
do.f (jCicoh and Mazumdail l2010alfb[ ). 



After the end of inflation, the inflaton starts coher- 
ent oscillations around its minimum. The frequency 
of oscillations are determined by the frequency of os- 
cillations, uj ~ iTT-eff ^ Hinf. Dur ing this epoch 
the inflaton can decay perturbativelv (lAlbrecht et~al. , 
19821: iDolgov and Kirilovalll990tlKorb and Turned. ll98S : 



Turner, 19831) . Averaging over many oscillations re- 
sults in a pressureless equation of state where {p) = 
((/)^/2 — V{(j))) vanishes, so that the energy density starts 
evolving like a matter domination (in a quadratic poten- 
tial) with = pi{ai/a)^ (subscript i denotes the quanti- 
ties right after the end of inflation). For potential the 
coherent oscillations yield an effective equation of state 
similar to that of a radiation epoch. If represents the 
total decay width of the inflaton to pairs of fermions. This 
releases the energy into the therma l bath of relativistic 
particles when H{a) = yf{l/Mll)pi{ai/af''^ « F^,. The 
energy density of the thermal bath is determined by the 
reheat temperature T^, given by: 



(96) 

where denotes the effective relativistic d.o.f in the 
plasma. However the inflaton decay products need to 
thermalize, which requires acquiring kinetic and chemi- 
cal equilibrium. 



If the inflaton coupling to the matter field is large, 
a completely new channel of reheating opens up due 
to t he coherent na t ure o f the i n flaton field, pr o posed 
bv (iKofman et all Il994 Il997l: fShtanov et all (l995t 
ITraschen and Brandenbergeii 199ol ). known as vreheat- 
ing. Let us first consider a simple toy model with inter- 
action Lagrangian 



1 



Ant = ~^g X 



2,,2 /2 



(97) 



where x is another scalar field, in a realistic set-up x 
could take the role of the SM Higgs. We can neglect 
the effect of expansion provided that the time period of 
preheating is small compared to the Hubble expansion 
time H^^, this is reasonable in many cases. 

The quantum theory of x particle production in the 
external classical infiaton background begins by expand- 
ing the quantum field x into creation and annihilation 
operators and aj^ as: 

X(i,x) = j^l^ J d'k {xUt)ake''''' + Xk{t)aie-^''-) , 

(98) 

where k is the momentum. Since the equation of motion 
for X is linear it can be studied simply mode by mode in 
Fourier space. The mode functions then satisfy: 

Xk + (fc' + ml+ gH^ sin^ {m^t)) Xk ^ , (99) 

where $ is the amplitude of oscillation in (p. This is the 
Mathieu equation which is written in the form 



Xk + {Ak-2qcos2z)xk = 0, 



(100) 



where the dimensionless time variable is z = m^t and 
a prime now denotes the derivative with respect to z. 
Comparing the coefficients, we find 



A, = 



2q 



q = 



9^ 



(101) 



The growth of the mode function corr esponds to particle 
production ( Birrell and Davie^ |1982| ). It is well known 
that the above Mathieu equation Eq. (|100p has instabil- 
ities for certain ranges of k: 



Xk oc exp(^fc2;) , 



(102) 



where Hk is called the Floquet exponent. For small val- 
ues of (7, i.e. g <C 1, resonance occurs in a narrow in- 
stability band about k = m, known as a "narrow reso- 
nance" band ( Traschen and Brandenberge r. 'I99tf) . The 
resonance is much more efficient if q S> 1 (IKofman et all 
I1994] . I1997D . In this case, resonance occurs in broad 
bands, i.e. the bands include all long wavelength modes 
/c — !■ 0, known as broad resonance. This can be un- 
derstood by studying the condition for particle produc- 
tion in the WKB approximation for the evolution of x 
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field which is violated. In the WKB approximation: 
which is valid as long as the adiabaticity 

condition 



dujl/dt < 2ujI 



(103) 



is satisfied. In the above, the effective frequency ujk is 
given by 



(104) 



By inserting the effective frequency Eq. (jl04l) into the 
condition Eq. (I103p and following some algebra, the adi- 
abaticity condition is violated for momenta 



e < 



3v^ 



gTO0<l> 



(105) 



For modes with these values of k, the adiabaticity condi- 
tion breaks down in each oscillation period when (j) is close 
to zero. The particle nu mber does not increase smo othly, 
but rather in "bursts" ()Kofman et ami994 ITqQTI ). 

The above analysis of neglecting the expansion of the 
uni verse is self-consisten t. However, as discussed in detail 
in ( Kofman et all Il997l ). the expansion of space can be 
included. The equation of motion for x becomes 

Xk + 5Hxk g'^'^itf sir? {m^t)j Xk = 0. 

(106) 

The adiabaticity condition is now violated for momenta 
satisfying: 



(107) 



Note that the expansion of space makes broad resonance 
more effective since more k modes are red-shifted into the 
instability band as time proceeds. The detailed analysis 
yields the same expression for the resonance band except 
for the exact value of the numerical coefficient of the 
first term on the r.h.s.. Broad parametric resonance ends 
when q < 1/4. 

In principle, it is also possible to excite the fermions 
non-perturbatively, in spite of the fact that the occupa- 
tion number of any fixed state cannot be gr eater than one 
(because of the Paul i exclus i on principle) (IBaacke et al. , 
19981: iGiudice et a/]. Il999at ICreene and Kofmanl . Il99£ 



2OO0I: IPeloso and Sorbol. l2000l). and higher spin ±3/2 
gravitinos (iGiudice et a ?. 1999b': 'K allosh et o L. '2000a''b'; 
iMaroto and Mazumdan . ,2000; .Nilles et a^.l . l2001K c'). 

a. Tachyonic prehetaing: It is possible that effective 
frequency of certain mode can be negative. For example 
in a symmetry breaking potential: V{4)) = j\{(tP' — rfY , 
for small field values, the effective mass of the fluctua- 
tions of (/> is negative and hence a "tachyonic " resonance 
will occur, as studied in (jFelder et al ]. l200lh . For small 
field values, the equation for the fluctuations (pk of is 



The modes with k < m grow with an exponent which 
approaches /i^ = 1 in the limit — )■ 0. Given initial 
vacuum amplitudes for the modes (f)k at the intial time 
t = of the resonance, the held dispersion at a later time 
t will be given by 



kdk 
An ' 



2t^/^mfrk^ 



(109) 



The growth of the fluctuations modes terminates once the 
dispersion becomes comparable to the symmetry break- 
ing scale. 

Tachyonic preheating also occurs in hybrid inflation 
models, see Eq. (|39|) . In this case, it is the fluctuations 
of ip w hich have tachyonic form and which grow exponen- 
tially ( Felder et all 12001 ) . Preheating in hybrid infla tion 
was first studied in ( Garcia-Bellido and Lindd . Il998l) us- 
ing the tools of broad parametric resonance. 

b. End of preheating: In the above analysis we have 
neglected the back-reaction of the produced x and 
particles on the dynamics. The presence of x parti- 
cles changes the effective mass of the infiaton oscilla- 
tions. This back-reaction effect is negligible as long as 
the change Am^ in the square mass of the infiaton is 
smaller than m^. In the Hartree approximation, the 
change in the infiat on mass due to y pa rticles is given 
by Am^ = (|Kofman et all Il997l ). Another im- 

portant condition is that the energy in the x particles 
should be sub-dominant. Therefore, p^^ ^ ((^x)^) — 
fc^(x^) ^ m'^{(p'^), It was found that is smaller than 
the potential energy of the infiaton field at the time ti 
as long as the value q at the time ti is larger than 1, 
i.e. q{ti) > 1. This is roughly speaking the same as 
the condition f or th e effectiveness of broad resonance 
(jKofman et a/.l . Il997t) . 



2. Thermalization 

Neither the perturbative decay of the infiaton nor the 
preheating mechanism produce a thermal spectrum of de- 
cay products. In a full thermal equilibrium the energy 
density p and the number density n of relativistic parti- 
cles scale as: p ^ T"*^ and n ^ T^, where T is the tem- 
perature of the thermal bath. Thus, in full equilibrium 
the average particle energy is given by: (-E)cq = {p/n)^ 



which obeys the scaling, {E) 



oq 



,1/4 _ T. 



0. 



(108) 



a. Perturbative reheating and thermalization: If the 

infiaton decays perturbatively, then right after the infla- 
ton has decayed completely, the energy density of the uni- 
verse is given by: p ~ 3 (r^Mp)^ , {E) « 7710 ^ p^/*. 
From conservation of energy, the number density of de- 
cayed particles is: n sa {p/m^) p?!'^. Hence, pertur- 
bative decay results in a dilute plasma that contains a 
small number of very energetic particles. A local thermal 
equilibrium requires re-distribution of the energy among 
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different particles, kinetic equilibrium, as well as increas- 
ing the total number of particles, chemical equilibrium. 
Therefore both number-conserving and number- violating 
reactions must be involved. 

The most important processes for kinetic equili- 
bration are 2 — > 2 scatterings with gauge boson 
exchange in the i-channel. Due to an infrared 
singularity, these s catterings are very efficie n t even 
in a dilute plasma (AUahyerdi and Mazumdail . l2006bl : 
[Davidson and Sarkail . 20001 ) . Chemical equilibrium is 
achieved by changing the number of particles in the re- 
heat plasma. In order to reach full equilibrium the total 
number of particles must increase by a factor of ricq/n. 



where 

«3/4 



i p/m and the equilibrium value is: Ucq 
This can be a very large number, i.e. Uea/n 



0(10 ' ^). It was recognized in ([AUahverdi and Drees 
2002; Davidson and Sarkar . 2000l). se e also ( AUahverdi 



2000 ; iJaikumar and MazumdaJ " 2004) that the most rel 
evant processes are 2 — >■ 3 scatterings with gauge- 
boson exchange in the t— channel. When these scatter- 
ing become efficient, the number of particles increases 
very rapidly, and f ull thermal equilibr i um is established 
shortly after that ()Enqvist and Sirkkal IToOSl) . 

b. Non-perturbative preheating and thermalization: 

In this case the occupation numbers of the excited quanta 
are typically very high after the initial stages of preheat- 
ing, Once the occupation numbers of the resonant modes 
become suffici ently large, re-scattering of_the fluctua- 
tions begins (iKhlebnikov and Tkachevl Il996l Il997allbl : 



iMicha and TkachevI, l2003l 12004) which terminates the 
phase of exponential growth of the occupation num- 
bers. The evolution of the field fluctuations evolves to 
a regime of turbulent scaling which is characteri z ed by 
the spectrum n{k) ^ fc"^/^ (|Micha and Tkachevl . l2003l 
[2004) . which is non-thermal (for a thermal distribution 
we would expect n(fc) ~ k~^). The phase of turbulence 
ends once most of the energy has been drained from the 
inflaton field. At this time quantum processes take over 
and lead to the thermalization of the spectrum. 



3. Calculation of Th within MSSM 

In the case of MSSM infl ation, the inflaton coupl ings 
to MSSM d.o./are known (|Allahverdi et al\ . \20nyi] . It 
is therefore possible to track the thermal history of the 
universe from the end of inflation. When the MSSM 
inflaton passes through minimum, i.e. (/) — 0, the entire 
gauge symmetry gets restored and all the do./ associated 
with the MSSM gauge group become massless, which is 
known as the point of enhanced gauqe symmetry. 

These are the massless modes which couple to the in- 
flaton directly, for instance the d.o.f corresponding to 
SUi2)w X C/(l)y, or that of SU{3)c x U{1)y- At VEVs 
away from the minimum, the same modes become heavy 
and therefore it is kinematically unfavorable to excite 
them. The actual process of excitation depends on how 



strongly the adiabaticity condition for the time depen- 
dent vacuum is violated for the inflaton zero mode. 

a. Couplings for LLe inflaton: Let us illustrate this 
with LiL2ez flat direction as an inflaton. The infla- 
ton non-zero VEV completely breaks the SU{2)w x 
U{1)y symmetry. This results in four massive real 
scalars . whose masses are obtained from the D- 
terms (jAllahverdi et all l2011bl) 



y 3 ^5^0' (x? 



xl) 



\gWxl- (110) 



Here gw, gv are the SU{2)w and U{1)y gauge couplings 
respectively, and denote the inflaton, see Eq. ([5T|) . The 
X's are Goldstone bosons from breakdown of SU{2)w x 
U{1)y. They are eaten by the Higgs mechanism and 
give rise to longitudinal components of the electroweak 
gauge flelds. In the unitary gauge, they are completely 
removed from the spectrum. The x particles decay to 
squarks, the Higgs particles, and the L3, ei, 62 sleptons 
with the decay rates given by: F^^ — T^.^ = T-^^ — 
3 5^(^/87rA/6 , r^, = (9 gf.0/167rV^). Note that the 
decay rate is proportional to the VEV of the inflaton, 
which sets the mass of x fields. Couplings of the infiaton 
to the gauge fields are obtained from the flat direction 
kinetic terms ( AUahverdi et all l2011b[ ) 



Cd^-^<I>'(2W+^^W^ 



(111) 



where W+ = {Wi - iW2)/V2 , W" = {Wi + iW2)/V2, 
where Wi,^ and are the SU{2)w and U{1)y gauge 
flelds respectively. The gauge fields decay to (s)quarks, 
Higgs and Higgsino particles, and L3, ei, 62 (s)leptons 
with the total decay widths: = F^y- = T\y^ = 

(3 g^<j>/8TrV6), Tb = (9 (7|,0/167rV2). Couplings of 
the i nflaton to fermions can also be found in a similar 
way (jAllahverdi et all l2011bl) . 



b. Instant preheating and thermalization: The flelds 
that are coupled to the inflaton acquire a VEV-dependent 
mass that varies in time due to the inflaton oscillations. 
For illustration, we first focus on the xi scalar, which are 
produced every time the infiaton goes through zero. The 
Fourier eigenmodes of xi have the corresponding energy 



(112) 



where is the bare mass of the x field. The 

growth of the occupation number of mode k can be 
computed exactly for the first zero-crossing, Uk^x ~ 

exp —ny/G^k"^ + iri^)/ {gw<}>o) , where the infiaton near 
the zero-crossing is given by 00 = (2y(0))i/2, from the 
conservation of energy, where (p is the amplitude of the 
infiaton oscillations, (j) ~ (j)o/\/3 ^ 10^^ GeV, where (po 
is the infiection point for inflation, Eq. (|84l) . Note that 
after a few oscillations, (f>o ~ rn^(j), since the expansion 
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rate during the inflaton oscillations is negligible by virtue 
of ~ 100 GeV and H{t) < Hi^f - 1 GeV. The total 
number density of particles thus produced follows 



(fik 

W) 



2\/27r3 



exp 



(113) 

where q = ((jf^0Q/24TOp 1. This expression corre- 
sponds to the asymptotic value and assumes there is no 
perturbative decay of the produced x particles. How- 
ever, immediately after adiabaticity is restored, t > t^, — 



71' 



-1/4 



Xi particles decay into lighter particles (i.e. 



those particles that have no gauge coupling to the in- 
flaton). In the case of LiL2e^ inflaton these are the 
(s)quarks, Higgs(ino), and L3, ei, 62 (s)leptons. 

Thus the fraction that is transferred from the inflaton 
to Xi's, and through their prompt decay into relativistic 
squarks, at every inflaton zero-crossing, can be computed 
analytically, and they are given by. 



2 

' 0.0067 g^e 0.0166 5y( 



(114) 

The total number of rf.o./ coupled to the LiL2e3 inflaton 
is 32 (4 from scalars, 4 x 3 = 12 from gauge fields, and 
4 X 4 = 16 from fermions). Therefore the fraction of the 
inflaton energy density that is transferred to relativistic 
quarks and squarks, see Eq. ()114|) . for ~ qw ^ 0.6 , 
has to be multiplied by (l-|-3-f4)=8 ( Allahverdi et all 
l2011bD : 



Proi//30 ^ 10.6% (per zero — crossing). 



(115) 



Note that this fraction is independent from the ampli- 
tude of oscillations. The draining the inflaton energy is 
quite efficient, nearly 10% of the infiaton energy density 
gets transferred to the relativistic species - but not all 
the SM do./ are in thermal equilibrium after one oscilla- 
tion. It takes near about 120 oscillations to reach the full 
chemical and kinetic equilibrium via processes requiring 
2 o 2 and 2 o 3 interactions. However due to the hi- 
erarchy between Hinf ~ 10~^m^, this happens within a 
single Hubble time after the end of in flation. One can 



estima te the flnal reheat temperature ([Allahverdi et ali . 
l2011bD 



1/4 1/4 
Po 



2 X 10* GeV . 



(116) 



where = 228.75 and po = (4/15)TO^(/)g, see Eq. 



III. MATTER-ANTI-MATTER ASYMMETRY 

If (p)reheating can provide a thermal bath where all 
the SM quarks and leptons are excited, it is then an im- 
portant question to ask - why the present day galax- 
ies and intergalactic medium is primarily made up of 
baryons rather than anti-baryons? 



The baryon abundance in the universe is denoted by 
fib = Pb/Pc, which deflnes the fractional baryon density 
Pb with respect to the critical energy density of the uni- 
verse: pc = 1.88 /i^ X 10~^^ g cm~'^. The observational 
uncertainties in the present value of the Hubble constant; 



Ho = 100 h km • ■ Mpc 



coded in h = 0.73 ( Kessler et ali l2009l ). It is useful to 
write in terms of the baryon and photon number densities 



/i/3000)Mpc- 



are en- 



rib 



2.68 X IQ-'^Qbh^ 



(117) 



where Ub is the baryon number density and is for 
anti-baryons. The photon number density is given by 
= (2C(3)/7r^)T'^. The best present estimation of 
the baryon density comes from BB N, which is ba s ed on 
SM physics with 3 neutr ino species (|Cvburt et all 120081: 
iFields and Sarkail I2006D 

0.019 < Qbh^ < 0.024 (95%CL) , (118) 
5.1 X 10"^° < ?7 < 6.5 X 10"^" (95%CL). (119) 

The observational data on D and "^He are consistent with 
each other and the expectations from the BBN analysis, 
but both prefer slightly higher value compared to the ^Li 
abundance Li/H\p = (1.7 ± 0.06 ± 0.44) x 10"!°, which 
is smaller than D and '^He by at least ^ 4.2cr. The ''Li 
abundance is primarily measured in the stellar systems 
such as globular clusters. 

From the acoustic peaks of the CMB the baryon 
fraction can be deduced. The WMAP data imply 
Vlbh'^ = 0.02273 ± . 00062 or r/ = 6.23 ± 0.17 x 
10-1" ()Komatsu et all l201ll) . The WMAP data relies 
on priors and the choice of number of parameters, it is 
possible to yiel d baryon abundance as low as Q,bh? = 
0.0175 ±0.0007 (|Hunt and Sarkail I2OI0I) . In spite of sys- 
tematic uncertainties the WMAP data is consistent with 
that of theoretical predictions from BBN. The major un- 
resolved problem is the Li abundance, stellar Li/ H mea- 
surements are inconsistent with both WMAP and BBN 
data, and this could be an useful probe of new physics at 
BBN, see for a review ( Jedamzik and Pospelovl i20Q9t ). 

Often in the literature the baryon asymmetry is given 
in relation to the entropy density s — 1.8(/*n^, where 
(7* measures the effective number of relativistic species 
which itself a function of temperature. At the present 
time sa 3.36, while during BBN 5* « 10.11, rising up 
to 106.75 at T » 100 GeV. In the presence of super- 
symmetry at T ^ 100 GeV, the number of effective rela- 
tivistic species are nearly doubled to 228.75. The baryon 
asymmetry defined as the difference of baryon and anti- 
baryon number densities relative to the entropy density, 
is bounded by 

58.4(7.2) X 10-" < < 9.2(9.2) x 10"" , (120) 
s 

at (95 %CL), where the numbers are 

CMB (jKomatsu et al. I, l201lh . and BBN 
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bounds (jCvburt et al\ . 120081: iFields and Sarkail . l2006l) . 
respectively. 

If at the beginning 77 = 0, then the origin of this 
small number can not be understood in a CPT invariant 
universe by a mere thermal decoupling of nucleons and 
anti-nucleons at T '-^ 20 MeV. The resulting asymmetry 
wou ld be too small by at lea st nine orders of magnitude, 
see (iKolb and Turneilll988D . Therefore it is important 



to seek mechan isms for generating baryon asyrnmetry , 
for reviews se e ("Dine and Kusenko( 120041: iRiottol . Il998l: 
iRubakov and Sliaposhnikov. ,1996^ 



A. Requirements for baryogenesis 



As pointed out first by Sakharov (jSakharovi Il967() . 
baryogenesis requires three ingredients: (1) baryon num- 
ber non-conservation, (2) C and CP violation, and (3) 
out-of-equilibrium condition. All these conditions are be- 
lieved to be met in the very early universe. 

c. Baryon number non-conservation: In the SM, 

baryon number B is vio l ated by n on-perturbative instan- 
't Hooftl . 



ton processes 



11976; 



Due to chiral anoma- 



lies both baryon number Jp and lepton number cur- 
rents are not conserved ( Adleil . Il969l : I Bell and Jackiwl 
I1969D . However the anomalous divergences come with 
an equal amplitude and an opposite sign. Therefore 
B ~ L remains conserved, while B + L may change 
via processes which interpolate between the multiple 
non-Abelian vacua of SU{2). The probability for the 
B + L violating tr ansitio n is however exponentially sup- 
pressed ( ^t Hooftl . Il976ai lbh. but at finite temperatures 
when T ^ Mw , baryo n violating transitions are in fact 



copious (lMantonl . lT983h . 

The B violation also leads to proton decay in 
GUTs (jLangackeii Il98ll) . The dimension 6 operator 
{QQQL)/A generates observable proton decay unless 
A > 10^5 GeV. In the MSSM the bound is A > 10^6 GeV 
because the decay can take place via a dimension 5 op- 
erator. In the MSSM superpotential there are also terms 
which can lead to AL = 1 and AB — 1. Similarly there 
are other processes such as neutron-anti-neutron oscilla- 
tions in SM and in SUSY theories which lead to AB = 2 
and AB = 1 transitions. These operators are constrained 
by the measurements of th e proton lifetime, which yield 



the bound Tp > 10^^ years ( Nakamura et a/.l . l2010f ). 



d. C and CP violation: The maximum C violation oc- 
curs in weak interactions while neutral Kaon is an ex- 
ample oi CP violation in the quar k sector which has a 
relative strength 10"^ (Nakamu ra et all [2010). CP 
violation is also expected to be found in the neutrino 
sector. Beyond the SM there are many sources for CP 
violation. An example is the axion proposed for solv- 
ing th e strong CP problem in QCD ( Peccei and Quinnl . 



e. Departure from thermal equilibrium: If B- 

violating processes are in thermal equilibrium, the in- 
verse pr ocesses will wash out the pre-existing asymmetry 
(Anfc)o (IWeinberd . IT979h . This is a consequence of S- 
matrix unitarity and CPT-theorem. However there are 
several ways of obtaining an out-of-equilibrium process 
in the early universe. Departure from a thermal equilib- 
rium cannot be achieved by mere particle physics consid- 
erations but is coupled to the dynamical evolution of the 
universe. 

1. Out-of-equilibrium decay or scattering: The 

condition for out-of-equilibrium decay or scattering is 
that the rate of interaction must be smaller than the ex- 
pansion rate of the universe, i.e. T < H. The universe in 
a thermal equilibrium can not produce any asymmetry, 
rather it tries to equilibrate any pre-existing asymmetry. 

2. Phase transitions: They are ubiquitous in the 
early universe. The transition could be of first, or 
of second (or of still higher) order. First order tran- 
sitions proceed by barrier penetration and subsequent 
bubble nucleation resulting in a temporary departure 
from equilibrium. The QCD and possibly electroweak 
phase transitions are examples of first order phase 
transitions. The nature and det ails of QCD phase 
transition is still an open Rebate (iKarsch e-t all (20011: 



transition is still an open ae bate l|Karscn et all l/UUlt 
iRaiagopal and Wilczekl . Il993l) . Second order phase tran- 
sitions have no barrier between the symmetric and the 
broken phase. They are continuous and equilibrium is 
maintained throughout the transition. 

3. Non-adiabatic motion of a scalar field: Any 
complex scalar field carries C and CP, but the symme- 
tries can be broken by terms in the Lagrangian. This can 
lead to a non-trivial trajectory of a complex scalar field 
in the phase space. If a coherent scalar field is trapped 
in a local minimum of the potential and if the shape of 
the potential changes to become a maximum, then the 
field may not have enough time to readjust with the po- 
tential and may experience completely non-adiabatic mo- 
tion. This is similar to a second order phase transition 
but it is the non-adiabatic classical motion which prevails 
over the quantum fluctuations, and therefore, departure 
from equilibrium can be achieved. If the field condensate 
carries a global charge such as the baryon number, the 
motion can charge up the condensate . This is the ba- 
sis fo r the Affleck-Dine baryogenesis ( Affleck and Dinj . 
I1985D . 



B. Sphalerons 

At finite temperatures B + L violation in the SM 
can be large due to sphaleron transitions between de- 
gener ate gauge vacua with diffe r ent Chern-Simo n s num - 
bers (jKlinkhamer and Mantonl 11984 iMantonl . Il983[) . 
Thermal scattering produces sphalerons which in effect 
decay in B + L non- conserving way s below 10^^ GeV 
( Bochkarev and Shaposhnikovl . Il987 ). and thus can ex- 
ponentially wash away B + L asymmetry. The three im- 
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portant ingredients which play important role are follow- 
ing. 

f. Chiral anomalies: In the SM there is classical con- 
servation of the baryon and lepton number currents Jj^ 
and J'l, but because of chiral anomaly ( at the quan- 
tum level) the current s are not conserved (lAdleil . Il969l: 
iBell and Jacki^ . Il969t ) . Instead (|'t Hooftl . Il976bl ). 



Sir 

9,J'l = -^N,Wrm,. + ^N,F'^f^F^p, (121) 



Wtt' 



where Ng is the number of generations, a2 and ai (Wi 



and F^^) are respectively the SU{2) and C/(l) gauge cou- 
plings (field strengths). Note that at the quantum level 
B + L is violated, but S — L = is still conserved. 

g. Gauge theory vacua: in the SU{2) gauge group, 
the vacua are classified by their homotopy class {0„(r)}, 
characterized by t he winding number n which label s the 
't Hooftl . Tl976at iPolvakovi 119771 ). 



so called 0-vacua 
gauge invariant quantity is the difference in the winding 
number (Chern-Simons number) 



N, 



a2 



cs 



fj.l' ■ 



(122) 



In the electroweak sector the field density WW is related 
to the divergence B + L current. Therefore a change 
in B + L reflects a change in the vacuum configuration 
determined by the difference in winding number 



A{B + L) 



d^xwrwa. 



'2NgN< 



cs ■ 
(123) 

For three generations of SM leptons and quarks the min- 
imal violation is A(_B -f L) = 6. Note that the proton 
decay p e+7r° requires A{B + L) — 2, so that despite 
_B-violation, proton decay is completely forbidden in the 
SM. The probability amplitude for tunneling from an n 
vacuum at i — )• — oo to an n + Ncs vac uum at t — » -|-oo 
can be estimated by the WKB method (|'t Hooftl . Il976al ) 



P{Ncs)b+l ~ exp 



-A7tNcs \ 
a2{Mz) J 



10 



-162Afcs 



(124) 



The baryon number violation rate is negligible at zero 
temperature, but as argued at fi nite temperatu r es the 
situation is c ompletely different ( Kuzmin et all Il985t 
lMantonl . [l983h . 

h. Thermal tunneling: below the critical temperature 
of the electroweak phase t ransition, the sphale ron rate is 
exponentially suppressed ( Carson et aLl . [l990l ): 



critical temperature the rate is however unsuppressed. 
Since the Chern-Simons number changes at most by 
ANcs ^ 1, one can estimate from Eq. (|122|) that 
ANcs - 9¥%hWf ^ I ^W, ^ {l/92lsph)- There- 
fore a typical energy of the sphaleron configuration is 
given by Esph l^hi^^^)^ ^ (l/sl^sp/i). At tem- 
peratures greater than the critical temperature there is 
no Boltzmann suppression, so that the thermal energy 
(X T > Esph- This determines the size of the sphaleron: 
Isph > l/sl^ Based on this coherence length scale one 
can estimate the baryon number violation per volume 
^ 'sph' ^^"^ P'^^' time ^ hph- On dimensional grounds 
the transition probability would then be given by 



- sph 



{l/l%^t) ^ Kia^TY 



(126) 



where k is a constant which incorporates various un- 
certainties. However, the process is inherently non- 
perturbative, and it has been argued that damping of the 
magnetic field in a plasm a suppresses the sph aleron rate 
by an extra power of a2 ( Arnold et all . Il997l ). with the 
consequence that Ygph ~ o^T^ . Lattice simulations with 



hard thermal loops also give Y sph ^ C'(10)q|T'^ ( Moord . 
I1999D . 

i. Washing out B + L: Assuming that in the early uni- 
verse the SM d.o.f are in equilibrium, the transitions 
ANcs — +1 and ANcs = ^1 are equally probable. 
The ratio of rates for the two transitions is given by 
Tgph +/Ysph - = exp(— A_F/T), where AF is the free en- 
ergy difference between the two vacua. Because of a finite 
B + L density, there is a net chemical po tential hb+l - 
Therefore one obtains ( Bochkarev and Shaposhnikovl 
I1987D 



duB^ 



dt 



= r 



sph 



sph 



,f ^sph 



(127) 



It then follows that an exponential depletion of ub+l due 
to sphaleron transitions remains active as long as 



r 



sph 



T3 



> H 



T <a: 



Mp 

2^72 

9* 



10^^ GeV . 



(128) 



This resuh imply that below T = 10^^ GeV, the 
sphaleron transitions can wash out any B + L asymmetry 
being produced earlier in a time scale r ~ {T^ /NgTsph). 
This seems to wreck GUT baryogenesis based on B — L 
conserving groups such as the minimal SU{5). 



C. Mechanisms for baryogenesis 



2.8 X IO^kT' 



V47r/ 



Esph(T) 

B{y/i) 



(125) 



where k is the functional determinant which can take the 
values lO"'' < n < 10"^ (|Dine et aUll992l) . Above the 



There are several scenarios for baryogenesis, the main 
contenders being GUT baryogenesis, electroweak baryo- 
genesis, leptogenesis, and baryogenesis through the decay 
of a field condensate, or Affleck-Dine baryogenesis. Here 
we give a brief description of these various alternatives. 
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1. GUT-baryogenesis 

This model relied on out-of-equilibrium decays of 
heavy GUT gau ge bosons X,Y — » qq, and X, F — » q l, 



for reviews see (|Dolgovl Il992l: iKolb and Turned . Il988l ). 

The decay rate of the gauge boson goes as Fx ~ axMx, 

1/2 

where Mx is the mass of the gauge boson and is 
the GUT gauge coupling. Assuming that the universe 
was in thermal equilibrium at the GUT scale, the decay 
temperature is given by 



Td 



-1/4 1/2 

9* <^x 



(MxMp)^^^ 



(129) 



which is smaller than the gauge boson mass. Thus, at 
T « Td, one expects nx ~ nx ~ n^, and hence the 
net baryon density is proportional to the photon num- 
ber density ns — ABn^. However below Tjj the gauge 
boson abundances decrease and eventually they go out- 
of-equilibrium. The net entropy generated due to their 
decay heats up the universe with a temperature which we 
denote here by Tj^. Let us naively assume that the en- 
ergy density of the universe at To is dominated by the X 
bosons with px ~ Mxnx, and their decay products lead 
to radiation with an energy density p = (7r^/30)5*T^, 
where ~ 0(100) for T > Mqut- Equating the expres- 
sions for the two energy densities one obtains 



nx ~ — Q* — — 
30^ Mx 



(130) 



Therefore the net baryon number comes out to be 
nB ABnx 3 Tr 



B = — 

s 



(j-jfi n^ 



—AB. (131) 

4 Mx 



Trh is determined from the relation w H^{Tj 
(7rV90)g*T4/M2. Thus, 



B : 



7* ^^'^axMp\ 



1/2 



M 



X 



AB. 



(132) 



Uncertainties in C and CP violation are now hidden in 
AB, but can be tuned to yield total B ^ 10^^° in many 
models. 

Above we have assumed that the universe is in thermal 
equilibrium when T > Mx- This might not be true, since 
for 2 o 2 processes the scattering rate is given by F ~ 
a^T, which becomes smaller than H at sufficiently high 
temperatures. Elastic 2 — > 2 processes maintain thermal 
contact typically only up to a maximum temperature ~ 
10^^ GeV, whi le chemical equilibrium is lost already at 
r~ 10^2 GeV (|Elmfors et a/J . ll994l : lEnqvist and EskolaL 
|1990D . 



2. Electroweak baryogenesis 

A popular baryogenesis candidate is based on the elec- 
troweak phase transition, during which one can in prin- 
ciple meet all the Sakharov conditions. There is the 



sphaleron-induced baryon number violation above the 
critical temperature, various sources of CP violation, and 
an out-of-equilibrium environment if the phase transi- 
tion is of the first order. In that case bubbles of bro- 
ken SU{2) X U{1)y are nucleated into a symmetric back- 
ground with a Higgs fie ld prof i le that cha nges through 
the bubble wall (K uzmin et all Il987l [l985l ). 

There are two possible mechanisms which work in dif- 
ferent regimes: local and non-local baryogenesis. In the 
local case both CP violation and baryon number vio- 
lation takes place near the bubble wall. This requires 
the velocity of the bubble wall t o be greater than the 
speed of the sound in the p l asma ( Ambiorn et aZ.I . [l990l : 
iTurok and Zadroznvl Il990l [ 199ll) . and the electroweak 
phase transition to be strongly first order with thin bub- 
ble walls. 

In the non-local case the bubble wall velocity speed is 
small compared to the sound speed in the plasma. In this 
mechanism the fermions, mainly the top quark and the 
tau-lepton, undergo CP violating interactions with the 
bubble wall, which results in a difference in the reflection 
and the transmission probabilities for the left and right 
chiral fermions. The net outcome is an overall chiral fiux 
into the unbroken phase from the broken phase. The fiux 
is then converted into baryons via sphaleron transitions 
inside the unbroken phase. The interactions are taking 
place in a thermal equilibrium except for the sphaleron 
transitions, the rate of which is slo wer than the rate at 
which the b ubble sweeps the space ( Cohen et all Il993t 
iNelson et a t. 1992). 

For a constant velocity profile of the bubble, v^, the 
net baryon asymmetry is generated by: 



nB 



- sph 



dt pb , 



(133) 



where /is is the chemical potential, which determines the 
tilt in the free energy of the sphaleron transitions, and 
numerically it is equivalent to: pB = piz — Vwt)/[(2N + 
5/3)r^]. Here p determines the profile of the bubble, and 
N denotes the number of Higgs doublets. The net baryon 
asymmetry can be calculated by following Eq. ([126^ : 



nB 
s 



100 



I? J- 



(134) 



where r is the transport time of the scattered fermions 
off the bubble wall, and = dzp{z). For the max- 
imum wall velocity ^ l/"\/3 and a typical: tT ^ 
10 — 1000 for top quarks, the maximum baryon asym- 
metry is given by: nB/s ~ lQ~^Fzl{vii,T^). The func- 
tion Fz also takes into account the CP phase, in the 
favorable scenario one would expect Fz/{vwT^) ~ 10~^. 
The d etails of the transport equations can be fou nd in 
Refs. (|Kainulainen et aLl . [20qTI: iNelson oZI . Il992l ). 



One great challenge for the electroweak baryogenesis is 
the smallness of CP violation in the SM at finite temper- 
atures. It has been pointe d out that an additional Higgs 
doublet ( McLerran et al\ . Il99ll : iTurok and Zadroznvl 
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Il99l[ ) would provide an extra source for CP violation 
in the Higgs sector. However, the situation is much im- 
proved in the MSSM where there are two Higgs dou- 
blets Hu and Hd, a n d two important sources of CP 
violation ( Ellis et oZI . Il982l ). The Higgses couple to 
the charginos and neutralinos at one loop level leading 
to a CP violating contribution. There is also a new 
source of CP violation in the mass matrix of the top 
squarks which caii give r ise to considerable CP violation 
(|Huet and Nelsonl . Il996l ). 

Bubble nucleation depends on the thermal tunneling 
rate, and the expansion rate of the universe. The tunnel- 
ing rate has to overcome the expansion rate in order to 
have a successful phase transition via bubble nucleation 
at a given critical temperature Tc > Tt > Tq. The effec- 
tive potential for the Higgs at finite temperatures can be 
computed, which takes the form: 



eff( 



7r03 + (A/4)04 (135) 



The order parameter is given by the ratio of {(j){Tc) /Tc) ^ 
7/ A, which has to be larger than one for first order phase 
transition. For Tc ~ 100 GeV, one obtains the condition 
for the sphaleron energy (R ubakov and Shaposhnikovl 
119961 : IShaposhnikovl [T987h 



Esph{Tc) 



> 7 log 



Esph{Tc) 



Tr 



-91og(10)-Mog(K). (136) 



which implies ( Bochkarev et aLl . Il99lh 



Tr. 



> 45 for K = 10"^ . 



In terms of the Higgs field value at T^, 

HTc) _ 9 EsphjTc) 1 EsphjTc) 
Tc ~ 47rB(A/52) Tc ^ 36 



(137) 



(138) 



where g is gauge couphng of S'C/(2)l, and B ^ 1.87. 
Then the bounds in Eqs. (|137p translate to 



^jTc)/Tc > 1.3, 



(139) 



which implies that the phase transition should be 
strongly first order in order that sphalerons do not wash 
away all the produced baryon asymmetry. This result is 
the main constraint on electroweak baryogenesis. 

Lattice studies suggest that in the SM the phase tran- 
sition is str ong ly first order only b e low Higgs mass Wff ^ 
72 G eV (jKaiantie et all Il996t iRummukainen etcdl . 
Il998[ ). Above this scale the transition is just a cross-over. 
Such a Hig gs mass is clearly exclu ded by the LEP mea- 
surements ( Nakamura et aLl . l201 0). thus excluding elec- 
troweak baryogenesis within the SM. However, this opens 
up a possibility to include new physics beyond the SM. 

a. Electroweak baryogenesis induced by new physics: 

it was pointed out that by modifying the SM 
Higgs self-interactions, especially the cubic term, it 



is poss ible to enhance th e first order phase tran- 
sition ( Anderson and Halll . I1992[K_ One such ex- 



ample has been cons i dered in ( Groiean et all 120051 : 
iMohapatra and Zhan 3. I1992D where non-renormalizable 
contribution to the Higgs potential has been considered 
of type: 



y($) = A (^$t$ _ !^ j +J_(^<i>t$_!i.j (140) 

where $ is the SM electroweak Higgs doublet and A is the 
scale of new physics which induces the corrections below 
the energy scale of A ~ Ojl) TeV. At zero temperature 
the CP-even scalar state can be expanded in terms of its 
zero-temperature VEV, (</>) = — 246 GeV, and the 
physical Higgs boson iJ: $ = 0/%/2 ^ jH + vo)/V2. 

The finite temperature effects are taken into account 
by adding a thermal mass to the potential: 



2\ 3 



Vjcl>,T)^cT^^y2 + Vj^,0) 



(141) 



where c is given in the high-temperature expansion of the 
one-loop thermal potential: 

c^^(4y^ + 3g' + g'^ + 4^ - 12^) , (142) 
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where g and g' are the 5[/(2)l and f7(l)y gauge cou- 
plings, and yt is the top Yukawa coupling. Note that 
there is no trilinear term in the effective potential. 

The critical temperature Tc at which the minima 7^ 
and (p = are degenerate is given by 



16cA2w4 



(143) 



The VEV of the Higgs field at the critical temperature 
in terms of mn, A and vq is 



(144) 



From Eqs. (|143p and (|144p . one finds that for any given 
mH there is an upper bound on A to make sure that 
the phase transition is always first order (w^ > 0), and 
there is a lower bound on A to make sure that the T = 
minimum at 7^ is a global minimum jTc > 0). 
These two combinations chart out a region where the 
phase transition is indeed first order: 



< A < V3^ (145) 



where rUc = voy^j4yf+3g^ + g^^)/8 w 200 GeV. In or- 
der to ensure that the thermal mass correction is positive: 
c > — >■ A > V^VQ/^ymff + 2m^. For these ranges of 
A the ratio Vc/Tc > 1, ensuring a successful sphaleron 
transition for the Higgs mass uih > 115 GeV. One nice 
aspect of this model is that the non-renormalizable scale 
A can also be constrained from the precision electroweak 
obser vable, which can be tested in near future by the 
LHC (jGroiean et a/.l . l2005l) . 
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3. Electroweak baryogenesis in MSSM 



4. Electroweak baryogenesis beyond IVISSiVI 



In the MSSM the ratio $(Tc)/Tc can be increased 
by virtue of the scalar loops which can make the cu- 
bic term in the temperature dependent Higgs potential 
large; Veffiip,T) = {-^i^ + aT^)^^ „ + 
In particular the right handed stop tjf coupling to the 
Higgs with a large Yukawa coupling. This leads to 
a strong first order phase transition - as the ratio of 
^{Tc)/Tc ^ 7/A > 1, where 7 determines the order 



paraineter |Carena et al . 1996 : Cline a l. 1998, 20Q3 



iLaind . Il996l : iLaine and Rummukainenl . IT998() . 

The finite temperature cubic term is given by: jT(p^ ~ 
(r/47r)[m~ (v3,T)]3/^, where the lightest right handed 
stop mass 



CT^ + O.lbMl cos(2^) + mM 1 - ^ 



(146) 

where At = At — fj,/ tan(/3) is the stop mixing parameter, 
At is the trilinear term in the MSSM superpotential, and 
fi is the soft-SUSY breaking mass parameter for the right- 
handed stop. The coefficient 7 of the cubic term jTip^ 
in the effective potential reads 



IMSSM ~ 7SM + 



Z/2 



1 - . (147) 



and can be at least one order of magnitude larger than 
7SM- The implications for the particle spectrum are: 

• A light right-handed stop: 120 GeV < < 
170 GeV < TOf 

• A heavy left-handed stop: tuq^ > 2 TeV. 

• A light SM-like Higgs: mn < 120 GeV, for 5 < 
tan/3 < 10. 

The present LEP constraint o n the lightest CP-eve n 
Higgs mass is mjj > 115 GeV (jNak amura et all . I2OIOD . 
Note that within MSSM, the lightest Higgs mass is 
bounded by: mfj < M| cos^ 2/3. Hence, even an MSSM- 
based electroweak baryogenesis may be at the verge of 
being ruled out. 

MSSM also provides new CP violating complex phases 
in the Higgsino sector, i.e. arg{iiAIi^2) > 10~^, with 
fJ; Ml, 2 < 400 GeV. The CP-violating phases are 
also constrained by the electric dipole moments. To 
match the observational limit on \de\ < 1.6 x lO^'^^ 
e cm (|Regan et aZ.I . |2002[ ). one requires first and second 
generation sfermion masses greater than 10 TeV. while 
the 2-loop electron dipole moment contribution comes 
out to be: |de| > 2 x 10"^^ c cm. 

The definitive test of the MSSM based electroweak 
baryogenesis will obviously come from the Hig g s and 
the stop searches at the LHC ( Carena et all l2003t 
IChung et a/.l . l2009D . 



Some of these problems of MSSM can be resolved 
in uMSSM (next-to minimal SUSY SM), with the help 
of introducing an extra singlet in the MSSM superpo- 
tential: W = m^S + XSHuHd + Wmssm- The S 
field gets a VEV to explain the /x = A (5) -term, but 
it also generates a singlet tadpole - its contribution to 
the vacuum energy, 6V = t^S - (l/167r2)"(S'/Mp)p2^ 
can be suppressed with the help of di screte symmetries , 
or Z^, where ^ rUsoftMp (lAbel et all Il995t 
iPanagiotakopoulos and Tamvakisl Il999l) . As a result the 
soft-SUSY breaking Higgs potential becomes: 

Vsoft = ts{S+h.c.)+m^,\S\^+ax{SHuHd+h.c.)+VMSSM , 

(148) 

Note that the trilinear, a\SHuHd now contributes to the 
7-term at the tree level, indicating potentially stronger 
first order phase transition even without a light stop and 
for niH > 120 GeV. The CP phases are distributed in 
gaugino masses as well as in the singlet, but not in the 
tree level of a\. 

One can similarly proceed with 4 SM singlets, and the 
Higgs doublet as in the ca se oi U(l)' e l ectrow eak baryo- 
genesis di scussed in Ref. ( | Kang et al. 1, I2005D . for a re- 
view see (|Kang et all l2009l) . where the superpotential 
contains: 



W = hSH^Hd + XS1S2S3 + Wmssm ■ 



(149) 



It is assumed that the C/(l)' is broken at higher VEVs, 
such as 1 — 2 TeV, and then the electroweak symmetry 
is broken at lower scales. The singlets Si, S2, S3 have 
VEVs greater than those of S and, Hu and Hd. The 
mass of Z' bosons are AIz' ~ C>{1) TeV. The tree level 
Higgs potential can now cont ain CP violating contrib u- 
tions from the phases A, 1^2 (|Kang et aLl . l2009l . l2005h : 

3 

Vsoft = Vmssm + ml\S\^ + ^\ 

i=\ 

-2A/i|5||ff°||-ffd|cos/33-2AAA|5i||52||53|cos/34 
-2m\s, \S\\Si\ cos /3i - 2m^<;^ \S\\S2\ cos ^2 
-2|m|^sJ|^i||^2| cos(-/3i + /32 + 7) (150) 

The potential can yield strong first order phase transition 
without large stop masses, and the new contributions 
to electron dip ole moments can be tam ed by tuning the 
Yukawa sector (|Kang et aLl . l2009l . l2005t) . 



5. Thermal Leptogenesis 

At temperatures 10^^ Qgy > T > 100 GeV, the P + L 
is completely erased by the sphaleron transitions, a net 
baryon asymmetry in the unive rse can stiU be generated 
from a non- vanishing B — L (iFukugita and Yanagidal . 
119861 : iHarvev and Kolbl . Il98ll : |Luty|, [199^ , even if there 
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were no baryon number violating interactions. The lep- 
ton number violating interactions can produce baryon 
asymmetry, a process whi ch is known as l e ptogen e- 
sis, for recent review s, see (|Buchmuller et all l2005allbl : 
iDavidson et al. 1. 120081) . 

The lepton number violation requires physics beyond 
the SM. The most attractive mechanism arise s in 5*0(10) 
which is left- right symmetric (for details, see (jLangacker , 
and has a natural foundation for the see-saw 
mechanism (iGell-M ann and S lanskv . 198rt 'Min kowski 
119771: iMohapatra an d ScnjanovlJll980tlYanagida|, \197^ 
as it incorporates a singlet right-handed Majorana neu- 
trino Nn with a mass Mn. A lepton number violation 
appears when the Majorana right handed neutrino de- 
cays into the SM lepton doublet and Higgs doublet, and 
their CP conjugate state through 



Nr^ H + l, Nr-^ H + l, 



(151) 



where Y}vi is the abundance of the right handed Majo- 
rana neutrino Ni and k is a thermal wash-out factor, 
which takes into account that the scatterings such as 
£H o £H tend to wash out any lepton asymmetry being 
created. 

In order to process the total lepton asymmetry 
into baryons, we need to know the che mical poten- 
tials (jKhlebnikov an d Shaposhnikovi |l988| ) 

i i 

(156) 

where i denotes three leptonic generations. The Yukawa 
interactions establish an equilibrium between the differ- 
ent generations (/i/i = /i/ and — ^q, etc.), and one 
obtains expressions for B and L in terms of the number 
of colors TV = 3, and the number of charged Higgs fields 
Nh 



where {H) I is the SM (Higgs) lepton. The relevant L 
violating interaction is then given by 

C D ^{Mn)uN,N, + y,jN~ijiT2H* + h.c. , (152) 

where i, j — 1,2,3. The above interac- 
tion is also responsible for generating the observed 
neutri no masses via t he c a nonical seesaw mecha- 
nism (Akhmcdov et al!, '2003': 'Buchmu ller et Ml, l2002t 
[Buchmuller and Pluma chcr, 199 g, 2000.), as required by 
the neutrino oscillation data (IGonzalez-Garcia et all . 



I2OIOD . This mass turns out to be rrii, « lyp-y^/M^v with 
V = 174 GeV, what implies right-handed neutrino mass 
scale of AfAf ~ O(10") GeV for \y\ ~ 1 and ~ 0.1 eV. 

Assuming a normal hierarchy in the heavy right 
handed neutrino sector. Mi <C A/2, M3 (correspond- 
ing to A^i, N2, N3). The CP asymmetry can be es- 
timated from the Ni decay, the asymmetry is gener- 
ated through the interference between tree level and 
one-loop d iagrams, which is given by (ICo vi et al, 1996; 
Flanz et aL . A995t iFukueita and Yanagida . ,1986; .Luty. . 



f lanz et ai.^ aijiJOC li'uKi 
19921 : iPlumacheil I1998D 



r(A^i ^ iH) - r(iVi ^ iH) 

T{Ni ^ IH) + T{Ni IH) 

-— y 



(153) 



lm{{yy^)uff{M^/Ml) , (154) 



1,2,3 

where / is a function which represents radiative cor- 
rections. In the case of SM, f{x) = \/x[{x — 2)/(a; — 
1) + (x - l)ln(l + 1/x)], and in the case of MSSM, 
fix) = Vi[2/(x - 1) + ln(l + 1/x)]. 

Let us take an example of the SM where the CP phase 
can be labeled by, |e| = 3Mi / {16ttv'^) ^/Am^^ sin 5, 
where Am^^m is the atmospheric ma ss scale of light neu- 



trinos ( Gonzalez-Garcia et aZ.I . [2010l) and 6 is the effective 
CP violating phase. The total lepton asymmetry is then 
given by 



Vl = \e\YNiK 



(155) 



L = 



UN^ + 9NNh 
6N + 3Nh 



(157) 



together with a relationship between B and B — L 
( Khlebnikov and Shaposhnikovi [l988) 



B^[^^l±^]iB-L). 
22N + 13Nh ' 



(158) 



The final asymmetry is then given by _B = (28/79) (B — 
L) in t he case of SM and B = (8/23)(B - L) for the 
MSSM (jKhlebnikov and Shaposhnikovi 119881) 

The baryon asymmetry based on the decays of right 
handed neutrinos in a thermal bath has been com- 
puted within MSSM (|BuchmuUer et ali l2002l l2005at 
iGiudice et al\ . I2004D . where besides the right handed 
neutrinos the right handed (s)neutrinos also participate 
in the interactions. The decay of a RH (s)neutrino 
with mass Mi results in a lepton asymmetry via one- 
loop self-energy and vertex corrections, see Eq. (|153p . 
If the asymmetry is mainly produced from the decay of 
the lightest right handed states, and assuming hlevai- 
chical right handed (s)neutrinos Mi <C M2,M^, we will 
have ( Davidson and Ibarral |2003 ) 



77 ~ 3 X lG-i"«: ^"'3 



nil 



Ml 



, (159) 



V 0.05 eV J VlO^ GeV 

for 0(1) CP- violating phases {rrii,^ < m^^ 
are the masses of light mostly light handed neutrinos). 
Here k is the efficiency factor accounting for the de- 
cay, inverse decay and scattering processes involving 
the rig ht handed sta tes (iBuchmuller et all . l2002l l2003l 
|2005ariGiu dice et ali. \20m. 

A decay parameter K can be defined as 



K 



Ti 



H{T = Ml) 



(160) 



where Fi is the decay width of the lightest right handed 
(s)neutrino. UK < 1, the decay of right handed 
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states will be out of equilibrium at all times. In this 
case the right handed states, which are mainly produced 
via scatterings of the left handed (s)leptons off the top 
(s)quarks and electroweak gauge/gaugino fields, never 
reach thermal equilibrium. The cross-section for produc- 
ing the right handed (s)neutrinos is oc T^^ (Mi), when 
T > Ml (< Ml), and hence most of them are produced 
when T ^ Mi. The efficiency factor reaches its maxi- 
mum value for k ~ 0.1 when m^-^ =^ 10^^ eV. For larger 
values of m^^ it drops again, because the inverse decays 
become important and suppress the generated asymme- 
try. Producing sufficient asymmetry then sets a lower 
bound , Ml > 10^ ^eV ([BuchmuUer et al. I, I2OO2I I2OO3I 



I20n5al : iGiudice" et all . I2OO4I ). Successful thermal leptoge- 
nesis therefore requires that > 10^ GeV. 

a. Resonant leptogenesis: If the mass splitting be- 
tween, say Ml, M2, is comparable to their decay 
widths, the CP asymmetry resonantly gets enhanced, 
see Eq. (|153l) . For example, let us consider A'^i and 
N2 . The dominant contribution to the CP asymme- 
try arises in the mixing of N 1 and N2, and it is given 
bv (iPilaftsis and Underwoodl . |2004 120051: iPlumacheil 
119981) 



Im(?;ty)2^ 



M|)MiM2 



87r{y^y)ii (Mf 



M|)2 



(M2r2-Miri)2 ■ 

(161) 

Now, assuming Mi - M2 and Mi - M2 ~ Fi - r2 There- 
fore, a large L asymmetry can be produced even if the 
initial abundance of A^i and N2 is small. 

b. Flavored leptogensis: So far we have assumed that 
all the leptonic flavors, i.e. r, fi,e, behave alike in a 
thermal bath. Especially in a non-SUSY case, where we 
can imagine a thermal bath of SM d. o.f with a tempera- 
ture < 10^^ GeV, the t- Yukawa interactions are in ther- 
mal equilibrium, while temperatures below 10^ GeV the 
muon Yukawa interactions are faster than the expansion 
rate of the universe and the leptogenesis rate. Since the 
wash out factor k is inversely proportional to the lepton 
violating interaction rates, so each of the fla vored symme- 
tries is subject to its own wash out effect ( Abada et all . 
l2006al lbl). The above Eq. (|T55]) gets modified to: 



(162) 



where Uf corresponds to the number of active flavors 
participating in thermal interactions. For temperatures 
ranging 10^ GeV < T < lO^^ GeV the flavored leptoge- 
nesis c an enhance the net bar yon asymmetry by a factor 
2 or 3 (|Davidson et aLl . l2008l) . 



c. Dirac leptogenesis: It is possible that the decay of a 
heavy particle accompanied by the CP distributed the 
lepton number equally between left handed and right 



handed particles with the net lepton number zero. A spe- 
cific example will be when the decay gives rise to a neg- 
ative lepton number in left-handed neutrinos, and a pos- 
itive lepton number of equal magnitude in right-handed 
neutrinos. If the observed neutrinos are Dirac in nature 
with a small Yukawa couplings h ^ 10^^^, then the left 
and right handed neutrinos will not come to thermal equi- 
librium before the electroweak scale, 7? ~ F ^ T^/Mp ~ 
h^T. Since the sphalerons interact with the left-handed 
neutrinos, violating B + L and conserving B — L, part of 
the lepton number in left handed neutrinos get converted 
into baryon number. Ar lower temperatures, the universe 
contains a total positive baryo n number, t otal positive 



lepton number, and B — L — (jPick et a, 



t otal p c 
xLl. l2OO0h . 



d. Leptogenesis via scattering: If there exists a 
shadow world similar to the SM sector, but hidden, and 
the only mediator is the heavy singlet neutrinos TV, then 
it is possible to realize leaking the lepton numbe r from 
the hidden to the visible sect or ( Bastero-Gil et aLl 120021 : 
iBento and Berezhianil l200l[ ). Let us consider a simple 
interaction between hidden (lepton doublet I' and Higgs 
(f)') and visible sector fields (SM lepton doublet I, and the 
SM Higgs 0) via 

h^al^Naq^ + KANa<t>' + ^M.feA^.iVb + H.C. (163) 

where the Yukawa interactions are given by hia and h^a ■ 
After integrating out the heavy neutrinos N with a mass 
Ma = gaM, where M being the overall mass scale and 
ga are order one real constants, we obtain an effective 
dimensional 5 operators 



hhd)(h H liluchd) 

2M * ^'^'^^ M ''^^ 



2M ^ 



H.C. , (164) 



with coupling constant matrices of the form A — hg~^h^ , 
A' — h'g~^h''^ and D = hg~^h''^. Let us suppose 
that the reheat temperatures in both hidden, T^, and 
visible sector, Tr, are below M. The only way the 
two sectors can interact via the lepton number violating 
scatterings mediated by the heavy neutrinos N which 
stay out of equilibrium, since Tn ^ M. The CP 
phase can be obtained in l(j> o £'(f>' and icf) o £'0'. 
the net asymmetry is given by Act = 3JS'/327r^Af^, 
where J = lmTT[{h'^ h')g-^{h'< h)g-\h^ h)* g-^] is the 
CP- violation parameter and S is the cm. of energy 
square. The final B — L asymmetry of the universe is 
given by (jBento and Berezhianil l200lD 



B -L = 



ns-L 



AHs 

10^2 GeV 
M 



J R 

4 



109 GeV 



(165) 



where s is the entropy density, and for Yukawa constants 
spreading in the range 0.1 — 1 can achieve the right lepton 
asymmetry. 
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e. Non-thermal 

vario us scenarios 

sis ( AUahverdi and Mazumdar 



1999 



1991 



leptogenesis: There exist 
of non-thermal leptogene- 



2000bl: iGiudice et all Il999at 



20031: lAsaka et al 



Lazarides and Shafi 



Muravama et al. . 1 19931 ) which can work for 
Tr < Mn- One classic example is when the right 
handed sneutrino, a scalar field, with mass Afjv, domi- 
nates the energy density of the universe and decays into 
the SM leptons and Higgs to reheat the universe and 
simultaneously creating the lepton asymmetry. The CP 
asymmetry can be created again from the interference 
between a tree level and one-loop quantum corrections, 
which yields the net asymmetry: 



riL Pn 

V ^^n— 



• (166) 



A similar expression can be us ed if any sneutrino con- 



densate de c ays aft er inflation (iBerezhiani et all 12001 



Mazumdaii l2004allbl: iMazumdar and Perez-Lorenzanal 



2004al lbl: Postma and Mazumdaii 2004 ). in which case 
Tr is replaced by the decay temperature of the sneu- 
trino condensate, i.e. Tp. The right handed neutri- 
nos and sneutrinos could also be excited non-thermally 
during preheating if they couple to the in flaton, which 
would generate non-thermal leptogenesis (iGiudice et al\ . 
ll999aD . 

f. Soft leptogenesis: In a perfect SUSY preserving 
limit the mass and the width of the right-handed neu- 
trino and sneutrino would be the same. Let us consider 
a single generation, where the mass is Mm, and their 



width is given by F = F-^Mat 47r 



mM% 



Y^v^/Mn, where v ~ 174 GeV is the Higgs VEV, and Yn 
is the Yukawa coupling. However, in a realistic scenario 
we would expect soft SUSY brea king terms which would 
be relevant for soft-leptogenesis dAUahverdi et al. 1. 120031 
iD'Ambrosio et ad I2Q03I: iGrossman et aLl . I2003ll2005[ ): 



BMn ~ ~ 



AYLNH + h.c. 



(167) 



This model has one physical CP violating phase given by: 
(j) = arg(^i?*). The soft SUSY breaking terms introduce 
mixing between the sneutrino N and the anti-sneutrino 
in a similar fashion to the — and — sys- 
tems. The mass and width difference of the two sneutrino 
mass eigenstates are given by 



Am = \B\ 



AT 



2|A|r 



M^ 



N 



(168) 



The CP violation in the mixing is responsible for gener- 
ating the lepton-number asymmetry in the final states of 
the N decay. This lepton asymmetry is converted into the 
baryon asymmetry through the spha leron process. The 
baryo n to entropy ratio is given by I D'Ambrosio et~al\ . 
I2Q03D : 



^ = -10-3 a 



AT\B\ 



[A\B\ 



r2 



Msin0, (169) 
Mn 



where the efficiency parameter a depends on the mech- 
anism that produces the right-handed sneutrinos. In a 
thermal production, the largest conceivable value could 
be of ord er a ^ 0.1 for the light neutrino mass rrii, ^ 
10-3 eV (|D'Ambrosio et all |2003| ). It may be slightly 
challenging to fix the parameters to obtain the right lep- 
ton asymmetry either making |i3|/r or r/|i?| small. The 
above requir ement gives a non- t rivial constraint on the 
param eters ( D'Ambrosio et al\ . 120031 : iGrossman et al\ . 
I2Q03D : 



A ~ 10^ GeV, Mn < 10* GeV, B < IGeV, 0-1. 

. 

Small value of \B\ cannot be obtained in gravity me- 
diated SUSY breaking scenarios, but it might be pos- 
sible to_aiTaiige_with^^ mediated SUSY break- 
ing ( Grossman et a/.l . 12004 



6. Affleck-Dine Baryogenesis 

As we discussed already, withi n MSSM there ex- 
ists cosmoloqical ly flat directions ( Dine et all . Il996bl : 
iGherghetta et all Il996l ). Field fluctuations along 
such flat directions are smoothed out by infla- 



tion ( Enqvist and Mazumdaii l2003l) . which effectively 
stretches out any gradients, and only the zero mode of 
the scalar condensate remains. Baryogenesis can then 
be a chieved by the perturbative decay of a conden- 
sate (jAllahverdi and Mazumdaii . l2Q07al l2008l) that car- 
ries baryonic cha rge, as was first p ointed out by Affleck 
and Dine (AD) (|Affleck and Dinel . |T985). As we will 
discuss, the flat direction condensate can get dynami- 
cally charged with a large B and/or L by virtue of CP- 
violating self-couplings. 

In the original version (lAffleck and Dind . I1985D 
baryons were produced by a direct decay of the con- 
densate. It was however pointed out that in t he case 
of gauge mediated SUSY breaki ng (iKusenkol . Il997bt 
iKusenko and Shaposhnikovl . 1 19981 ), and in the case of 
gravity mediated SUSY breaking (lEnqvist et 

Ml, I2000I: 

lEnqvist and McDonaldl . Il998l [19991, l2000|) , that the AD 
flat direction condensate in most cases is not stable 
but fragments and ev entually fo r ms no n-topological soli- 
tons called Q-balls (IColemanl I1985D . In gauge me- 
diated SUSY breaking scenarios these Q-balls can be 
made a long lived dark matter ca n didat e (|Kusenko et all 
119981: iKusenko and Shoemakeii l2009l). For reviews 



see (Dine and Kusenkol 2004 lEnqvist and Mazumdaii . 
120031) . 

Since, SUSY is broken by the finite energy density of 
the inflaton, the AD condensate receives corrections in 
the case of F-term inflation. Let us consider a generic 
superpotential for the AD field given by Eq. (1751) . then 
the effective potential for the AD field will be given by 
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FIG. 4: The dynamical motion of AD condensate 
(j) = (f>i + i<j)2 for a gravity mediated case with d = 4 
(soHd) and d — Q (dashed) with the initial condition 
9, = -tt/IO. 



(iDine et all [l995bl Il996bl) 



2| J,|2 



aXriH- 



dAL 



d-3 



h.C. 



d-3 



h.C. 



(171) 



The first and the third terms are the Hubble-induced and 
low-energy soft mass terms, respectively, while the second 
and the fourth terms are the Hubble-induced and low- 
energy A terms. The last term is the contribution from 
the non-renormalizable superpotential. The coefficients 
|C/|, a, Ad ~ 0{1), and the coupling A « l/{d- 1)!. 
Note that low-energy term has a mass dimension. The 
a, yl-terms in Eq. (|17ip violate the global U{1) symmetry 
carried by (j>. If \a\ is 0(1), the phase of (0) is related to 
the phase of a through n9 + 9a — t^', otherwise 9 will take 
some random value, which will generally be of 0(1). This 
is the initial CP- violation which is required for baryoge- 
nesis/leptogenesis. The AD baryogenesis is quite robust 
and can occur even ii i presence of positively large H ubble- 
induced corrections ( Kasuva and KawasakiL [20061 ). 

At large VEVs the first term dictates the dynamics of 
the AD field. If C/ < , the absolute minimum of the 
potential is (/> = and during infiation the condensate will 
evolve to its global minimum in one Hubble time. On the 
other hand if C/ > 0, the absolute value of the AD field 
settles during inflation to the minimum given by \(t>\ ~ 
{HjMp^^y/'^^'^. After the end of inflation the minimum 
of the condensate evolves from its initial large VEV to 
its global minimum </) = 0, note that the dynamics is 
non-trivial when the condensate starts oscillating when 
H{t) ~ ~ O(IOO) GeV. The dynamics of the AD 
condensate is non-trivial as shown in Fig. 

If inflation is driven by D-term, one does not get the 
Hubble induced mass correction to the flat direction so 
that C/, a — 0. Also the Hubble induced a-term is ab- 
sent. However the Hubble induced mass correction even- 
tually dominates once D-term induced inflation comes to 



an end. 

The baryon/lepton number density is related to the 
dynamics of the AD field by 



(172) 



where /3 is corresponding baryon and/or lepton charge of 
the AD field. The equation of motion for the AD field is 
given by 



(j) + 3H(t) + dV{(j))/d(j)* = 0. 
The above two equations give rise to 

'dV{<j>) 



(173) 



nB.L + 5HnB,L — 2/31m 



= 2/3t^I-(«0') • (174) 



dM: 



The net baryon and/or lepton number can be obtained 
by integrating the above equation 

a\t)nB,L{t)^ma\^^ f a^i'mt')]" sm{0) dt' , 



(175) 

Note that "a" introduces an extra CP phase which can be 
parameterized by sin(5). Note that the asymmetry is not 
governed by the Hubble induced A term, the amplitude of 
the oscillations will be damped and so the A-term, which 
is proportional to a large power of 4> will become gradu- 
ally negligib le. The net baryon and / or lepton asymmetry 
is given by ( Dine et all . Il996b[ ) 



nB,L{tosc) 



^ 3(d- 3) "^'^'^° sin 261 sin 5, 

2(d--2) ,frf-3N2/(rf-2) 

^ 3(^-3) "^"^ ("i^A/p ) (176) 



where sin 5 ^ sin 20 « £>(!). When the inflaton de- 
cay products have completely thermalized with a reheat 
temperature Tr, the baryon and/or lepton asymmetry is 
given by 



nB.L 



1 



Tr 



6{d - 3)"^ M^m^, 



{m^M^-'f^'-'^ (177) 



where we have used H{tosc) ~ TT^(p, and s is the entropy 
density of the universe at the time of reheating. For 
d = 4, the baryon-to-entropy ratio is 



nB,L 



1 X 10 



-10 



and for d = 6 
nB.L 



10- 



-10 



X 13 



x/3 



1 TeV\ / Tj 



) ^109 GeVy ' 



1 TeV 



1/2 



R 



100 GeV, 



(178) 



(179) 
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where we have taken the net CP phase to be ^ 0{1). 
The asymmetry remains frozen unless there is additional 
entropy production afterwards. 

The lepton asymmetry calculated above in 
Eqs. (I178I179P can be transformed into baryon number 
asymmetry via sphalerons ns/s = (8/23)ni/s. AD 
leptogenesis has important implications in neutrino 
physics also, because in the MSSM, the LH^ direction 
is lifted by the d — A non-renormali zable operator which 
also Rives r ise to neutrino masses ( Asaka et al\ . l2000al : 
iDine aUllQQebl ): 



2M, ^ ' 



2(H„ 



(LiHu) , 



(180) 



where we have assumed the see-saw relation i — 
{Hu}"^ /Mi with diagonal entries for the neutrinos Ui, i ~ 
1,2,3. The final hb/s can be related to the lightest neu- 
trino mass since the flat direction moves furthest along 
the eigenvector of LiLj which corresponds to the small- 
est eig envalue of the neutrino mass matrix ( Asaka et all 
l2000aHDine et al\ . \l996\Ji ) . 



riL ^^_io a /^"*3/2 

— w 1 X 10 X p \ — — 

s \ 



108 GeV 



10-^ cV 



(181) 

where rrii, i denotes the lightest neutrino. Similarly 
the d = 6 case corresponds to the flat direc tion 
udd (|Enavist and McDonaldl . l2000l: iMcDonaldl Il997l) . 



7. Baryogenesis below the electroweak scale 

At temperatures below the electroweak scale, the 
sphaleron transitions are rather inactive. If the uni- 
verse reheats below the electroweak scale then it is a 
challenge to generate the required baryon asymmetry. 
The leptogenesis based scenarios are hard to imple- 
ment below the electroweak scale. However within SM 
there exists a possibility of rea lizing cold elect r oweak 
baryogenesis as discu ss ed in ( Cornwall et al\ . 2001 



Cornwall and Kusenkol l2000l: lEngvist et aLl. l2010b : 



Garcia-Bellido aU 1 19991: iKrauss and Troddenl . Il999 : 



Tranberg all l2010l : iTranberg and Smitl . I2003l l2006[ ) 

SUSY further opens a door to realize baryon asymmetry 
at temperatures even close to the B BN via R-parity vio- 
lating interactions (Clinc and Rabv. ll991 : Kitano et all . 



l2008riKohri all i2009i : iScherrer et ami99lh . Here we 



will discuss both the scenarios. 

a. Cold electroweak baryogenesis: There are mecha- 
nisms to obtain cold electroweak baryogenesis where it 
is assumed that the SM d. o.f are not in thermal equilib- 
rium. Moreover in a cold environment the CP-violation 
in SM is much larger than at the electroweak tempera- 
tures, and most of the baryon asymmetry is produced at 
the initial quench when the Higgs field is rolling down 
the potential. Baryon production essentially stops after 
the first few oscillations, after which the coherent Higgs 



field will start decaying, thereby reheating the universe. 
However, the hurdle is to obtain this fast quench without 
the presence of strong first order phase transition. 

There are couple of possibilities of realizing cold ini- 
tial condition and out of equilibrium condition. There 
could be a very low scale of inflation which might not 
be responsible for generating the seed perturbations, 
or the universe could be simply trapped in a vacuum 
where the SM d.o.f are not even excited. The out-of- 
equilibrium condition can be obtained during the co- 
herent oscillations of the scalar fields. In order to re- 
alize this idea, we would require a scalar field coupled 
to the SM Higgs, a^H^. During the coherent oscilla- 
tions, it is possible to ha ve the baryon number violat- 
ing sphaleron transi ti ons (ICornwall and Kusenkol '2000'; 



iGarcia-Benido et all \l99^ iTranberg and Smir r2003.) 



This can happen since the Higgs oscillations can ex- 
cite the electroweak gauge bosons from the time de- 
pendent vacuum fluctuations with a very large occu- 
pation number, similar to case of preheating. These 
long wavelength fluctuations of the gauge fields are re- 
sponsible for overcoming the sphaleron barriers which 
leads to the baryon number violation. Furthermore, 
there could be extra sources of CP-violations during 
the oscillations as pointed out in ([Cornwall et al. j_ 20011 : 
ICornwall and Kusenkol l2000l : ITranberg et all \20m . 



b. R-parity violation and baryogenesis: The current 
limits on some of the i?-parity violating interactions are 
poorly understood. Let us now consider a scenario where 
B and L are violated within MSSM, with a superpotcn- 
tial: 

^ = l^iLiHu + XijkLiLjCk + Xijf.LiQjdk + Xijj.Uidjdk , 

(182) 

where Li = (i^j,eO, = {u^.di), Ha = iK^f^^V: 
Hd = {h^,h'^)'^, etc are SU{2)l doublets and m^, 
are SU{2)l singlet quarks. In Ea. (|182p . the first three 
terms violate lepton number by one unit (AL = 1), 
while the last term violates baryon number by one unit 
{AB = 1). For the stability of proton we assume that 
Xijk — X'^jj. = 0. This can be accomplished if there exists 
any conservation of lepton number, which then forces /i'j 
to be zero. However, the electric dipole moment of neu- 
tron gives ( Barbier et all [20051 ) 



Im(A^'i2A:i'32) <0.03 



o.or 



M 
TeV 



(183) 



and the non-observatio n of n—n oscillation g ives an upper 
bound on X^^^ to be ( Barbier et aLl [20051) 



|A'i'u.| < (10- 




(184) 



While A332 is hardly constrained and can be taken to be 
as large as 0(1). Let us consider that a scalar field 4> 
decays to MSSM d.o./ right before BBN primarily into 
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gauge bosons and gauginos via i?-parity violating cou- 
plings A^^fc. 

Let us assume that the gauginos are heavier than the 
quarks and squarks. As a result their decay to a pair 
of quark and squark through one loop quantum correc- 
tion gives rise to a net CP violation. The magnitude of 
CP violation in the decay; g — > ti'^ can be estimated 
as (jCline and Rabv, 1991): 



r (.g ^ ti-) - r (g tt) ^ X-,, Im {Al^,mg) 



T^tot 

a 



I6n 



(185) 



where ^323 is the trilinear SUSY breaking term and 
we also assume a maximal CP violation. As a result 
the decay of gauginos produce more squarks (antisqarks) 
than antisquarks (squarks). The baryon number violat- 
ing {AB = 1) decay, induced by A323 of squarks (anti- 
squarks) to quarks (antiquarks) then gives rise to a net 
baryon asymmetry. Note that the decay of squarks (anti- 
squarks) are much faster than any other processes that 
would erase the produced baryon asymmetry. Hence the 
B-asymmetry can simply be given by: 



VB 



4 m0 



(186) 



where Eg ~ 0.5 is the branching ratio of the decay of (j) 
to gg, and in the above equation s is the entropy density 
resulted through the decay of (f). For Tn/m^ ^ 10"'' and 
^ 10^ GeV. Therefore a reasonable CP violation of 
order e ~ 0.01 — 0.001 could accommodate the desired 
baryon asymmetry o f 0(1O~^°) close to the temperature 
of r ^ 10-1 MeV (|Cline and RabvL Il99ll: iKohri" et all 
l2009( ). 



IV. DARK MATTER 

There is a conclusive evidence that a considerable frac- 
tion of the current energy density is in the form of a non- 
baryonic dark matter. The dynamical motions of astro- 
nomical objects suc h as r otation curves for spiral galax- 
ies (jBegeman et a^]. ll99l[ ). velocity dispersion of individ- 
ual galaxie s in galaxy clusters , large x-ray temperatures 
of clusters (|Flores et al. I. |200<D . bu lk flows and the pecu - 
liar motion of our own local group ( Dressier et aU[l987l ). 
all implies the presence of a dark matter. The mass 
of galaxy clusters inferred by their gravitational lensing 
of background images is als o consistent with th e large 
dark-to-visible mass ratios ()Bolton et al I |2006[ ). Per- 
haps the most compelling evidence, at a statistical sig- 
nificance of 8(7 comes from the two colliding clus t ers of 
galaxies, known as the Bullet cluster ( Clowe et al ]. 120061 : 
iMark evitch et al.. . ,2004) . It was found that the spatial 
offset of the center of the total mass from the center of the 
baryonic mass peaks cannot be explained with an alter- 
ation of gravitational force law. Furthermore, the large 
scale structure formation from the initial seed perturba- 
tions from inflation requires a significant non-baryonic 



dark matter component (jAbazaiian et al. I, |2009^ 
terms of the critical density, p, 
1.88 X 10"29g 
lOO/i kmsec-^Mpc 
from WMAP a nd large scale s t ructu re data is ^dm 
Pdm/Pc ^ 0.22 (jKomatsu et adl201lh . 



In 



and with Hubble constant Hq = 
the dark matter density inferred 



a weakly in- 
undiscovered. 



The dark matter is assumed to be 
teracting massive particle (WIMP), yet 
There are many w ell mot i vated particle physics candi- 
dates, e .g. (Berto ne et all. l2005t lJungman et all . \l99(k 

l2008li 



iKusenkol 120091 : iTaoso et al 



all of which arise 
from beyond the SM physics. The dark matter is as- 
sumed to be stable on the scale of cosmological structure 
formation. By virtue of new symmetries, for example 
R-parity conservation in SUSY allows the lightes t SUSY 
particle (LSP) t o be absolutely stable ( Ellis et aZl . ll984al : 
iGoldbe rg, '1983), or in the case of extra dimensions, the 
Kaluza-Klein (KK) parity leaves the li ghtest KK particle 
(LKP) stable ([Servant and Taitl . l2003l ). 



In many cases some of these symmetries which pro- 
tect the dark matter particle from decaying are broken 
by sufficiently suppressed higher-dimensional operators, 
such that the dark matter might as well have a finite 
life time comparable to the age of the universe. In the 
context of SUSY grand unification, operators with mass 
dimension 6 are expected to make SUSY dark matter 
unstable, with a time-scale 



87r 



mx 



TeVV 



mx 



M, 



GUT 



2 X 1016 GeV 
(187) 

where Mgut ~ 10^^ GeV, and rax is the dark matter 
particle. The lower dimensional operators would yield 
much shorter time scale, as it would lead to dark mat- 
ter decay long before the structure formation. Within 
SUSY one compelling candidate could be the gravitino 
with R-parity weakly broken in the hadronic sector, 
yield ing the required ba ryon asymmetry also in the pro- 
cess (|Kohri et aLl . l2009l ). 



The widely accepted lore is that after radiation-matter 
equality, when the universe becomes matter dominated, 
the density perturbations in the dark matter begin to 
grow, and drive the oscillations of the photon-baryonic 
fluid around the dark matter gravitational potential 
wells. Immediately after the epoch of recombination 
the baryons kinematically decouple from photons, which 
then free-stream through the universe; the baryons on 
the other hand slowly fall into the potential wells created 
by the dark matter particles, eventuall y becomiiig light 
emitting galaxy, see for more details ( Dodelsoiil . 120031: 
iKolb and Turneil . Il988t [Peeblesl I1994D . There are three 
broad categories of dark matter which have been central 
to our discussion. 
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A. Types of dark matter 

1. Hot Dark Matter 

If the dark matter particle is colhsionless, then they 
can damp the fluctuations from higher to lower density 
regions above the free-streaming scale. This hot dark 
matter consists of particles which are relativistic at the 
time of structur e formation and therefo re lead to large 
damping scales ( Bond and Szalavl Il983[ ). 

The SM neutrinos are the simplest examples of hot 
dark matter. In the early universe they can be decoupled 
from a relativistic bath at T ^ 1 MeV, leading to a relic 
abundance today that depends on the sum of the flavor 
masses: 



90 cV 



(188) 



Various observational constraints combining Ly-a for- 
est, CMB, Sup e r Nova e and Galaxy Cluste rs data leads 
to (|FogU et all 120081: ISeliak et all l2006l ): J2m>^ < 
0.17 eV (95 % CL). Similar limits can be applied to 
any generic hot dark matte r candidate, such as ax- 
ions (iHannes tad et al. j_ 20ip[l or to hot st erile neutri- 
nos (iDodelso ~t all 120061 : iKusenkol . l2009l). The free- 
strea ming length for neutrinos is ( Kolb and Turneil . 
I1988D : 



A 



FS 



20 



30 eV 



Mpc. 



(189) 



For instance, the universe dominated by the eV neutri- 
nos would lead to suppressed structures at 600 Mpc scale, 
roughly the size of supercluster. Furthermore, hot dark 
matter would predict a top-down hierarchy in the forma- 
tion of structures, with small structures forming by frag- 
mentation of larger ones, while observations show that 
larger galaxies have formed from the mergers of the ini- 
tially small galaxies. 



2. Cold Dark Matter 

The standard theory of structure formation requires 
cold dark matter (CDM), whose free-streaming length is 
such that only fluctuat ions roughly below t h e Earth mass 
scale are suppressed (iBertschingeii 2006 : Green et i 



2001: Loeb and Zaldarriaga, 



2004 1 120051: [Hofmann et al 
20051 ). The CDM candidates are heavy and non- 
relativistic at the time of their freeze-out from thermal 
plasma. The current paradigm of ACDM is falsifiable 
whose predictive power can be used to probe the struc- 
tures at various cosmological scales, such as the abun- 
dance of clusters at z < 1 and the galaxy-galaxy corre- 
lation functions have proven it a successful and widely 
accepted cosmological model of large scale structure for- 
mation. 

The N-body simulations based on ACDM provide a 
strong hint of a universal dark matter profile, with the 



same shape for all masses, and initial power spectrum. 
The halo density can be parametrized by: 



p{r) = 



Po 



{r/R,)-> [l + (r/i?,)«](^-'^)/" 



(190) 



where po and the radius Rs vary from halo to halo, the 
parameters a, /3 and 7 vary slightly from one profile to 
other. The four most popular ones are: 

• Nav arro, Frenk and White (NFW) pro- 
file ([Navarro et all Il997[ ). where a = 1, /3 = 
3,7 = 1, and = 20 Kpc. 



• Moore profile (jMoore et all Il999[ ). where a 
1.5, ^ = 3, 7 = 1.5, and Rs = 28 Kpc. 



• Kra profile ( Kravtsov et all Il998( ). where a — 
2, 13 = 3, 7 = 0.4, and Rs = 10 Kpc. 



• Modi fied Isothermal profile ( Bergstrom et all 
Il998( ). where a = 2, /3 = 3, 7 = 0, and i?^ = 
3.5 Kpc. 

Amongst all the four profiles, the scales where devia- 
tions are most pronounced (the inner few kiloparsecs) are 
also the most compromised by numerical uncertainties. 
The power-law index value, 7, in the inner most regions is 
part of the numerical uncertainties and still under debate, 
as all four simulations provide different numbers. The 
simulations hint towards a cuspy profile, as the density in 
the inner regions becomes large, while from the rotation 
curves of low surface brightness (LSB) galaxies point to- 
wards uniform dark matter density profile with constant 
density cores (Gentile et all l2004[ ). In our own galaxy 
the situation is even more murky, as the observations 
of the velocity dispersion of stars near the core suggests 
a supermassive black hole at the ce nter of our Galaxy , 
with a mass Msmbh « 2.6 x lO^M© (jGhez et all\im^ . 
Many galaxies have been found to host supermassive 
blackholes of 10^ — lO^M©. It has been argued that 
if supermassive blackhole exists at the galactic center, 
the accretion of dark matter by the blac k hole would 
enhance the dark matter density ( Peeblesl Il994l) . To 
alleviate some of these problems, da rk matter with a 
strong elastic scattering cross section ( Dave "oil . 1200 it 
ISpergel and Steinhardtl. 2000l). or large annihilation cross 
sections ( Kaplinghat et a/T2000( ) have been proposed. 

There are further discrepancies between observations 
and numerical simulations. The number of satellite ha- 
los as predicted by simulations exceeds the number of 
obser ved Dwarf ga l axies i n a typical galaxy like Milky- 
Way (|Klvpin et all Il999l: iMoore et all Il999l ). However 



recent hydrodynamical simulations with ACDM, includ- 
ing the supernovae induced outflows suggest a fall in the 
dark-matter density to less than half of what it would 
otherwise be within the central Kpc. 
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3. Warm Dark Matter 



Besides hot and cold dark matter, the early uni- 
verse can also provide warm dark matter (WDM) can- 
didates whose velocity dispersion lies between that of 
hot and CDM. The presence of WDM reduces the 
power at small scales due to lar ger free-streaming 
length compared to that of a CD M ( Bode et"aH . 1200 1[ 
ISommer-Larsen and Dolgovi l200l[) . 

The origin of WDM can be found within ster- 
ile states. For instance, the see-saw mechanism 
for th e active neutrino masses from the SM s inglet 
states (iGell-Mann and Slanskvl. Il980l: [Mink owski. Il977t 



iMohapatra and Senianovid . ll981 



lYanagidal Il979f ) would 
naturally generate masses to the active m{ui^2,3) ^ 
y'^{H)^ /Mn, and sterile neutrinos m{va) ^ Mn {a > 3) 
in Eq. (I152p . if we take ij = I, - --n + 3. The typi- 
cal mixing angles in this case are: 9ai ~ y'^^{H)'^ / . 
In order to explain the neutrino masses from atmo- 

2 is sufficient, how- 



spheric and solar ne utrino data, n = 
ever f or pul sar kicks ( Kusenkj . 120061: 
1996, 19 9^, supernovae explo sion 



Kusenko and Segre , 
Fryer and Kusenkoi 



(II 

2006i: ,Hidaka and Fulleil . l2006l l2007l), as well as sterile 



neutr i no as a dark matt e r candidate (lAbazaiian et al. 



20011: lAsaka et all 12005': Dod elson and W idro 
2002; Petraki and Kusenk 



Dolgov and Hansenl. . 

Shi and Fulleil . I1999D . we require at least n 



coj, 



1994; 



20081 



3, so 



total 6 sterile Majorana states, for a review on all these 
effects, see (.KusenkO; . .2009^ . The presence of such extra 
sterile neutrinos is a lso supported by Vg — » Ve oscillations 
observed at LSND ('Aguilar et a/.v200ll) , and the r ecent 
results by MiniBoone ( Aguilar- Arcvalo et al\ . [2OI0I) . 

A sterile neutrino with a KeV mass can be an ideal 
WDM candidate which can be produced in the early 
universe by oscillation/conversion of thermal active neu- 
trinos, with a momentum distrib ution significantly sup- 
prcsscd from a thermal spectr um (jAbazaiian et al, 2001; 
[Podclson and Widrow ^ 1994). A typical free-streamini 



scale is given by, see ( Abazaiian and KoushiappaE 1200 



B. WIMP production 

1. Thermal relics 

At early times it is assumed that the dark matter parti- 
cle, denoted by X is in chemical and kinetic equilibrium, 
i.e. in local thermodynamic equilibrium. The dark mat- 
ter will be in equilibrium as long as reactions can keep 
X in chemical equilibrium and the reaction rate can pro- 
ceed rapidly enough as compared to the expansion rate 
of the universe, H{t). When the reaction rate becomes 
smaller than the expansion rate, then the particle X can 
no longer be in its equilibrium, and thereafter its abun- 
dance with respect to the entropy density becomes con- 
stant. When this occurs the dark matter particle is said 
to be "frozen out." 

The equilibrium abundance of X relative to the en- 
tropy density depends upon the ratio of the mass of 
the particle to the temperature. Let us define the vari- 
able Y = nx/s, where nx is the number density of 
X with mass mx, and s = 2t{^ g^T^/A'b is the en- 
tropy density, where counts the number of relativistic 
d.o.f. The equilibrium value of Y , Yeq oc cxp(— x) for 
X = mx/T ^ 1, while Yeq ~ constant for a; ^ 1. 

The precise value of Yeq can be computed exactly 
by so lving the Boltzmann equation ( Kolb and Turneii . 
I1988D : 



fix + ZHnx = -{(Tv){nx - {n"^) ) 



(192) 



where dot denotes time derivative, a is the total annihila- 
tion cross section, v is the velocity, bracket denotes ther- 
mally averaged quantities, and n*^' is the number density 
of X in thermal equilibrium: 



n-^" =g(mT/27r)^/^e-"^/^, 



(193) 



where T is the temperature. In terms of Y = nx / s and 
X = mx/T^ and using the conservation of entropy per 
comoving volume [sa^ = constant), we rewrite Eq. (jl92p 
as: 



Afs ~ 840 Kpc h"^ 



1 KcY\ (<p/T> 
3.15 



(191) 



where is the mass of the sterile flavor eigenstate, 0.9 > 
(p/T) > 1 is the mean momentum over tempera- 
ture of the neutrino distribution and ranges from 1 (for a 
thermal) to ~ 0.9 (for a non-thermal) distribution. Very 
stringent bounds on the mass of WDM particles have 
been obtained by different groups. Typically, the bounds 
range from m, ^ ^0 " 20 KeV (95 % CL) {mwDM > 
2 — 4 KeV), see ( Kusenkoi |2009() . It is quite plausible to 
imagine a mixed dark matter scenario, where more than 
one species contributed to the total dark matter abun- 
dance. If there is a fraction of sterile neutrinos or WDM, 
then the above bounds can even be relaxed. 



dY _ {<7V)S „ (ye9)2) 



dx 



Hx 



(194) 



In the case of heavy X, the cross section can be expanded 
with respect to the velocity in powers of v'^ , (av) = a + 
b{v'^) + 0{{v'^)) + ... a + 66/x, where x ^mx/T and 
a, b are expressed in GeV~^. Typically a 7^ for s-wave 
annihilation, and a = for p-wave annihilation. We 
can rewrite Eq. (jl94p in terms of a new variable: A = 

Y - 

A' = - /(a;)A(2r"« + A) , (195) 

where prime denotes d/dx, and 



f{x) = :^TOxMp(a + 6b/x)x-^ . 



(196) 
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One can find a simple analytic solution for Eq. (|195l) for 
two extreme regimes 

A-^A' = , ^ » ^ , A' » Y'"' . (198) 



Integrating the last equation for (a;/,oo), and using 
A{xf) >• Aoo, we find 



45 



(199) 



2. Coannihilating WIMPs 

If there are N particles, Xi [i — 1, . . . ,N) with the 
lightest one, Xi, which have nearly degenerated masses 
TOi, such that mi < m2 < ■•• < m-jv-i < 
and internal d.o./ (statistical weights) gi. The next to 
lightest dark matter particle will be -/V2. In this case 
the abov e calculation of relic density, Eg. (fT^ ) gets 
modified fBinetruv et aLl . Il984 [G"riest and Seckei ll99lt 
[Servant and Taitl . 1200311 



dn 

— = -3iJn 
dt 



N 



ill J 



cq cq\ 



(203) 



In terms of the present energy density, px = mxnx = 
TTT'XSoYoo, where sq = 2889.2 cm~^ is the present entropy 
density, the relic abundance of dark matter particle in 
terns of the critical energy density is given by: 



where n = X^i^i number density of the relic 

particle, since all other particles decay much before the 
long-lived Xi. The total annihilation rate for Xi — Xj 
into a SM particle is given by: 

= Y^^iXiX, ^ Xsm), and (204) 

X 

= ^{Pi-p,Y-m1myEiE^, (205) 

is the relative particle velocity, with pi and Ei are the 
four- momentum and energy of particle i. One requires 
to define a thermal averaged (aijVij), which is defined 
by: 



nxh^ 



1.07 X 10 



9 



Mp 



GeV"^ , (200) 



where the freeze-out temperature is defined by solving 
this equation A(x/) — cY'^'^{xf) iteratively for early and 
late time solutions, for c ^ 0(1) 



In 



c(c -f2)V45 



gast mxMp{a + 6b/xf) 
27r5/2 



1/2 , — 
9* 



(201) 

An approximate order of magnitude estimation of the 
abundance can be written as: 



rixh^ 



3 X 10-2'^ cm^ s^i 0.1 pbarn 



(202) 



For a WIMP interacting with a heavy gauge boson, would 
naturally yield an upper bound on mx- From the above 
Eq. (|202p . on dimensional grounds, ft ^ l/{av) ~ I/ti^ 
for {av) ~ Q^{mx /M'^Y , where M is the mass of the 
new gauge boson. For mx ~ Af ~ 1 TeV, abun- 
dance of the dark matter becomes of order of unity, 
Vtxh? ^ for a ^ 0.1. Actually, the unitarity 

bound li mits the dark matter mass to be below mx < 
300 TeV (|Griest and Kamionkowskii Il990[ ). For a real- 
istic scenario a ^ 0.01, the unitarity bound would yield 
mx < 3 TeV. 

Note that the above non-relativistic expansion of 
(av) ss a + 6b/x may not hold universally. When a 
mass of second particle becomes nearly degenerate with 
the d ark matter particle X as in the case of coanni hila- 
tion (|Binetruv et all Il984 iGriest and Secket ll99lD . or 
the cross section is strongly varying function of the cen- 
ter of m ass energy as in the ca se of a resonant anni- 
hilation ( Griest and Seckel Il99ll) . In the latter case, tr, 
gets a boost by resonant annihilation when mx ~ mAl^^ 
where X annihilates with an exchange of particle, A, with 
a mass mA- 



f d^pid^pjfjj 



(206) 



where fi are distribution functions in the Maxwell- 
Boltzmann approximation. 

Typically, when the scattering rate of particles off SM 
particles in a thermal background is much faster than 
their annihilation rate, then in the above Eq. (jl92l) . (av) 
is replaced by: 



N 

(cToffw) = y^o- 



,^^{1 + Aif/^{1 + A,)3/2e-(^-+^.) 

9eff 



(207) 

where A^ = (m^ - mi) /mi, and geff = Ylf 9i{^ + 
Ai)^/^ exp{—xAi). In the case of co-annihilation, the 
freeze-out temperature is determined by 



= In 



c(c + 2)\/45 



geff mx Mpjaeff + Qbeff/xf) 
27r5/2 



1/2 , — 
9* Jx] 



(208) 

where a-eff and &e// are the coefficients of the Taylor 
expansion of (Jeff- The relic abundance for A'^i is now 
given by 



3 X 10-2^ cm3 s-i 

gV^x'j^iJa + ih/xf) 



where (209) 



la = Xf I aeffX dx , h ^ 2x 

Ixf 



f I beffX ^dx. 

Xf 
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In any realistic framework there are many particles which 
interact with the dark matter candidate X. They all 
eventually decay into X, and at the time of freeze- 
out the density of all heavy particles is exponentially 
suppressed except when there is mass degeneracy oc- 
curs between heavy particles and the X. The details 
of coannihilation has been studied extensively within 
SUSY (Edsjo and GondolQ, [1993) , and pubhcly available 
numerical codes include coann ihilations with all SUSY 
particles (jGondolo et al. I. I2004D . 



would require a very small fraction of the inflaton energy 
density to be transferred to the dark matter particle X, 
otherwise the universe would be dominated by the dark 
matter particles. 

For a singlet hidden sector dark matter, it is really a 
challenge not to overproduce them directly from the decay 
of the inflaton. If the inflaton sector belongs to the hid- 
den sector, then it is natural to have inflaton couplings to 
such hidden sector dark matter field. There are three pos- 
sible ways to obtain a small fraction of px{t*) / Pinfit*) 
in order to match the current observations. 



3. Non-thermal relics 

The dark matter particle, X, can also be created in 
an out of equilibrium condition, i.e. X must not have 
been equilibrium when it froze out. A sufficient condi- 
tion for non-equilibrium is that the annihilation rate (per 
particle) must be smaller than the expansion rate of the 
universe: nx{<Jv) < H. 

Let us assume that X were non-relativistic at the time 
of production and they were never in local thermody- 
namical equilibrium. The largest dark matter density 
will thus be determined by the largest freeze out temper- 
ature, which can be attainable in the universe. Assuming 
this to be the reheat temperature, and the universe 
follows a radiation domina tion, then the rat ios of energy 
densities will be given by (iKolb et aLl . [l998l) : 



Px{t(i) _ pxjtR) 
Prita) Pritu) 



Tr 
To 



(210) 



where Tq is the present temperature and corresponds 
to the present time, pj is the energy density in radiation, 
and px = nixnx denotes the energy density in the dark 
matter with the number density nx- If we further as- 
sume that X particles were created at time t = t^ < tn, 
sometime during the coherent oscillations of the inflaton 
and before the completion of reheating, then both the 
X particle energy density and the inflaton energy den- 
sity would redshift approximately at the same rate until 
reheating is completed. Therefore, 



PxjtR) 



Pxjt*) 
iMlH^{t,y 



(211) 



assuming that the inflaton energy density dominated the 
universe. Since, Vlx = px{to)/ pdto), where Pc{to) — 
3H^M^ and Hq = 100 km sec~^ Mpc ~\ then using 



Eq. ((2T0)) . one obtains (|Kolb et a/.l . ll998[ ) 



nxh^ 



Q,rh 



10' 



In 

To 

( 



VlO^GeVy p^nf{t* 



Mx \ nx{t,) 
Mp J 3MpH^{t,) 

\ Px{t*) 



(212) 



(213) 



where Vl^h? w 4.31 x 10^^ is the fraction of critical energy 
density in radiation today, and Tq ^ 2.3 x 10~^^ GeV. 
The above expression tells us that a non-thermal creation 



(a) Gravitational production 

The dark matter can be created from the transition 
of the equation of state of the universe from inflation to 
matter domination or radiation dom ination, due to non- 
adiabatic evolution of the vacuum (iGhung et all l200ll 
Il999bl l2000t iKolb et all Il998l) . The underlying mech- 
anism is similar to the metric fluctuations which seed 
the structure formation, except now the excitations can 
create massive particles. The gravitational production of 
dark matter is universal, and it can occur even if the dark 
matter coupling to the infiaton is vanishingly small. 

Let us consider a simple action for X field with a metric 
ds^ = dt'^ - a'^{t)dy? 
conformal time. 



■(77) [d?7^ — dx^] , where 77 is a 



S ■ 



d^x X 



y2 



{VXf 



m\X^ - ^RX' 



(214) 

where R is the Ricci scalar. Let us expand the X field in 
terms of creation and annihilation operators which obey: 

[afei'oL] = (^^^Hki -ka), and 



X 



(27r)3/2a(?7) 



akUk{vi)e 



zk-x 



4ufc(?/)e' 



(215) 

where the mode functions obey the identity UkU^ — 
u'f.ul. = i, and prime denotes derivative w.r.t. 77. The 
mode equation is given by: 

<{V) + [k^ + + (6^ - l)a"/a]uk{v) = 0, (216) 

The parameter ^ = 1/6, for conformal and ^ = for 
minimal coupling. Here we will consider the conformal 
coupling for simplicity. The number density of X parti- 
cles can be estimated by a Bogoliubov transformation: 



(217) 



where rjo = —00, and 771 — +0 0. The energy densi ty of 
produced particles is given by (jChung et aLl . ll999bD : 



f°° dk ~ 

Px{vi)=mxHlf Cairn))-' 7^k^\Pk\\ (218) 



27r2 



where the number operator is defined at 771. Assum- 
ing that the transition from inflation-radiation or matter 
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domination is smooth, the largest energy density can be 
obtained if nix/Hinf ~ 1. If 0.04 < mx/Hinf < 2. If 
Hinf ^ ~ lO^^GeV and is the mass of the infla- 
ton, then X particles produced gravitationally can match 
the density today of the order of the critical density pro- 
vided they are long lived. Such super heavy massive dark 
matter particle X is known as Wimpzillas! 



(b) Direct decay of the inflaton 



The dark matter can also be created from direct infla- 



ton decay if mx < iti^/2, with a rate Tx ~ h\m^i 
where hx is the interaction strength. The total inflaton 
decay rate is given by a/1/3 T^/Mp, while the in- 

flaton number density at the time of decay is given by 
ri0 ^ ^r/'^0- This constrains the overall coupling to 



hx < 327r 



1 Tr 10 



3 Mp mx 



(219) 



where mx is in units of GeV. This is required due to the 
fact that the produced X must not overdose the universe 
which, for Qx < 0.22 and Ho = 70 km sec"^Mpc"\ 
reads 



Tlx /n~f < 4 X 10 ^mx^ 



(220) 



when mx is expressed in units of GeV ( Allahverdi et all 
I2OO2I ). It is evident from the overclosure bound that hx 
needs to be very small. 



(c) Creation during reheating 

If the process of reheating is slow and not instanta- 
neous, then it is possible to create WIMP from the am- 
bient plasma which is in the process of acquiring ther- 
malization via scatterings. The Boltzmann equations 
for inflaton energy density, p0, radiation energy den- 
sity j_22-j_and_darkjnAttCT px a re given 



sity . Pr, ana aar K: matt er energy aensity px a re 
by (|Chung et aZl . Il999al: iKolb and Turned Il988[) : 



P0 -I- 3Hp^ + r^p0 = 
PR + AH PR - T^p^ - 



mx 



Px + 'iHpx 



mx 



iplf 



, 



= 

(221) 



where dot denotes time derivative, and thermal averaged 
cross section is given by: (cr|w|). The equilibrium energy 
density for the X particles, p^, is determined by the 
radiation teniperature, T = (30pff /•Tr^qa,)^/'^ . Following 
(|Chung et a/.l . ll999al : lKolb and Turneilll988[ ). it is useful 
to introduce two dimensionless constants, and ax, 

{(j\v\) = axm'^ , and 



defined in terms of T 



a^m^ 
R = prO^ ; 



X = PXTTI 



X 



^"3. With 



these parameters the Boltzmann equations are: 



$' = -ci 
R' = 
X' = 



^2 



ci 



^/<^x + R 

^-2 



C3 



^/^x + R 



R 

$ -HC2 

ix^- 



^/<^x^R 



[X^ - xl 



(222) 



where x = am^, prime denotes d/dx, and the constants 
Ci, C2, and C3 are given by 



Cl 



3 Mp 
, 

OTT m^f, 



C2 = Cl - 



mx 



C3 = C2- 



mx 



The equilibrium value of X is given in terms of the 
temperature T and Eq. (|193p : X^q — nxmx = 
{m\/ml){\/2^f/^x\T/Mxf'^eM-Mx/T). 

It is straightforward to solve the system of equa- 
tions in Eq. (|222l) with initial conditions at a; = 
XI, R{xi) = X{xi) = 0, and ^{xj) = $/ = 
(3/87r)(Mp/m^)(7?|/m^)a;J. At early time solution for 
R can be easily obtained: 



R ~ 0.4c 



(x^/^-xT) 



5/2\ ^1/2 
x/ $/ 



{H » 



(223) 



By maximizing the above equation, which is obtained 
at x/xi = (8/3)^/^ = 1.48, the largest temperature of 
the ambient plasma can be even larg er than the reheat 
temperature (jKolb and Turneil . Il988| ) 



Tmax/Tr - 0.77 (9/57r35, 



(224) 



From Eq. (|223p when x/xj > 1, T scales as a^^/^, which 
implies that entropy is created in the early-time regime. 
For t he choices oi m^, a d,, 5*, and ax, and Vlxh? = 
0.22 l|Chung et aLl . ll999aD : 



200 

9* 



^'V^V (225) 
V mx J 



(d) Non-perturbative creation of darl< matter 

The dark matter can be created non-perturbatively 
during the coherent oscillations of the inflaton. 

a. Superheavy dark matter during preheating: If the 

dark matter couples to the inflaton directly, then it is 
more efficient to excite them from the coherent oscilla- 
tions of the inflaton during preheating. One of the most 
interesting applications of preheating is the copious pro- 
duction of particles which have a mass greater than the 
inflaton mass m^. Such processes are impossible in per- 
turbation theory and in the theory of narrow parametric 
resonance. 

Following Eq. (|97p . let us suppose that the dark mat- 
ter X is coupled to the inflaton with an interaction term: 
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{\/2)g'^X'^4P' . During the broad resonance regime, as we 
have discussed in section Hl.F.li superheavy X-particles 
with mass mx 3> can be produced. The momentum 
dependent frequency uJk{t) violates the adiabatic condi- 
tion of time dependent vacuum, see Eq. (jl03p , when 



(226) 



where <f> is the amplitude of the inflaton oscillations. 
The maximal range of momenta for which particle pro- 
duction occurs corresponds to 0(t) = 0*, where 0* ~ 



1/2y (r7i00)/.g. The maximal value of momentum for 
particles produced at that epoch can be estimated by 
^max + '^x ~ The resonance becomes effi- 

cient for gm^(t) ^ 4rn^. Thus, the inflaton oscillations 
may lead to a copious production of superheavy particles 
with mx ^ if the amplitude of the field (p is large 
enough, gcf) > 4m|- / (jKofman et all ll997|) . 

b. Dark matter from the fragmentation of a scalar 
condensate: Let us assume that coherent oscillations of 
a scalar condensate in a potential U{4>) has a frequency 
which is large compared to the expansion rate of the uni- 
verse. The equation of state is obtained by averaging, 
p/ p = / p) — 1, over one oscillation cycle T. The re- 

sult is: p = (7 - l)p, where 7 = (2/T) /J^ (1 - U{(j))/p) dt 
(iTurned . [19831) . l^'or the case U ^ rn^cjp', one finds 7 = 1, 
so that one effectively obtains the usual case of pressure- 
less, non-relativistic cold matter. When the motion of 
the condensate is not simply oscillatory, such as in the 
case of a rotating trajectory with a phase, one can gener- 
alize the above calculation by integrating over the orbit 
of the condensate. Let us consider the potential 



(227) 



one finds that 7 = (1 -|- a;)/(l -I- a;/2) , -p = x/{2 + x). 
There arises a negative pressure whenever x < 0. This is 
a sign of an instability of the scalar field under arbitrarily 
small perturbations. The quantum fluctuations in the 
condensate grow according to when effective mass of the 
scalar field is much larger than the expansion 



4 = -K\^5^ 



(228) 



li K — 2x < quantum fluctuations of the condensate 
at the scale, A = 27r/|k|, will grow exponentially in time: 



5(j)^{t) = J0(O)exp (-Kk^t) 



(229) 



In reality the onset of non-linearity sets the scale at 
which the spatial coherence of the field can no longer 
be maintained and the condensate fragments. The 
initial fluctuations in the condensate owes to the in- 
flationary perturbations. If the condensate carries a 
global charge, due to charge conservation the energy- 
to-charge ratio changes as the the condensate frag- 
ments. This is what happens in the case of MSSM. 



The AD condensate which was responsible for generating 
the ba ryon asymmetry at the first in st ance could frag- 
ment ( Kasuva and Kawasakil l2000allbl : iKusenko et al\ . 
I2OO9I ). The ground state of these fragmented lumps is 
a non-t( 3pological soliton with a f i xed charge, called the 
Q-ba ll (|ColemanL Il985l : iKusenkol Il997hl : iLee and Panel 
I1992I ). In gauge mediated SUSY breaking scenarios the 
Q-ball s can absolutel y stabl e and can be a candidate for 
CDM (iKusenko et ai . 1998: iKusenko and Shaposhnikovl 
119981: IKusenko anTshoemaked. 120091 ). The Q-balls can 
also be formed fr o m the fragmentation of the infla- 
ton ( Engvist et all . l2002al fb[). The slow surface evapo- 
ration of a Q-ball will also create SUSY LSP, which we 
will discuss below. 



C. Candidates 

In this subsection we will discuss some of the dark mat- 
ter candidates which are well motivated, and the chal- 
lenges they face. 



1. Primordial blackholes 

The primordial b lackholes (PBH) c a n be created in 
the e arly universe ( Carr and HawkinS Il974l : iHawkinj . 
Il97l[ ) , and they can survive the age of the universe with 
a typical lifetime of an evaporating blackhole which is 
given by: 



1017 sec 



M 



10^° grams 



(230) 



If the initial mass M w 10 ^''^ g, the blackhole will be evap- 
orating now, for heavier blackholes the Hawking evap- 
oration is negligible, and they can be a CDM candi- 
date. When M ~ 10^ g, the blackholes would decay 
at the time of BBN. The PBHs are formed from the 
collapse of order one perturbations, S = Sp/p ~ C(l), 
inside the Hubble patch. The detailed numerical simu- 
lation suggest Sr. ^0.7 ( Jedamzik and Niemeveii Il999l : 
iNiemever and Jedamzikl Il998l Il999[ ). In a radiation 
dominated epoch the blackhole mass is bounded by the 
Hubble mass 



H 



10% ( i^^) , 



(231) 



where T is temperature of the thermal bath during radi- 
ation. In spite of novelty in this idea, the detailed calcu- 
lations suggest that it is hard to form primordial black- 
holes just from the collapse of sub-Hubble over densed 
regions - one requires more po wer on small scales n < 
1.25 - 1.30 in Vr(k) = Ak""-^ (ICarr and Lidseyl . Il993l: 
iDrees and ErfanI I2OIII: iGreen and Liddld . Il997l ). while 
the CMB data points towards n ^ 0.96. It was shown in 
a hy drodynamical simulation ( Jedamzik and Niemeveii . 
Il999( ) that primordial blackholes can also be produced 
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in a first order phase tran sition, and during preheat- 
ing (iGreen and Mahkl . l200lD . 

The abundance of PBH contains many uncertainties, 
as the details of the initial gravitational collapse and 
the initial number density of npg^ depends on many 
physical circumstances. These uncertainties can how- 
ever be encoded in terms of the ratio determined by 
the initial time, U; PpBni^i) / Pi^i) = Mnp^^^U)/ p{ti) = 
AMupgjj /3Tis{Ti), by assuming p = 'isT/4, 



scale (|Nakamura et all I2OIOI : iRaffeltl . 120081: ISikivid . 
l2000f ) must be large. 



fa/N > 6 X 10^ GeV , 



(235) 



and the axion mass must be very small, nia < 0.01 eV. 
The cosmological constraints on /a > 2 x lO*" GeV 
(toq < 0.3 eV) arises from BBN, accounting for the 
current bound on the relativistic species, N'^^^ 



3.lt^;2 (jlocco et al ]. l2009t ). Another interesting bound 



P'{M) w 7.98 X 10^29 ( ^ 



3/2 



3ff (^0) 



1 Gpc" 



(232) 



where corresponds to present time. In t erms of this 
fraction 0' , the PBH ab undance is given by (ICarilll975t 
IGreen and Liddld . I1999D 



^ ^ o a ^ ' 



0.5 



13' (M) 
1.15 X 10- 



The value of I3'{M) can be constrained from 
^PBH < ^CDM, which yields /3'(M) < 2.04 x 
10-i8(r2ci3M/0.25)(M/10i5 g)i/2, for mass M > lO^^ g. 
A tighter constrain on /?' arises from a range of astro- 
physical observations, such as BBN, CMB anisotropy, 

■,15 



arises from isocurvature perturbations from CMB. At 
best one can allow less than 10% of the total pertur- 
bations to arise from sources other than the inflaton 
fluctuations-the axions being massless during inflation 
can account for such fluctuations, which limits fa > 
10^^ — 1 0^^ GeV, however, it depends on the scale of 
inflation (jBeltran al\ . 120071 : ISteffenl . l2009f ) . 

The axion life time depends on the axion-photon 
interaction, which gives a long life time compared 
to the age of the universe, i.e., Tq = F" 
64^/(.g2^,m3) ~ 10^°( f^A^-y iQi" GeV)^ 
mu/rud - 0.56 (Ko lb and Turneij . Il988l ). 

Axion is massless for T > 1 GeV > Aqcd and 
it acquires mass only through instanton effects for 
T < Aqcd- For fa/N < 3 x 10^ GeV (correspond- 
ing to TUa > 0.2 eV), the axion is a thermal relic 



>77 
S 



for 



and 7-ray backgrounds for M < 10^^ g, th e bound ^^^^^ decouples after the quark-hadron transition, Tf < 



weakens for larger mass blackholes ( Carr et al. 



2. Axions 

The ax ions were introduced to s olve the strong CP 
problem ( Peccei and Quinnl . Il977al lbh which requires a 
new global chiral symmetry U(1)pq that is broken spon- 
taneously at_rtwPeccei^uini^^ /„ (for reviews , 
see (iKiml . Il987l: IRaffeltl . Il990l: ISikivid . 120081: iTurneil . 
119901) ). The corresp onding ps e udo-Nambu- Go l dston e 
boson is the axion a (|Weinbersj . Il978l: IWilczekl . Il978[) . 
which couples to the gluons 



o n — 
fa/N 32^2 ^M-^ 



(234) 



/ 

150 MeV. The axion is kept in thermal equilibrium with 
TTTT -f-)- TTa. The relic thermal abundance is given by 
^ah"^ = 0.077 (10/g*(r/)) {ma/10 eV), where 5* denotes 
the number of effectively massless degrees of freedom. 

In an opposite limit, when fa/N is very large, the ax- 
ions are never in thermal equilibrium, and in particular 
when Tfi < fa the PQ symmetry is never restored. The 
main production mechanism is due to the coherent os- 
cillations of the axion due to the initial misalignment 
angle Qi of the axion. At T ~ Aqcd, the axion ob- 
tains a temperature dependent effective mass and oscil- 
late coherently around its minimum when ma{Tosc) — 
H{T). These oscillations of the axion condensate be- 
haves as cold dark matter ( Abbott and Sikivid . Il983l : 
iDine and Fischleij . ll983l : lPreskill et al\ . \l98^ with a relic 



density that is governed by the initial misahgnme nt angle 



where N is the color anomaly of the PQ symmetry de- 
pends on the interactions. This interaction term compen- 
sates the vacuum contribution in the QCD Lagrangian 
£e = 9(g2/327r2)G ^^G°^'-, in a way that 9^6-1- 
Arg (detM) < 10"^ (jNakamura et aLl . l2010[) . in order to 
match the electric dipole moment of the neutron. The 
dynamical solution yields when (a) = —Qfa/N, at which 
the effective potential for the axion has its minimum. 

The axion can interact via heavy quark while all other 
SM field s do not carry a ny PQ charge, i n which case 
N = I (jKiml . [l979: Shif man et all Il980[ ). The axion 
can directly couple to the SM, and at the lowest order it 
will induce no n-renormalizable coupling wit h the gluons, 
where TV = 6 (|Dine all Il98ll : IZhitnitskvl . [1980). 

The axion searches, various astrophysical and 
cosmological observations suggest that the PQ 



-TT < 8i < TT ( Beltran et 



f th e mitial misal i gnme nt 
al ]. l2007l : ISikivi3 . l2008f l: 



Uah' ^ 0.15 e fm e? (y^I^) (236) 

with f = 0(1) parametrizing theoretical uncertainties 
related to details of the quark-hadron transition and 
/(Of) accounting for anharmonicity in Qi - /(Of) — > 1 
for Of -> 0. For lO^^GeV < fa/N < lO^^ Qgy, this 
"misalignment mechanism" can provide the correct dark 
matter abundance. 



D. SUSY WIMP 

The most general gauge invariant and renormalizable 
superpotential would also include baryon number B or 
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lepton number L violating terms, with each violating by 
one unit: Wal=i = ^X'^''LiLjek+X"^''L,Qjdk+fi"L,Hf, 
and Wab=i — \\"^^^Uidjdk^ where i — 1,2,3 repre- 
sents the family indices. The chiral supermultiplets carry 
baryon number assignments B — +1/3 for Qi, B — —l/Z 
for Ui,di, and B = for all others. The total lepton num- 
ber assignments are L — +1 for Li, L = —\ for e^, and 
L = for all the others. Unless A' and A" terms are very 
much suppressed, one would obtain rapid proton decay 
which violates both B and L by one unit. 

There exists a discrete Z2 symmetry, which can forbid 
baryon and lepton n umber violating terms, known as im- 
parity ( Favetl . ]T979l) . For each particle: 
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with Pr = +1 for the SM particles and the Higgs bosons, 
while Pfj = —1 for all the sleptons, squarks, gauginos, 
and Higgsinos. Here s is spin of the particle. Besides for- 
bidding B and L violation from the renormalizable inter- 
actions, i?-parity has interesting phenomenological and 
cosmological consequences. The lightest sparticle with 
Pa — —1, the LSP, must be absolutely stable. If elec- 
trically neutral, the LSP is a natura l candidate for dar k 
matter (jPimopoulos and" Hali Il988t lEllis et aLl . ri984aD . 
The advantage here is that their cross sections are gov- 
erned by the SM gauge group - and therefore the dark 
matter paradigm is embedded within a visible sector. 
However, there are some exceptions which we will dis- 
cuss first. 



1. Gravitino 

The gravitino is a spin-3/2 supersymmetric partner 
of the graviton, which is coupled to all the sectors 
universally, e.g. visible and hidden sectors, with a 
Planck suppressed interaction. In this respect gravitino 
is truly a singlet dark matter candidate, if it happens 
to be the LSP. The gravitino mass can vary (^3/2 ^ 
0(100) GeV - 0(1) KeV) depending on the details of 
the SUSY breaking schemes, for instance, in gauge and 

fravity- mediated SUSY br eaking scenarios ( Dine et al\ . 
996a, 1995a; Gi udice and R attazzi. 199^. Indeed, with- 
out considering the SUSY breaking mechanisms and the 
SUSY breaking scale, we can treat the gravitino mass as 
a free parameter. 



Production: Gravitinos with both the helicities can 
be produced from a thermal bath, they are never in ther- 
mal equilibrium. They are produced mainly through 
the scatterings-within MSSM there are many scatter- 
ing channels which include fermion, sfermion, gauge and 
gaugino quanta all of which have a cross-section c>c 1/Mp. 
The thermal abundance is given by (up to a logarith- 
mic correction) for = 228.75 as i n the case o f MSSM: 



{Tr/10^° GeV) X 10-12 (lElUs et al\ . ll984bD . 

and y±i/2 ~ [1 + Af|(Tfl)/12m2/2](^fl/10"' GeV) x 



±3/2 




10' 



10= lO" 10^ 

Gravitino Mass (GeV) 



FIG. 5: The bound on reheat temperature Tr with 
respect to an unstable gravitino mass 7713/2, where 
neutralino is the LSP with a mass 117 GeV ( indicated 
by the shaded light-orange region in which 77x3/2 < m.^o 
for (7771/2, mo) = (300, 141) GeV, Aq = 0, tan/3 = 30.' 
The thermal relic density is given by: Vt^oh^ = 0.111. 
Above the dotted hne labeled as r^LSP, the Xi density 
from decays of thermally produced gravitinos e xceeds 
0.118. see (|Kawasaki et"all . l2008( ). 
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10^12 l|Bolz et aLl . l200lD . where Mg is the gluino mass. 

Note that for Mg < 7773/2 both the helicity states have 
essentially the same abundance, while for Mg ^ 7773/2 
production of helicity ±1/2 states is enhanced due to 
their Goldstino nature. The net abundance 173/2 /i^ can 
be approximated by the convenient expression for the 
universal gaugi no masses Mi, 2. 3 = ?77i /2 at Mqut and 
TO3/2 < Mi,2,3 ijPradler and Steffenl . [20081) 



1^3/2/72 



0.32fi^)f-"^^/^ 
V 7773/2 /VI 



■(- 



Tr 



TeV/ VlO^GeV/ 
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Thermally produced gravitinos have a negligible free- 
streaming velocity today. However, the gravitinos cre- 
ated from decays can be warm or hot dark matter. 

Besides thermal production, gravitino can be pro- 
duced non-thermally from the decay of the NLSP 
(next-to- lightest SUSY particle). Obviously different 
NLSP's gi ve slightly different abundances. For sneutrino , 
see Refs. jArina and Fornengol 120071: lEllis et all l200g : 



Feng et al. 



2004) 



fo r stop NLSP, see (iBerger et al. 
20081 : iDiaz-Cruz et aLl . [20071: iKang et a/.l. 120081 ) A sim- 



ple a pproximation yields ( Pospelov et aLl . [20081 : ISteffenl . 
120061) 



0.7 X 10' 



mr 
'1 

1 TeV 
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where the total slepton number density is given by assum- 
ing an equal number density of positively and negatively 
charged slepton's. The NLSP's can have a long lifetime 
Tr . For a slepton NLSP, one finds in the limit mi — >■ 0, 



487r7772/2^^P 

777^ 
h 



1 - 



'3/2 



9 

777= 



(240) 
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which holds not only for a charged slepton NLSP but 
also for the sneutrino NLSP. Similar expressions for the 



lifetimes for the neutralino NLSP (iFeng et ali. 1200 
and the stop NLSP can be found in toiaz-Cruz et ali 
I2OO7I ). If the NLSP decays into the gravitino LSP 
after BBN, the SM particles emitted in addition to 
the gravitino can affect the abundances of the pri- 
mordial light elements. Also the presence of a long- 
lived negatively charged particle, champ, i.e. 
can l ead to bound s t ates t hat catalyze BBN reac - 
tions (|De Ruiula et all . Il990l : iDimopoulos et all Il990[ ). 
The new acceptable limits on champs: 100{qx/e)^ < 
mx < I0^{qx/e)'^ Te V , virtually ruled o ut an y low 
scale SUSY champs (IChuzhov and Kolbl . I2009D . It 
was suggested in Ref (iPospelovl 1200/1) that bound- 
state formation of champ with He'^ can lead to a sub- 
stantial production of primordial Li ^, which puts the 
constraint 77- 



< 5 X 10^ sec. (iBird et all l200g : 



- • — -*t — ^ I I — I I "J I ' ' 

Hamaguchi et all I2007t IPospelovl I2007t IPospeloy et al. 



20081 : iPradler and Steffenl . l2008t iTakavamal l2008l)" 



Uncertainties: The main uncertainties on gravitino 
abundance arise from the hidden sectors. If the in- 
flaton sector is embedded within a hidden sector then 
there are many more sources of gravitino production- 
the inflaton could deca y directly into gravitino during 
reheating or preheating (Frey et aLV2006VGiudice et all, 
1999b; Kallosh ei at, 200 0a;b; Maroto and Mazumdar, 
2OOOI : iNilles efai] . l2001at iNilles and Pelosoi l200"il ). the 



inflaton couplings to other hidden sectors can similarly 
excite gravitinos giving rise to large uncertainties in their 
total abundance. 



r^In flaton 7 2 
"3/2 ^ 



Hidden j 2 
3/2 " 
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All these contributions can easily overproduce gravitinos, 
i.e. $13/2/1^ ^ 1, especially the last two sectors are largely 
unconstrained by particle physics. These uncertainties 
can be minimized if the last phase of inflation occurs 
within MSSM, as discussed in Sect. ??, then the only 
predominant source of gravitino production arises from 
the decay of the MSSM inflaton and from the MSSM 
thermal bath. 

Another solution has been put forward ~ for high scale 
and hidden sector models of inflation - since the flat di- 
rections of MSSM can be displaced from their minimum 
during inflation, the flat direction VEV at early times 
would generate time dependent masses to the MSSM 
fields which are coupled to the flat direction. As a result, 
the inflaton might not even d ecay into all the MSSM d.o.f 
due t o kinematical b locking (jAUahverdi and Mazumdar . 
I2OO5I l2006bl l2007bD . Furthermore, the flat direction 
VEV also generates masses to gauge bosons and gauginos 
which participate in scatterings, therefore delaying the 
actual thermalization process and lowering the reheating 
temperature, i.e. Tr. Both these effects address the ther- 
mal gravitino overpro duction problem without any nee d 
of extra assumptions (jAUahverdi and Mazumdail 2005f) . 



Unstable gravitino: An unstable gravitino decays 
to particle-sparticle pairs, and its decay rate is given 
by r3/2 ^ ml/^/4M^. If 7713/2 < 50 TeV, the graviti- 
nos decay during or after BBN, which can ruin its suc- 
cessful pre dictions for th e prirn o rdial abundance of ligh t 
elements (jCvburt et all l2003t iKawasaki eTdl . l2005t) . 
If the gravitinos decay radiatively, the most stringent 
bound, (713/2/5) < 10~^^ — 10~ ^ ^, ar ises for 7713/2 — 
100 GeV - 1 TeV llCvburt et all I2003D . On the other 
hand, much stronger bounds are derived if the graviti- 
nos mainly decay through the hadronic modes. In par- 
ticular, for a hadronic branching ratio ~ 1, and in the 

same ma ss range, (713/^ / s) < 10~^^ — 10~^^ will be 

required ( Kawasaki et all l2005l . |2008[ ). This puts con- 
straint on reheat temperature of the universe, i.e. Tr, at 
which these unstable gravitinos are produced, see Fig. [5] 
An intriguing possibility arises if R-parity is broken. The 
gravitino LSP with 7^3/2 ^ 1 GeV can still be a long- 
lived cold dark matter, while evading the bounds on 
from Fig. [Sj The gravitino in this case cannot decay into 
hadrons which is kinematically suppressed, and the three- 
body dec ay life time is typic ally larger than the age of the 
universe (|Kohri et aL l. l2009l) . For a GeV mass gravitino, 
the present day free-streaming velocity is < 10~^ km/s, 
which corresponds to that of a cold dark matter. 



Detection: The direct detection of gravitino will 
be impossible at the LHC, their production will be 
extremely suppressed. If the NLSP is long lived 
(quasi-stable) as in the case of stau then they would 
penetrat e the collider d etecto r in a way sirn il ar to 
muons (iDrees and Tatal . Il990l: iFeng and Moroil . Il998l: 
iNisati et al . 1997t) . If the produced staus are slow, 
then from the associated highly ionizing tracks and 
time of flight measure ments one can determine their 
mass ((iniritall,[200i). This might give some indirect 
handle on gravitinos. 



2. Axino 

The axino, a, is a superpartner of the axion, is an- 
other example of a gauge sing l et dark matter candi- 
date dKiml, 119841: iKim and Niiiesl . Il984t INilles and Rabvl . 
Il982( ). It interacts extremely weakly since its c ouplings 
are s u ppres s ed by th e PQ scale fa ^ 10* GeV (jRaffeltl . 
I2OO7I . I2OO8I: ISikivid . l2008l ). and its mass can range 
from eV to GeV. In the hadronic axion model ( Kiml . 
I1979I : IShifman et aZI . [19801) in a SUSY setting, the ax- 
ino couples to MSSM field indirectly via loops of heavy 
(s)quarks. Typically a decouples early at a tempera- 
ture Tf > 10^ GeV, below this temperature, they are 
mainly created from scatterings of MSSM fields in an 
kine tic equilibrium. The the r mal abundance is given 
by (|Brandenburg and Steffe^ . |2004 IChoi et all 120081: 
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FIG. 6: The plot shows the reheatmg temperature 
with respect to the axion mass for fa/N — 10^^ GeV. 
The gray band indicates n^h^ = 0.113t^ |^i^. Thermally 

produced axinos can be c lassified as hot, warm, and 
cold dark matter (|Covi et a t. 200 1,) as indicated in the 
plot. The plot is taken f rom (iBrandenburg and Steffenl . 
l2004f l. 
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FIG. 7: On the left hand panel the co-annihilation band 

is shown for (mp — plane for tan = 40. The 

constraints are shown in Figure, see (iNath et aLl . l2010h . 

On the right hand panel the overlapping contours 
between MSSM inflation and neutralino dark matter for 
different values of (n^. 6h) within 95%c.l. are shown for 
the same parameter region l|Allahverdi et al. l. l2007aD . 



ICovi et a/.l . [200l 

nah^ ~ 5.5gf(Tj^)ln 

(0.1 Gev) 



1.108 \ /lO" GeV\ 
9s{Tb.)) V fa/N J 
Tr 
104 GeV 
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with the axion-model-dependent color anomaly N of the 
PQ symmetry breaking scale fa- Thermally p roduced ax- 
inos c an be hot, warm, and cold dark matter ( Covi et all . 
I2OOID as shown in Fig. H 

Non-thermal production of a has many uncertainties. 
The a can be created from the decay of the NLSP, direct 
decay from the inflaton or moduli, or any other hidden 
sector. The expression for the final abundance is similar 
to that of Eq. (|24ip . where 3/2 — !• a. In this regard a 
faces similar challenges as that of a gravitino dark matter. 

The a LSP is inaccessible in direct /indirect dark 
matter searches if R-parity is conserved. Their direct 
production at collider is strongly suppressed. Never- 
theless, quasi-stable f's could ap pear in collider de- 
tectors (and neut r ino t elescopes ( Ahlers et all l2006t 
lAlbuauergue et all . l2007t) ) as a possible signature of the 
SLSP. 



3. Neutralino 

In the MSSM the binos B (superpartner of B), winos 
W (superpartner of W) and Higgsinos (i?° and H^) 
mix into 4 Majorana fermion eigenstates, called neu- 
tralinos with 4 mass eigenstates: XiJ X2' Xsj xI^ oi'" 
dered with increasing mass. The LSP is thus denoted 
by Xi = NnB + N12W + NuHl>^ + NuHl The gaug- 



fn = Nf^ + N^^, are determined by the mixing ma- 
trix, A^, which diagonalizes the neutralino mass ma- 



trix (jJungman et ad Il996l: Ikane et aZ.I . Il994l) . 

The x5's were in thermal equilibrium for primordial 
temperatures oi T > Tf ^ m^o/20. At Tf, the annihi- 
lation rate of the (by then) non-relativistic Xi's becomes 
smaller than the Hubble rate so that they decouple from 
the thermal plasma, see Sect. I IV. B. II 

It is easy to work with a limited set of parameters, the 
mSUGRA model is a simple model which contains only 
five parameters: 



Too, ™i/2, tan/3 and sign{ii). 



(243) 



Too is the universal scalar soft breaking parameter, TO1/2 
is the universal gaugino mass, Aq is the universal cu- 
bic soft breaking mass, measures at Mgut, and tan/3 — 
{H^)/{H^) at the electroweak scale. 

The model parameters are already constrained by dif- 
ferent experimental results, (i) the ligh t Higgs mass 
bound of M/jO > 114 GeV from LEP (jBarate et al\ . 
I2OO3D . (ii) the h S7 branching ratio bound of 1.8 x 
10~4 < B{B — > As7) < 4.5x10"'' (we assume here a rela- 
tively broad range, since there are theoretical errors in ex- 
tract ing the branching ratio from the data) (jAlam et al\ . 
llQQSl ). (iii) the 2cr bound on the d ark matter re lic den- 
sity: 0.095 < r^cDM/i^ < 0.129 fiKomatsu et all 12009). 
(iv) the bound on the lightest chargino mass of M^± > 



104 GeV from LEP (|Nakamura et al\ . \201(i ) and (v)'the 
muon magnetic moment anomaly a^, where one gets a 
3.3ct deviation from the SM from the experimental re- 
sult ([Bennett et al\A2004 \. 

The allowed mSUGRA parameter space, at present, 
has mostly three distinct regions: (i) the stau-neutralino 
ino fraction, fc = N-^i + ^12' and Higgsino fraction, (fi — % J ), coannihilation region where Xg the lightest 
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SUSY particle (LSP), (ii) the Xo having a dominant Hig- 
gsino component (focus point) and (iii) the scalar Higgs 
(A", H°) annihilation funnel (2M^i ~ M^o_jyo). These 
three regions have been selected out by the CDM con- 
straint. There stills exists a bulk region where none of 
these above properties is observed, but this region is now 
very small due to the existence of other experimental 
bounds. The allowed parameter space for the neutralino 
dark matter (blue region) for tan(/3) — 40 is shown in the 
left panel of Fig. H 

Detection: In general the observable signals for SUSY 
at LHC are: n leptons-l- m jets + flT (missing transverse 
energy), where either n or m could be 0. The existence of 
missing energy in the signal will tell us the possibility of 
dark matter candidate. There are SM backgrounds, e.g. 
W and Z bosons decaying to neutrinos providing fq^. 
The clean signal for SUSY would be jets + pY, without 
isolated leptons. One of the key analysis for mSUGRA is 
to measure the M^ff which is the sum of the transverse 
momenta of the four leading jets and the missing trans- 
verse energy: Me// = Pt.I +?'T,2 +Pt,3 +Pt,4 + One 
has to further measure the masses (squarks, sleptons), 
and tan /3, and the mixing matrices which lead to the cal- 
culation of the relic density. One particularly favored pa- 
rameter space is the coannihilation region where the stau 
and the neutralino masses are close for smaller values of 
toq. The mass difference. Am, governs the relic abun- 
dance due to the Boltzmann suppression factor e~^*^/^, 
see section IIV.B.2I Therefore measuring AM directly 
gives handle on measuring the relic abundance at th e 
LHC, see for a detailed discussion in ( Nath et aLl . |2010[ ). 



MSSM inflation and dark matter: After consid- 
ering all these bounds it was found that there exists an 
interesting overlap between the constraints from MSSM 
inflation and the abundance, see the right hand panel 
of Fig. [7] for tan j3 = 40. The constraints on the parame- 
ter space arising from the inflation appearing to be con- 
sistent with the constraints arising from the dark matter 
content of the universe and other experimental results. 

It is also interesting to note that the allowed region 
for Uid2d^ as an MSSM inflaton with a mass m^, re- 
quired by the CMB observations for A = 1, see Fig. [3] 
in Sect. III. E. 21 lies in the stau- neutralino coannihilation 
region which requires smaller values of the SUSY par- 
ticle masses ( Allahyerdi et al\. 2007a l). Similar analysis 



were performed in (jBalazs et al\ . 1200 5. 2004), where the 



authors studied the overlap between MSSM parameters 
for the electroweak baryogenesis and the x5 dark matter 
abundance. 



4. Sneutrino 

The lightest right handed (RH) sneutrino N can be a 
good dark matter candidate when the SM gauge group 
is augmented to SU{2,)c x SU{2)l x U{1)y x U{1)b-l, 
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The solid lines from left to right 



are for = 0.094 and 0.129 respectively. The Z'-ino 
mass is equal to the Bino mass since the new U{1) gauge 
coupling is the same as the hypercharge gauge coup ling. 
The plot is taken from ()Allahverdi et a/.l . ]20071: 



with a superpotential (jAllahverdi et all l2007bl . l2010b[ ) 



W = Wmssm + Wb-l + IjNH^L . 



(244) 



The model introduces new gauge boson Z', two Higgs 
fields H[, H2, and their superpartners, the Yukawa cou- 
pling is denoted by y. The spontaneous breaking of 
U{1)b-l will generate Majorana neutrino masses, or if 
y « 10~^^, then Dirac neutri no masses, or a mixtur e 
of both Dirac and Majorana ( AUahverdi et all l2011at ). 
If the right handed sneutrino is the LSP then it pro- 
vides another compelling candidate which is well moti- 
vated from particle theory an d can be embedd e d with 
least u nknown uncertainties (AUahverd i al. I, l2007bl: 
lArina a nd Forncngo, 2007; Lcc et al, 200^. 

Scatterings via the new U{1) gauge interactions bring 
the RH sneutrino into thermal equilibrium. In order to 
calculate the relic abundance of the RH sneutrino, we 
need to know the masses of the additional gauge boson Z' 
and its SUSY partner Z' , the new Higgsino masses, Higgs 
VEVs which break the new U{1)b~l gauge symmetry, 
the RH sneutrino mass, the new gauge coupling, and the 
charge assignments for the additional U{1). The primary 
diagrams responsible to provide the right amount of relic 
density are mediated by Z' in the ^-channel 

By assuming that the new gauge symmetry is broken 
around TeV in Fig.[8l we show the relic density values for 
smaller masses of sneutrino where the lighter stop mass 
is < 1 TeV. The smaller stop mass will be easily acces- 
sible at the LHC. By varying new gaugino and Higgsino 
masses a nd the ratio of the VE Vs of new Higgs fields, the 
WMAP (|Komatsu et al\ . \200^ allowed values of the rehc 
density, i.e., 0.094 — 0.129 is satisfied for many points. In 
the case of a larger sneutrino mass in this model, the 
correct dark matter abundance can be obtained by anni- 
hilati on via Z' pole ( AUahverdi et all l2007bl : [Lee et all 
I2007D . 

Detection: Since the dark matter candidate, the RH 
sneutrino, interacts with quarks via the Z' boson, it is 
possible to see it via the direct detection experiments. 
The detection cross sections are not small as the inter- 
action diagram involves Z' in the t-channel. The typical 
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cross section is about 2xl0~® pb for a Z' mass around 
2 TeV. It is possible to probe this model in the upcom- 
ing dark matter detection experiments. The signal for 
this scenario at the LHC will contain standard jets plus 
missing energy and jets plus leptons plus missing energy. 
The jets and the leptons will be produced from the cas- 
cade decays of squarks and gluinos into the final state 
containing the sneutrino. 



5. Stable and evaporating Q-ball, LSP dark matter 

The AD condensate fragments to fo r m Q- 
balls, for reviews s e e (iDine and Kusenkol 120041 : 
lEngvist and Mazumdail . l2003l) . and a finite size Q- 
ball has a minimum of energy and it is stable with 
respect to decay into free quanta if C/(0)/(/)2 — min, for 
(/-o >0. 

For a s ufficiently large Q, the energy of a soliton is then 
given by (|Colemanl . ll985l : lKusenkd . ll997btlLee and Panel 
I1992D : E = \vQ\ < 7710 1 Q I, which ensures its stability 
against decay into plane wave solutions with ~ (/Sq 
inside and ~ outside the soliton, where x) 



(fi(x). The value of i' was computed in ( Kusenkol 



1997al lbl). Note that the global U{1) symmetry is thus 



broken inside the soliton by the VEV, however, remains 
unbroken outside the soliton. The most crucial piece is 
the presence of a global U{1) charge of the Q-ball which 
actually prevents it from decaying and makes the soliton 
stable. 

In the context of gauge mediated SUSY 
breaking the AD potential takes the form 

Uiif) ^ mllog(l + \ip\ym l) (iDvah et al\ . 

Il998l iKusenko and Shaposhnikwl Il998[ ). where 
7770 O(TeV), represents the SUSY break- 

ing scale. The profile of the Q-ball is given by 
ip{r) ~ exp(— 77i0r), whe re the energy, radius and the 
VEV of a Q-ba ll goes as ("Enqv ist and Mazumdail l2003t 
IKusenko and Shaposhnikov. 19981 ) 



)l/4 



7710 



,^0 



Q"" (245) 



This allows for the existence of some entirely sta- 
ble Q-balls with a large baryon number Q ^ B (B- 
balls). Indeed, if the mass of a B-ball is Mb ~ 
(1 TeV) X B3/4^ ^Yien the energy per baryon number 
[Mb/B) - (1 TeV) X jg iggs than 1 GeV for 

B > 10^2. Such large B-balls cannot dissociate into 
protons and neutrons and are entirely stable - thanks 
to the conservation of energy and the baryon num- 
ber. If they were produced in the early universe, 
they would exist at present as a for m of dark mat- 
ter ( Kusenko and Shaposhnikovl 1 19981). There are as- 
trophysical and terrestrial lilni its ( Kusenko et all Il998t 
IKusenko and Shoemakeil l2009l ). and direct searches for 
the Q-balls, which plac es a lower limit on Q > 
10^2 ^Arafune et a/.ll2000l) . 



In the gravity mediated case the _B-balls are not sta- 
ble, but they evaporate via surface evaporation. In this 
process AD condensate can generate the required baryon 
asymmetry and also create dark mater. The appropriate 
candidate will be the udd flat direction, lifted by 77 = 6 
operator, which carries the baryon number and the right 
dark matter abundance, see Eq. (|179|) . 

When a i?-ball decays, for each unit of B produced, 
corresponding to the decay of 3 squarks to quarks, there 
will be at least three units of i?-parity produced, corre- 
sponding to at least 3 LSPs (depending on the nature of 
the cascade produced by the squark decay and the LSP 
mass, more LSP pairs could be produced). Let N^^p > 3 
be the number of LSPs produced per baryon number 
and Jb be the fraction of the total B asymmetry con- 
tained in _B-balls. Then the baryon to dark matter ratio, 
rB = Pr/pdm, and the dark matter abu ndance are given 
by (|Enavist and McDonaldl . Il998l Il999t) . 



rB 



^LSP ~ 3/b 



(246) 



where 7n„ is the nucleon mass and m^gp is the LSP mass. 
It is rather natural to have < 1. 

The LSPs produced in i3-ball decays will collide with 
themselves and with other weakly interacting particles 
in the background and settle locally into a kinetic equi- 
librium. Thermal contact can be maintained until Tf ^ 
777^gp/20, and a rough freeze-out condition for LSPs (if 
they were initially in thermal equilibrium) will be given 



by: 



where (Toff is the LSP an- 



nihilation cross-section and the subscript / refers to the 
freeze-out values. The thermally averaged cross section 
can be written as (aeSv) = a + bT/m^o, where a and 
6 depend on the couplings and the masses of the light 
fermions ( Jungman et aLl [19961) . 

Assuming an efficient LSP production, so that fB = ^, 
one finds for the LSP density for b 



where Uf w {m^oTfy'^ exp[—m^o/Tf]. The LSPs pro- 
duced in i3-ball decays will spread out by a random walk 
with a rate v determined by the collision frequency di- 
vided by a thermal velocity Vth ~ {T /m^o)^/'^ . Since 
the decay is spherically symmetric, it is very likely that 
the LSPs have a Gaussian distribution. In terms of the 
density par ameter f^^o , the neutr alino abundance can be 
written as i Fuiii and Hamaguchi |2003) 



TOy? 10-^GeV2 lOOMeV / 10 \ 



Q,r;0 ~ 0.5 „ , , 1 , . 

^1 lOOGeV {av) Td \g*{Td) 



(247) 



where the decay temperature of the B-ball is given by 



Td « 21 



777^ 



100 GeV 



3/16 / 10 



20 



1/8 



100 



3/16 



GeV. 
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Below this temperature the annihilations of Xi ^-re neg- 
ligible. Similar analysis can be performed for other LSP 
candidates. For example, if the gravitino is an LSP, the 
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gra vitino abundance from the Q-ba lls decay will be given 
by ( Shoemaker and Kusenkol . |2009[ ): 



0.11 



VI GeV/ 



0.02 



(249) 



where Ng = 3 and /{, ~ 1. The above expression is 
valid for temperatures below lO'' GeV. Similar expression 
can be derived for the Q-ball s deca ying into axino dark 
matter ( Roszkowski and Setol . l2007l ). 



E. Detection of WIMP 

The direct detection of WIMP is the cleanest way to 
seek the identity of the dark matter, their detection is 
possible thro ugh elastic collision with the n uclei at terres- 
trial targets (jCoodman and Wittenl . Il985l) . This method 
is especially promising for detecting SUSY WIMP can- 
didates, such as neutralino or sneutrino. There are also 
ways of inferring WIMP in the sky by using the galaxy it- 
self as a detector in the indirect dark matter searches via 
studying the gamma rays, high energy neutrinos, charged 
leptons, proton, anti-proton background from the decay 
or annihilation of the dark matter particles. 



1. Direct detection of WIIVlPs 

Important quantity is the recoil energy deposited by 
the WIMP interaction with the nucleus of mass tojv in 
an elastic collision, Ej. — m^v^^l — cos0)/tojv, where 
is the WIMP nucleus reduced mass, 9 is the scattering 
angle in the dark matter-nucleus center-of-mass frame, 
and V is the velocity relative to the detector, and it is of 
the order of the galactic rotation velocity ~ 200 km/s. 
Typical recoil energies are Er ~ 0(1 — 100) keV. 

The differential rate for WIMP elastic scattering off 
nuclei is given by (iLewin and SmithI , Il996l) 



dR 

dE^. 



Po 
mw 



dv f{v) V 



da 
dE^ 
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where iV^ represents the number of the target nuclei, mx 
is the dark matter mass and po ^ 0.3 GeV/cm'^ is the 
local WIMP density in the galactic halo, v and f{v) are 
the WIMP velocity and velocity distribution function in 
the Earth frame, which we take it to be Maxwellian, and 
da/dEj^, is the WIMP-nucleus differential cross section. 
The velocity w^m = \/ {mNEr/2m'^), and v„iax is the 
escape velocity of the WIMP in the Earth frame, Vesc — 
5Utlt km/s. 

In fact, the Earth velocity with respect to the dark 
matter halo must be written as = vq (1.05-1-0.07 coswi) 
where 1.05 is the galactic velocity of the Sun and 
w = 27r/lyear. The 7% modulation is due to the 



rotat i on of the Ea.rth a round the Sun (jPrukier et al\ . 
119861 : iFreese" et al In the above expression, f{v) 

must be replaced by f{\v — iTeD- There also exists a 



forward-backward a symmetry in a directiona l signa l as 
first pointed out in (jCopi et all Il999t ISpergel Il988[ ). 

For a given momentum transfer q, the differential cross 
section depends on the nuclear form factor 



da 
d^ 



(251) 



where F{q) is a dimensionless form factor such that 
F(0) = 1, in which case ctq corresponds to the total cross- 
section. It is possible to estimate the parameters ao and 
F{q), for example in the case of neutralino WIMP, which 
is a Majorana f ermion therefore it on ly has axial and 
scalar couplings ( Jungman et a/.l . [19961) . 



a. Spin-dependent cross-section: The axial part of 
the neutralino-quark interaction is mediated via Z bo- 
son and squark exchange Cq^ ~ (Xj^^j^X) (97/^759). 
At the level of nutralino-nucleon interaction by con- 
sidering the nucleon matrix element {n\q'^'^j^\n), 
the effective Lagrangian is given by: -CJ^^f = 
2V2Gf a{n) i^l'^lb^) ("■7/i75'^)i where Gf is the Fermi 
constant and a„ is a dimensionless parameter. For a 
nucleus of spin J, with (Sp) and (Sn) being the aver- 
age spins "carried" by protons and neutrons respectively, 
the cross-section at zero momentum transfer is given 



by ( Jungman et aU [19961 ) 



^(<Z = 0) = ^G|A2j(J+l), 
dq"^ TTW-^ 



(252) 



where A = (ap(5p) -I- a„(5„))/J. Additional cor- 
rections to the form factor is required to take 
into account of the non-zero momentum trans- 
fer. There are many experiments which are sensi- 
tive to spin-dependent c ross section wi t h pur e pro- 
ton cou plings. DAMA (iBernabei et "oil . I2004D. PI- 
CAS SO (lArchambault et al\. \2Q0^. KIMS (iLee et al\ . 
l2007t ). Recentlv COUPP (' Behnke et all I2OIII) has set 
the best constraint on spin-dependent cross section down 
to 7 X 10-38 cm2 for a WIMP mass - 30 GeV. 

b. Spin-independent cross-section: The scalar part of 
the neutralino-quark interaction is mediated via Higgses 
and squark exchanges: Cqx = fq{Q<l){XX). To express 
the neutralino-nucleon coupling one needs the nucleon 
matrix element mq{n\qq\n) = rn(„), the effective interac- 
tion has the form: = f(^n) {XX){nn), where /„ con- 
tains the information about hadron physics, and typically 
it is the same for proton or neutron, i.e., /p ~ /«. In the 
case of right handed sneutrino also there exists no spin- 
dependent part as it has no axial-vector interactions, and 
the cross-section is dominated by the spin-independent 
part. One can define a single spin- independent WIMP- 
nucleon cross s ection gp = indep endent of the spin 



of the nucleon ( Jungman et all Il99 



epe 



da 



lF\q) 



[Zfj, + iA-Z)f„]'-^ 
dq'^ TTw^ 
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FIG. 9: Spin- independent elastic WIMP-nucleon 
cross-section a as function of WIMP mass m^^. The new 
XENONIOO limit at 90% CL is shown as the thick 
(blue) line together with the la a nd 2g sensitivity of 
this run (shaded blue band). From ( Aprile et al\ . [201l[) . 



In the spin-independant case, for low nuclear recoil en- 
ergies Fi^q) ^ 1, there is a coherence effect which boosts 
the WIMP-nucleus cross section by a factor A^m'^. As 
a result this technique is better suited to detect the 
WIMP candidate with heavy nucleus targets. Recently 
the most stringent limits on the elastic spin-independent 
WIMP-nucleon cross-section has been g iven by num- 
ber o f experiments, suc h as CDMS-II (Ahmed et all . 



l2010f ). EDEL WEISS-II (lArmengaud et all l201lh and 
XENONIOO (| Aprile et all l201ll ). These limits are 
shown in Fig. [HI The shaded gray area also shows 
the expected region of CMSSM for the WIMP mass 
and the cross section are indic ated at 68% and 
95% CL ()Buchmueller et all l201l[ ). The current re- 
sults also cove rs the 90% CL are as favored by Co- 
GeN T f green) (lAalseth all l201ll) and DAMA (light 
red) (jSavage etall\200^ . 

It should be noted th at CoGeNT (lAalseth et all 
|2011D . and DA MA/Nal (iBernabei et all I2004D and 
DAM A/LIBRA (jBernabei et all l2008l) collaborations 
have observed an annual modulation signal. The com- 
bined results of the latter group stands at greater than 
8(7 statistical significance. The modulation signal phase 
matches well with the ex pected annual s i gnal of WIMPS, 
and subsequent data ( Bernabei et all l2010f ) has in- 
creased the statistical significance of the modulation sig- 
nal. However the annual modulation claim has not been 
verified by any other experiments, especially the null re- 
sults from CDMS, XENONIO, and X ENONIOO data. 
The CDMS data (|Ahmed et all |2010D shows 2 signal 
events with 0.6 ±0.1 events expected as background, but 
due to low statistics this result has not provided sufficient 
evidence for the dark matter. 



2. Indirect detection 

The indirect dark matter detection depends on the 
nature of the WIMP. If the WIMP belongs to the vis- 
ible sector, or if it has some SM interactions, then their 



annihilation or decay would yield to the known parti- 
cle physics sp ectrum, for a recent review on indirect 
detection, see (jPorter et all l201l[ ). However the spec- 
trum would depend on where they are produced, their 
energy deposition, and what are their final states, e.g. 
7, e^, • • • etc. The signal could be a hard spectrum with 
a monochroma tic line if W IMPS annihilate directly into 
photons ( Berg strom et all [l99 8) , or a continuum spec- 
trum if they annihilate into a pair of intermediate parti- 
cles {{q = u,d,c,s,t,b), q, Z, g, W^jl"^). The former 
process is generically suppressed compared to the lat- 
ter. Most of these latter particles, i.e. W, Z, g decay 
into p, p, 7r°, and a tiny fraction of deuterium or anti- 
deuterium D/D. The tt'^s decay to gamma rays, while 
the TT^ decays produce . If the final states of decay or 
annihilations are e^s or fi^s, they dominantly produce 
a hard spectrum, with the fi^ decays into t'^t and 
i/e- If the final states have r^, they produce a softer 
spectrum and a strong neutrino signal. The can also 
decay hadronically to pions and thus can also produce a 
strong 7-ray signal. 

The source spectrum is generically given by 



E 



dE 



■B 
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where / denotes the annihilating final states, each with 
branching fraction Bf s with E being the energy of sec- 
ondary particles. The production rate per annihilation 
of species / is given by dNf/dE. For a decaying dark 
matter {av)/2mx can be replaced by T/Mx, where F is 
the decay rate. 

The fiux of such final states would depend on the an- 
nihilation or decay rate, which in turn would depend on 
the dark matter density cx p^c- Therefore, the natural 
sources to look at in the sky are the nearby galactic cen- 
ters - where there are large astrophysical uncertainties, 
dwarf galaxies - which have small astrophysical back- 
ground, and galactic centers - where the dark matter 
density is very large but distant sources would yield a 
local tiny flux. 

Gamma rays and neutrinos are perhaps the cleanest 
signals if they are produced as a result of WIMP an- 
nihilations or decays as they are undeflected by mag- 
netic fields and effectively indicating the direction to their 
source. The flux is given by the integral of the WIMP 
density-squared along the line of sight from the observer 
to the source, multiplied by the production spectrum 
(b^iE,^p)^Ji^)x<f^iE) 



^j{E,i,)^ 



dNf 



dE Airmj^ 



Bf,s / dsp'iris,^)), (255) 



where "/" denotes the flnal states and the coordinate s 
runs along the line of sight, E is the gamma-ray energy, 
and ip is the elongation angle with respect to the center 
of the source. The astrophysics related term is hidden 
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in (jBergstrom et all Il998l ) 
1 / 1 



- ^ Koba/ashi (1999) o *MS (2002) » FERMI (2010) 

- 4 CAPRICE (2000) ■ ATIC-1,2 (2003) 
-■j^HEAT (2001) X PPB-BETS (2003) 

- o BETS (2001) T HESS (2003-(f9) j + 5^ 



^.5 Kpc \ 0.3 GeV/cm-^ 
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where the integration is in the direction ip along the hne 
£. The above expression is also valid for neutrinos if the 
source is not far away from us. 

On the other hand the charged particles do not 
have the directional sensitivity, they lose it in their 
course of path in a random motion due to the in- 
terstellar magnetic field. The motion of are de- 
flected by the interstellar radiation field by synchro- 
ton radiation in presence of magnetic field. If pro- 
duced at energies > 100 GeV and if their sources are 
within few kiloparsecs, then they can reach the solar sys- 
tem. Furthermore, the cosmic rays from the primary 
and secondary products can also generate 7 rays dur- 
ing their course through the ISM, all these effects can 
be captured numerically in the publically available code 
GALPROP (jStrong and MoskalenkolllQQSHStrong et all . 
|2000LI2004^ ■ for a review see (|Strong et al\ . \200l\} . 

The main challenge is to disentangle the WIMP sig- 
nals from astrophysical backgrounds, the powerful dis- 
criminator is the spectral tilt in the power spectrum. 
It is quite possible to have a significant fraction of 
WIMP annihilation or decay into monoenergetic pho- 
tons, giving rise a distinctive line in the gamma-ray spec- 
trum ( Bergstrom et aLl . ll998l ). but the likely signal would 
be to have a relatively hard continuum spectrum with a 
bump or edge near the WIMP mass that is above the 
astrophysical background. 

In recent years there have been many new experiments 
which have propelled the indirect search for dark mat- 
ter research vigorously. For instance, the ATIC data 
- which shows a signifi cant bump in the electron flux 
around 300 - 800 GeV (|Chang et al\ . \200^ . where con- 
ventional astrophysical sources would have predicted a 
decaying power law spectrum, and the PAMELA - which 
shows a positron fraction wh i ch ha. s a rising slope above 
10 GeV (|Adriani et oil . l2009l[2010l ). but the anti-proton 
flux and the fraction matches well with the expectations 
of the astrophysical origins. The Fermi-LAT collabora- 
tion which has als o produced a i i electron spectrum fronr 
7 GeV to 1 TeV (|Abdo all 120091 : lAckermann et all 
I2OIOD does not conflrm the rising slope of the ATIC, see 
Fig. [ini rather the data match es well with the HESS 
at hi gher end of the spectrum ( Aharonian et all l2008l 
[200i. 

The rise in the e+ fraction in the PAMELA data 
may have its roots in the astrophysical objects such as 
nearby pulsars, Monogem at a distance of d = 290pc, and 
Geminga at a distance of c? = 160pc. Both are nearby 
objects to Earth which can contribute significantly to the 
flux, and can match both the PAMELA and Fermi- 



LAT data sets, see ( Grasso et all [2009). A dark mat- 
ter interpretation of these data sets requires an ad-hoc 
assumption on the leptohilic nature of the WIMP with 




10° 10' 10' 10' 

E (GeV) 

FIG. 10: Combined electron and positron spectrum as 
measured by Fermi-LAT for one ye ar of observations. 
together with other me asurements ( Abdo et aLl . l2009t 
lAckermann et all . I2OIOI) . The systematic errors for the 
measurement are shown by the grey band. The 
systematic uncertainty associated with the absolute 
energy scale is shown by the non-vertical arrow. The 
dashed line shows the background from secondary in 
cosmic rays from GALPROP. 



an annihilation cross-section {av) ^ 10 
than 2 — 3 orders of magnitude enhancement from ther- 
mally produced WIMPS at the TeV scale, se Eq. (1^ . 
Such an enhancement in the cross-section is possible 
at low energies and with low velocities of dark mat- 
ter provided that the range of the interaction respon- 
sible for annihilation is much shorter than the long- 
range attrac tive potential, known as a Sommerfeld en- 
hancement ( Arkaiii-Hamed et all l2009l ). see for a de- 
tailed derivatio n (ISlatveil. 201Cl[) i n the mass range of 
500-900 GeV (|Grasso eVam2009[ ). It is also possible to 
conceive a scenario of non-thermal dark matter produc- 
tion in order to evade the strict bound on cross-section, 



or decaying dark matter scenarios , i.e. (jArvanitaki et all 
120091: iFairbairn and Zupaiil . l2009l) . 

Finally, one would also expect a gamma ray signal from 
the galactic center, the typical signature for a WIMP 
annihilation will be a line at the WIMP mass, due to 
the 27s, or 7Z production channels. The Fermi-LAT 
collaboration has released a diffused galactic and extra 
galactic 7-ray background ( Abdo et a/.l. 1201 Obi) , however 
no lines were observed yet ( Abdo etall . l2010al ) . 



The high energy neutrinos are being another fron- 
tier for the indirect dark matter searches as the con- 
struction of new experiment IceCube with large volume 
is underway, in which case even Earth could be used 
as a detector to study th e nature of dark matter, e.g. 
( Albuquerque et a/] . 12004), or neutrino emission from 
the passage of Q-balls ( Kusenko and Shoemakeii |2009() . 
Searching for WIMPs in more than one type of experi- 
ments; direct, indirect, and the LHC will be necessary in 
order to stamp the origin of such elusive particles. 
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